BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,029	0	T	C	29.7%	307.0 / 152.7	212	intergenic (‑239/+69)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base T (78/59); new base C (58/0); total (150/62)
Fisher's exact test for biased strand distribution p-value = 7.30e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 217 total aligned reads displayed.

ATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATG  >  NC_000913/3705930‑3706127

 taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑ttt                                                                                                   >  1:421559/1‑100 (MQ=255)

 taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑ttt                                                                                                   <  2:488437/100‑1 (MQ=255)

  aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTc                                                                                                  >  2:1407566/1‑100 (MQ=255)

   cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCg                                                                                                 >  1:154900/1‑100 (MQ=255)

    aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGt                                                                                                >  2:1417621/1‑100 (MQ=255)

       ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAgg                                                                                             <  2:509674/100‑1 (MQ=255)

       gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAgg                                                                                             >  1:1916155/7‑100 (MQ=18)

        tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGc                                                                                            <  1:996075/100‑1 (MQ=255)

              cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATa                                                                                      >  1:1630597/1‑100 (MQ=255)

                   tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                                                                 <  1:820293/100‑1 (MQ=255)

                   tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                                                                 <  1:1726720/100‑1 (MQ=255)

                   tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                                                                 >  2:1836634/1‑100 (MQ=255)

                    gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGa                                                                                <  1:2240555/100‑1 (MQ=255)

                     cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAg                                                                               <  1:1779697/100‑1 (MQ=255)

                     cgcTTTGCTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTTC‑‑A‑GA‑AT‑TCTCGTAGGCCTGATAAGcag                                                                                 >  2:1713307/1‑96 (MQ=17)

                      gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                                                              >  1:2349384/1‑100 (MQ=255)

                      gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                                                              <  1:216799/100‑1 (MQ=255)

                       gcttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAgcg                                                                             >  2:2197671/11‑100 (MQ=17)

                        cttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAgcgc                                                                            >  2:965775/10‑100 (MQ=17)

                          tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCAt                                                                          >  1:1261428/8‑100 (MQ=14)

                          tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                                                                          >  2:757035/1‑100 (MQ=255)

                           ttgtaggCCGGATAAGGCGTTCGCGCCGCATCCGGCATAC‑AAAAGTTTGCAAATTCAATAAATTGC‑‑A‑GA‑AT‑GCATGTAGGCCTGATAAGCGAAGCACATc                                                                         >  2:2452989/7‑100 (MQ=12)

                             ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                       <  2:2337828/100‑1 (MQ=255)

                             ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                       <  2:421558/100‑1 (MQ=255)

                             gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                       >  2:1245512/5‑100 (MQ=17)

                             gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAg                                                                       >  2:2422358/5‑100 (MQ=21)

                             gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAg                                                                       >  1:5693/5‑100 (MQ=18)

                              taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                      >  1:2459294/4‑100 (MQ=18)

                              taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                     >  1:356265/4‑100 (MQ=21)

                               aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                     >  2:2019591/3‑100 (MQ=18)

                               aggCCGGATAAGGCGGTCACGCCGCATCCGGCCTAC‑GAAACCGCAAGAATTCAATATATTGC‑‑A‑GG‑AG‑CGGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                     >  2:1635350/3‑100 (MQ=18)

                                ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  2:2830079/2‑100 (MQ=17)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:2578198/2‑100 (MQ=18)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:2647303/2‑100 (MQ=18)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:1055559/2‑100 (MQ=17)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                    >  2:1327292/2‑100 (MQ=21)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                    >  1:2689107/2‑100 (MQ=25)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  2:1028915/2‑100 (MQ=18)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:1217599/2‑100 (MQ=17)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                                    >  2:2094932/2‑100 (MQ=18)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  2:600718/2‑100 (MQ=18)

                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                   >  1:2828885/2‑100 (MQ=21)

                                ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACTTTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGAGCATCAGGc                                                                     >  2:2493784/2‑99 (MQ=17)

                                 gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                   >  1:310663/1‑100 (MQ=255)

                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                  >  1:2079449/1‑100 (MQ=18)

                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  1:806694/1‑99 (MQ=17)

                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  1:2189986/1‑99 (MQ=17)

                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                   >  1:1572285/1‑99 (MQ=21)

                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCag                                                                   >  2:1508364/1‑99 (MQ=18)

                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                                                                  >  2:2659404/1‑97 (MQ=21)

                                  ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                  >  1:2310095/9‑98 (MQ=25)

                                  ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                  >  1:1874648/1‑100 (MQ=255)

                                  ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  1:2224785/1‑98 (MQ=17)

                                   cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                 >  2:720531/8‑97 (MQ=21)

                                   cgatctaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                                                                 <  1:1146744/93‑4 (MQ=21)

                                   cgatctaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                                                                 <  1:2597716/93‑4 (MQ=25)

                                    ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                >  1:2001385/7‑100 (MQ=18)

                                    ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                >  1:2747323/7‑100 (MQ=17)

                                    ggataaggcgTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                               >  1:2136497/7‑98 (MQ=18)

                                    ggataaggcgTTCACGCCTTATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                >  2:2629190/7‑100 (MQ=21)

                                    ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCACGTAGGCCTGATAATCGTAGCTCATCAAGCAgttt                                                                >  1:1106356/7‑100 (MQ=17)

                                    ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                >  1:1298466/7‑96 (MQ=18)

                                    ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                               >  2:2271105/7‑98 (MQ=21)

                                    ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TCATGTACGCCTGATACGCGTAGCGCATCAGGCAtttat                                                               >  2:1952617/7‑98 (MQ=21)

                                     gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                                                               >  1:1641059/6‑95 (MQ=21)

                                      aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                              <  2:2841885/100‑1 (MQ=255)

                                       taaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAattctgt                                                             >  2:1276733/4‑96 (MQ=21)

                                       tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                             >  1:610356/1‑100 (MQ=255)

                                        aaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatcttacg                                                            <  1:1327297/98‑9 (MQ=37)

                                         aggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactagatcg                                                           >  2:2597697/2‑91 (MQ=25)

                                          ggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgtt                                                          <  2:2578180/100‑1 (MQ=18)

                                          ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                          >  1:828002/1‑100 (MQ=255)

                                             gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                       <  2:1876493/100‑1 (MQ=255)

                                             gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                       <  2:1221602/100‑1 (MQ=255)

                                              ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                      >  1:134289/1‑100 (MQ=255)

                                              ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                      <  2:1460154/100‑1 (MQ=255)

                                              ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                      <  2:123995/100‑1 (MQ=255)

                                              ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                      <  2:1596941/100‑1 (MQ=255)

                                                aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                     <  1:1407571/99‑1 (MQ=255)

                                                aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                    >  1:1967762/1‑100 (MQ=255)

                                                 aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                   <  1:704248/100‑1 (MQ=255)

                                                     cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                               <  2:1127792/100‑1 (MQ=255)

                                                     cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                               <  2:769853/100‑1 (MQ=255)

                                                      ttATCCGGCCTAC‑AAGATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                              <  1:2846753/100‑1 (MQ=255)

                                                           cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                         <  2:2057089/100‑1 (MQ=255)

                                                             gCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                       <  2:1708506/100‑1 (MQ=255)

                                                                  c‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                  >  2:333693/1‑100 (MQ=255)

                                                                    aaaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                 >  1:2436113/1‑100 (MQ=255)

                                                                      aaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                               <  2:377670/100‑1 (MQ=255)

                                                                         cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                            >  2:1085040/1‑100 (MQ=255)

                                                                           tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                          <  2:1808652/100‑1 (MQ=255)

                                                                              aaaTTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                       <  1:1843325/100‑1 (MQ=255)

                                                                               aaTTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                      <  1:97915/100‑1 (MQ=255)

                                                                                  tCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                   <  2:271996/100‑1 (MQ=255)

                                                                                  tCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                   <  2:825323/100‑1 (MQ=255)

                                                                                   cAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                  <  1:1010727/100‑1 (MQ=255)

                                                                                   cAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                  <  2:2669874/100‑1 (MQ=255)

                                                                                                a‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg       >  1:1935516/1‑100 (MQ=255)

                                                                                                  gA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt      >  2:1771653/1‑100 (MQ=255)

                                                                                                    naT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTGGCCTTTGCCATCAGTCTTGAATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg     >  2:1580879/2‑100 (MQ=255)

ATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATG  >  NC_000913/3705930‑3706127

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: `‑`