D22 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	918	A→G	T213T (ACA→ACG)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	918	0	A	G	91.0%	414.5 / 18.1	178	T213T (ACA→ACG)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base A (15/1); new base G (157/4); total (173/5)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 196 total aligned reads displayed.

GTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAAC‑AACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGA  >  NC_007790/819‑1002

gTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaac‑a                                                                                      >  1:1390916/1‑100 (MQ=255)

  tAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAgac‑gac                                                                                    <  1:949061/100‑1 (MQ=255)

   aaGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAac‑gaca                                                                                   >  1:3275000/1‑100 (MQ=255)

   aaGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAac‑aaca                                                                                   >  1:1053711/1‑100 (MQ=255)

   aaGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAac‑aaca                                                                                   >  1:4830429/1‑100 (MQ=255)

    ggTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACTATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAg                                                                                  >  1:1642547/2‑100 (MQ=255)

       tctggtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAATAC‑GACAGatc                                                                               >  1:4351662/3‑98 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGGTGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACg                                                                               >  1:1967399/1‑100 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTCATTTGATTAATACAAGTGATGAAGATAAGAC‑GACAGatc                                                                               >  1:4204706/1‑98 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaac‑aagcagat                                                                               >  1:3137634/1‑94 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAGAT‑GACAGatc                                                                               >  1:1268927/1‑98 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAACAC‑GACAGatc                                                                               >  1:1828131/1‑98 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑AACAGatc                                                                               >  1:267243/1‑98 (MQ=255)

       tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GAGAGACg                                                                               >  1:2191025/1‑100 (MQ=255)

       tATGGTGGATTACCTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAACAC‑GACAGatc                                                                               >  1:3219173/1‑98 (MQ=255)

       gaTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACg                                                                               >  1:3224719/2‑100 (MQ=255)

       cttggtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaac‑acagatcg                                                                               >  1:530341/3‑93 (MQ=255)

       caTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAG‑GACAGatc                                                                               >  1:3120740/2‑98 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAGGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                              >  1:3457267/1‑100 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGTTGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑GACAGatcg                                                                              >  1:808183/1‑97 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaac‑gaccagatc                                                                              >  1:509112/1‑94 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAATAC‑GACAGatcg                                                                              >  1:195568/1‑97 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAGAC‑GACAGatcg                                                                              >  1:4737057/1‑97 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAac‑gacaagaga                                                                              >  1:2886136/1‑95 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑GACAGatcg                                                                              >  1:3553086/1‑97 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGAACAGatc                                                                               >  1:1874474/1‑98 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAgccca                                                                              >  1:3295375/1‑96 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                              >  1:4058479/1‑100 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                              >  1:4313446/1‑100 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                              >  1:4789539/1‑100 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                              >  1:573722/1‑100 (MQ=255)

        aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑AACAGatcg                                                                              >  1:3471113/1‑97 (MQ=255)

         tggtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAgc‑gacaagatcg                                                                             >  1:1683060/1‑94 (MQ=255)

         tggtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACgaa                                                                             >  1:1153999/1‑100 (MQ=255)

          ggtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACgaaa                                                                            >  1:2851411/1‑99 (MQ=255)

            tggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACTGACgaagat                                                                          >  1:2326862/1‑99 (MQ=255)

            acgATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAAC‑AACAGACgaagat                                                                          >  1:2957162/3‑99 (MQ=255)

             ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAGGC‑GACAGatcggaaga                                                                         >  1:5191940/1‑92 (MQ=255)

             ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGatcggaaga                                                                         >  1:150351/1‑92 (MQ=255)

               aTTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGatcggaagagc                                                                       >  1:2902291/1‑90 (MQ=255)

                 tacTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGaaaa                                                                     <  1:2225754/98‑1 (MQ=255)

                 tacTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGaaaa                                                                     >  1:4891305/3‑100 (MQ=255)

                          aaaagCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTc                                                            >  1:871799/1‑100 (MQ=255)

                             agcATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAAt                                                         >  1:2622231/3‑100 (MQ=255)

                             agcATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGAt                                                         >  1:2344990/3‑100 (MQ=255)

                                   aaaTTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAATTACAGc                                                   <  1:4782792/100‑1 (MQ=255)

                                         tAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTAcgcagatcgg                                             >  1:2037708/1‑91 (MQ=255)

                                          aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTgg                                            >  1:3878897/1‑100 (MQ=255)

                                          aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATATgg                                            >  1:1933850/1‑100 (MQ=255)

                                          aaaCCTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCaaagatc                                            >  1:1016498/1‑95 (MQ=255)

                                            gcctAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaa                                          >  1:4849241/4‑100 (MQ=255)

                                               tagataaTGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATCTGGaaaga                                       >  1:926568/7‑97 (MQ=255)

                                               tAGACNATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:1290492/1‑100 (MQ=255)

                                               tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCaagatcggaaga                                       >  1:3381636/1‑90 (MQ=255)

                                               tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:1190852/1‑100 (MQ=255)

                                               tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:4582547/1‑100 (MQ=255)

                                                cgACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                      >  1:4269807/2‑100 (MQ=255)

                                                aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GGCAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                      >  1:3075960/1‑100 (MQ=255)

                                                aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                      >  1:3502045/1‑100 (MQ=255)

                                                aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                      >  1:3927968/1‑100 (MQ=255)

                                                  aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTAGGAAtagatc                                    >  1:3849257/1‑95 (MQ=255)

                                                  aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATGTGGaagatcgg                                    >  1:4643995/1‑94 (MQ=255)

                                                   tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:4947558/2‑100 (MQ=255)

                                                   tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:4126783/2‑100 (MQ=255)

                                                   tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTCGGaagatcgga                                   >  1:615312/2‑93 (MQ=255)

                                                   ccaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaagatcgga                                   >  1:1615269/3‑93 (MQ=255)

                                                   cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:58136/1‑100 (MQ=255)

                                                   cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:4494282/1‑100 (MQ=255)

                                                    gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:3513393/1‑100 (MQ=255)

                                                    gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:879901/1‑100 (MQ=255)

                                                    gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:1127035/1‑100 (MQ=255)

                                                    gATGCCGAAAATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaagatcggaa                                  >  1:179832/1‑92 (MQ=255)

                                                    aaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaa                                  >  1:3821780/2‑92 (MQ=255)

                                                     ttgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:3999501/5‑100 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:934506/1‑100 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:1197030/1‑100 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:3185821/1‑100 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:471166/1‑100 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAAGTTGGaagatcggaag                                 >  1:2531838/1‑91 (MQ=255)

                                                     aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaag                                 >  1:3766346/1‑91 (MQ=255)

                                                      tgccGAAGGTGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                >  1:1557283/4‑100 (MQ=255)

                                                      tgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                >  1:3586584/4‑100 (MQ=255)

                                                      tgccGAAGATGGAAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                >  1:776867/4‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:356874/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:4436497/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:3639706/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:2979119/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:2231337/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:3806033/3‑100 (MQ=255)

                                                       gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                               >  1:3668742/3‑100 (MQ=255)

                                                             tgtGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                         >  1:4540774/3‑100 (MQ=255)

                                                             actGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                         >  1:496171/3‑100 (MQ=255)

                                                              aTGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAGGCAAAatta                        >  1:1370015/1‑100 (MQ=255)

                                                               tGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattat                       >  1:829016/1‑100 (MQ=255)

                                                                ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                      >  1:1327673/1‑100 (MQ=255)

                                                                ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                      >  1:1823001/1‑100 (MQ=255)

                                                                        gATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaa              >  1:2262278/1‑98 (MQ=255)

                                                                                taaGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaag      >  1:5254427/2‑96 (MQ=255)

                                                                                cAAGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGTCATCAgat      >  1:2399702/1‑98 (MQ=255)

                                                                                    tgatgaAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaaggagt  >  1:4715017/1‑92 (MQ=255)

GTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAAC‑AACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGA  >  NC_007790/819‑1002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`