BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,048	0	G	T	39.1%	307.7 / 255.8	208	intergenic (‑258/+50)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base G (58/65); new base T (77/2); total (141/67)
Fisher's exact test for biased strand distribution p-value = 7.68e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 216 total aligned reads displayed.

GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705950‑3706147

gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGa                                                                                                    >  2:1910813/1‑100 (MQ=255)

  ncTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:2018889/2‑100 (MQ=255)

  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  2:47467/100‑1 (MQ=255)

  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:804200/1‑100 (MQ=255)

   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  2:311916/100‑1 (MQ=255)

   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:2441741/1‑100 (MQ=255)

   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:1461731/1‑100 (MQ=255)

   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  1:1122282/100‑1 (MQ=255)

    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                                <  2:582409/100‑1 (MQ=255)

    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTCC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                                <  1:1510908/100‑1 (MQ=255)

      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑C‑‑‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2353187/8‑100 (MQ=17)

       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  1:446705/1‑100 (MQ=255)

        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑‑‑C‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:34661/6‑100 (MQ=18)

        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑C‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:1101117/6‑100 (MQ=18)

          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑G‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1399973/4‑100 (MQ=18)

          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑A‑‑‑C‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:2289317/4‑100 (MQ=18)

          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:1565302/4‑100 (MQ=21)

          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑‑‑T‑‑‑‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:281462/4‑100 (MQ=18)

           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑C‑‑‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1799599/3‑100 (MQ=17)

           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:1255125/3‑100 (MQ=17)

           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑T‑‑‑C‑‑‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                                                         >  2:2975724/3‑100 (MQ=18)

           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:856168/3‑100 (MQ=18)

           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑G‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                                                        >  2:2927278/3‑99 (MQ=21)

           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑C‑‑‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:37439/3‑100 (MQ=21)

           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:609550/3‑100 (MQ=17)

           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGC‑‑G‑GA‑GA‑A‑‑‑A‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  1:914955/3‑99 (MQ=18)

            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑C‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2546412/2‑100 (MQ=18)

            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑A‑AT‑T‑‑‑CAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  2:494779/2‑100 (MQ=21)

            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑A‑‑‑T‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:944018/2‑100 (MQ=21)

            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑C‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1886256/2‑100 (MQ=21)

            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑C‑‑‑T‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2050703/2‑100 (MQ=18)

            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑A‑‑‑C‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1840588/2‑100 (MQ=17)

            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                       >  1:1178203/2‑100 (MQ=21)

              ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                      >  1:2484979/9‑100 (MQ=18)

               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑C‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  1:2384252/8‑97 (MQ=21)

               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑C‑‑‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  1:358444/8‑97 (MQ=17)

               cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                                                                                    >  2:770288/8‑99 (MQ=21)

                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:843012/7‑100 (MQ=17)

                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑T‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:552389/7‑100 (MQ=18)

                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑T‑‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  2:2913916/7‑96 (MQ=17)

                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                   >  1:2757822/7‑98 (MQ=21)

                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑A‑AT‑TCAAC‑‑‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAattt                                                                                    >  1:1569875/6‑100 (MQ=21)

                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑A‑‑‑T‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                                                                                   >  2:1867264/6‑95 (MQ=25)

                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:2539529/100‑1 (MQ=255)

                   taaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑‑G‑GG‑CT‑C‑‑‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                                                 >  1:1130765/4‑100 (MQ=18)

                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑T‑‑‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctttac                                                                                 >  1:1623825/4‑93 (MQ=21)

                    aaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctcgca                                                                                >  2:127389/3‑92 (MQ=18)

                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  1:522768/1‑100 (MQ=255)

                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  2:2583071/1‑100 (MQ=255)

                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  1:983734/1‑100 (MQ=255)

                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  1:1662939/1‑100 (MQ=255)

                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtac                                                                      <  2:1448560/100‑3 (MQ=255)

                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                      <  1:317525/100‑1 (MQ=255)

                                   tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 >  1:3027080/1‑100 (MQ=255)

                                    aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  2:740095/100‑1 (MQ=255)

                                      ccGGCCTAT‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  1:2836218/1‑100 (MQ=255)

                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            <  2:2181151/100‑1 (MQ=255)

                                          ccTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          >  1:1753765/1‑100 (MQ=255)

                                           cTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         >  2:2090293/1‑100 (MQ=255)

                                           cTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         >  2:43908/1‑100 (MQ=255)

                                            tAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                        >  1:588983/1‑100 (MQ=255)

                                              c‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      <  1:1325530/100‑1 (MQ=255)

                                              c‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGAATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      >  2:616922/1‑100 (MQ=255)

                                                gaaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     <  2:2437148/99‑1 (MQ=255)

                                                 aaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  2:2396108/100‑1 (MQ=255)

                                                  aaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑C‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTGTa                                                     <  2:855359/98‑1 (MQ=255)

                                                        gCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:23375/100‑1 (MQ=255)

                                                          aaaTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  1:2852497/100‑1 (MQ=255)

                                                           aaTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCGTCAGGCATTTTTGCCTTTGCCATCAGTATTTTATGGCTTTTAATTATTCGATa                                          >  1:631117/1‑100 (MQ=255)

                                                            aTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                         <  2:1634567/100‑1 (MQ=255)

                                                             ttCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                        <  2:251794/100‑1 (MQ=255)

                                                              tCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:1753750/100‑1 (MQ=255)

                                                               cAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  1:1023507/100‑1 (MQ=255)

                                                               cAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      >  1:1987479/1‑100 (MQ=255)

                                                                  tACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGt                                   >  1:686296/1‑100 (MQ=255)

                                                                       tGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              <  1:210969/100‑1 (MQ=255)

                                                                        gC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  2:1646648/1‑100 (MQ=255)

                                                                        gC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:1046413/1‑100 (MQ=255)

                                                                        gC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:2064893/1‑100 (MQ=255)

                                                                        gC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:2899275/1‑100 (MQ=255)

                                                                        gC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:2499318/100‑1 (MQ=255)

                                                                         c‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  2:2755167/1‑100 (MQ=255)

                                                                         c‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  1:1722543/1‑100 (MQ=255)

                                                                            a‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:279152/100‑1 (MQ=255)

                                                                            a‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  1:1193885/100‑1 (MQ=255)

                                                                              gA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  2:1055707/1‑100 (MQ=255)

                                                                                 aT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        >  1:2707298/1‑100 (MQ=255)

                                                                                  t‑t‑‑‑t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                       >  2:569830/1‑100 (MQ=255)

                                                                                            tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                    <  2:1294885/100‑1 (MQ=255)

                                                                                            tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                    <  2:1293477/100‑1 (MQ=255)

                                                                                                gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  >  2:2712917/1‑100 (MQ=255)

                                                                                                gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:2136179/100‑1 (MQ=255)

                                                                                                   ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  1:1586448/1‑100 (MQ=255)

                                                                                                          tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        <  2:2349340/100‑1 (MQ=255)

                                                                                                              cagAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTATAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     <  2:1359953/97‑1 (MQ=255)

                                                                                                               ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcn   <  2:2018871/100‑2 (MQ=255)

                                                                                                               ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:438258/100‑1 (MQ=255)

                                                                                                               ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  2:804193/100‑1 (MQ=255)

                                                                                                                gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  <  1:2123702/100‑1 (MQ=255)

                                                                                                                gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCCTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  <  1:839631/100‑1 (MQ=255)

GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705950‑3706147

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