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	26.3%	513.6 / 204.1	324	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 (116/100); new base C (76/1); total (222/102)
Fisher's exact test for biased strand distribution p-value = 6.14e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 331 total aligned reads displayed.

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGT  >  NC_000913/3705931‑3706128

taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑At‑t‑‑‑tt                                                                                                    <  2:922424/100‑1 (MQ=255)

 aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTc                                                                                                   <  2:15635/100‑1 (MQ=255)

  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCg                                                                                                  <  1:1017674/100‑1 (MQ=255)

  cAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCGGACTAC‑AAAGTCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCg                                                                                                  <  2:3840933/100‑1 (MQ=255)

   aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGt                                                                                                 >  2:3597139/1‑100 (MQ=255)

      ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAgg                                                                                              <  2:2897720/100‑1 (MQ=255)

      ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAgg                                                                                              >  2:2446772/1‑100 (MQ=255)

            aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGAt                                                                                        <  2:2758982/100‑1 (MQ=255)

             cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATa                                                                                       <  1:1209518/100‑1 (MQ=255)

                aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGc                                                                                    >  2:1945300/1‑100 (MQ=255)

                  tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCg                                                                                   <  2:663491/99‑1 (MQ=255)

                    cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAg                                                                                >  1:3508771/1‑100 (MQ=255)

                      cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                                                                              <  1:949760/100‑1 (MQ=255)

                      cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                                                                              >  1:3380532/1‑100 (MQ=255)

                       tccgatctggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑AGG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAgcgc                                                                             <  2:3005078/91‑1 (MQ=18)

                          gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc                                                                          >  1:2772259/1‑100 (MQ=255)

                           tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                         <  1:2322353/100‑1 (MQ=255)

                           tgtATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                         >  2:3547255/3‑100 (MQ=255)

                            ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                        <  2:1476382/100‑1 (MQ=255)

                            ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                        <  1:1381135/100‑1 (MQ=255)

                            gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑AGA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                        >  1:584944/5‑100 (MQ=18)

                            gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑AGG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                        >  1:2372414/5‑100 (MQ=18)

                            gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAACATGCAAATTCAATAGATTGC‑AGAGAT‑T‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                        >  1:1999707/5‑100 (MQ=17)

                             taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑AGA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                       >  2:3441225/4‑100 (MQ=18)

                             tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                                       <  2:2006929/100‑1 (MQ=255)

                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:1768886/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑AAT‑CT‑T‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                                                                      >  1:2736501/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑AAT‑CT‑T‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                                                                      >  2:3266462/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCATGCAAATTCAATATATTGC‑AGA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                      >  1:3042469/3‑100 (MQ=17)

                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑ATG‑AG‑C‑‑‑CACGTATGCCTGATAAGCGTAGCGCATCa                                                                         >  1:198084/3‑97 (MQ=17)

                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                      >  1:1351576/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑AGA‑GA‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                      >  1:1300363/3‑100 (MQ=17)

                              aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                      >  1:1572263/1‑100 (MQ=255)

                               ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑AGA‑AA‑T‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  2:1903404/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑AGA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:471890/2‑100 (MQ=17)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑AGG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCagatcgga                                                                     >  1:3005102/2‑93 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑GGG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:2157786/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑AGT‑TT‑A‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                     >  2:452755/2‑100 (MQ=21)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  2:643599/2‑100 (MQ=17)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGC‑A‑A‑TA‑T‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                                    >  2:860427/2‑98 (MQ=21)

                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑AGA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                    >  1:2074816/1‑99 (MQ=17)

                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑AGG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                    >  1:745626/1‑99 (MQ=18)

                                 ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                    >  1:3341638/9‑99 (MQ=21)

                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                   >  1:3246393/9‑98 (MQ=18)

                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑AGG‑AG‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                  >  2:123381/9‑100 (MQ=21)

                                 ccggataaggcgTTCACGCCGCATCCGGCCGAC‑AAAAACCTGCAAATTCAATATATTGC‑AGG‑AT‑A‑‑‑CGTGTAGGCCTGATAAGCGTAGCGAATCAGGCaat                                                                   >  2:2550037/9‑98 (MQ=17)

                                  cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑AGA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                                                                  >  2:3734067/8‑97 (MQ=17)

                                  cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑ATG‑AG‑C‑‑‑CACGTATGCCTGATAATCGTAGAGCATCAGACaa                                                                    >  1:668337/8‑97 (MQ=17)

                                  cggataaggcgTTCACGCCACATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑AGG‑AT‑A‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                  >  1:1471993/8‑97 (MQ=18)

                                   ngataaggcgTTTATGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATACATTGCAATT‑CA‑A‑‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                   >  2:1643385/7‑97 (MQ=18)

                                   ggataagtcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  1:2041409/9‑100 (MQ=17)

                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑AGA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  2:51071/7‑100 (MQ=18)

                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑AGA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  1:3640349/7‑100 (MQ=18)

                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑AGA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  2:2240005/7‑100 (MQ=18)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑GGG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  1:3432884/7‑100 (MQ=21)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑AGT‑TT‑A‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                 >  2:2403686/7‑96 (MQ=25)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  1:2305955/7‑100 (MQ=17)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑AGG‑AA‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                                 >  2:2127903/7‑100 (MQ=18)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑AGG‑AA‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                                 >  1:3657859/7‑100 (MQ=18)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑AGG‑GA‑C‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                 >  1:2462539/7‑96 (MQ=21)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑AGG‑GA‑C‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                 >  1:613220/7‑96 (MQ=21)

                                   ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGC‑ACA‑GC‑T‑‑‑AGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                                                                 >  2:333001/7‑96 (MQ=37)

                                    gataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑ATG‑AG‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                >  1:2221857/6‑100 (MQ=17)

                                    gataaggcgTTCACGCCGCATCCGGCCTAC‑AAACGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                                                                >  2:2337833/6‑99 (MQ=18)

                                    gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                <  1:828571/100‑1 (MQ=255)

                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATTGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGTATTtttg                                                               <  2:2044357/100‑1 (MQ=255)

                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                               <  1:1150871/100‑1 (MQ=255)

                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                               <  2:2468697/100‑1 (MQ=255)

                                        cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                            <  1:958040/100‑1 (MQ=255)

                                        aggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑AGG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTttggcgt                                                            <  1:298377/99‑3 (MQ=18)

                                         ggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑AAT‑CT‑T‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagcttacatt                                                           <  1:2425370/100‑11 (MQ=25)

                                               aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                     >  2:760196/1‑100 (MQ=255)

                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                    <  2:1864006/100‑1 (MQ=255)

                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                    <  2:2710231/100‑1 (MQ=255)

                                                 ngCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                   >  2:1392847/2‑100 (MQ=255)

                                                   ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCTTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                  <  1:166732/99‑1 (MQ=255)

                                                   ccTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGC‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                 >  2:317812/1‑100 (MQ=255)

                                                       aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                             <  1:3532228/100‑1 (MQ=255)

                                                        tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                            >  1:497173/1‑100 (MQ=255)

                                                          cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                          >  2:1154046/1‑100 (MQ=255)

                                                             ccTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                       <  2:1022980/100‑1 (MQ=255)

                                                                aC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                    >  1:438354/1‑100 (MQ=255)

                                                                     aaTCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                <  1:3578929/100‑1 (MQ=255)

                                                                       tCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                              <  2:3676762/100‑1 (MQ=255)

                                                                          tGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                           <  2:2800608/100‑1 (MQ=255)

                                                                           gCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                          <  1:2511543/100‑1 (MQ=255)

                                                                           gCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                          <  1:2826393/100‑1 (MQ=255)

                                                                              aaTTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                       <  2:3019877/100‑1 (MQ=255)

                                                                               aTTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                      <  2:1990638/100‑1 (MQ=255)

                                                                               aTTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                      <  2:2630628/100‑1 (MQ=255)

                                                                                 tCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                    >  1:1822352/1‑100 (MQ=255)

                                                                                   aaTACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                  <  1:384416/100‑1 (MQ=255)

                                                                                        aTTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt             >  1:1548014/1‑100 (MQ=255)

                                                                                         ttGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTCTTAATTATTCGATAGAGACGGTTTc            >  2:219426/1‑100 (MQ=255)

                                                                                           gC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATCGAGACGTTTTcca          <  1:2631100/100‑3 (MQ=255)

                                                                                           gC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa          >  2:1079686/1‑100 (MQ=255)

                                                                                           gC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa          >  1:1403765/1‑100 (MQ=255)

                                                                                           gC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa          >  1:855937/1‑100 (MQ=255)

                                                                                               gA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt       <  1:2197738/100‑1 (MQ=255)

                                                                                                  aT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa     >  2:580635/1‑100 (MQ=255)

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑AGA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGT  >  NC_000913/3705931‑3706128

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