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	38.4%	428.4 / inf	276	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 (88/79); new base T (104/0); total (197/79)
Fisher's exact test for biased strand distribution p-value = 2.57e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 283 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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

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

       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  2:510925/8‑100 (MQ=18)

       tgtgTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                              >  1:2206763/1‑100 (MQ=18)

        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑T‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                             >  2:2345972/7‑100 (MQ=17)

         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  1:725937/100‑1 (MQ=255)

         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  2:1826380/1‑100 (MQ=255)

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

         cgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑‑G‑GG‑CT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  1:2512424/6‑100 (MQ=18)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  1:2099995/1‑100 (MQ=255)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  1:37113/1‑100 (MQ=255)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:620308/5‑100 (MQ=17)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑T‑‑‑T‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:780444/5‑100 (MQ=17)

           taggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:2269847/4‑100 (MQ=18)

           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑TC‑T‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:1472395/4‑100 (MQ=18)

           taggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GA‑GA‑T‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:409248/4‑100 (MQ=11)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑‑‑T‑‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:200442/4‑100 (MQ=18)

           taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          >  2:2439937/4‑100 (MQ=18)

            atgtcGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         <  1:2685208/96‑1 (MQ=255)

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

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

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

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

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1661117/3‑100 (MQ=21)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑G‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                                                        >  2:2281968/3‑99 (MQ=18)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑‑‑T‑‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  2:1160865/3‑99 (MQ=18)

             tGCCGGATGCGGCGTTAACGCCTTATCTGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        >  2:2147614/1‑100 (MQ=255)

             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:2650298/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:35164/2‑100 (MQ=21)

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

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

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1549417/2‑100 (MQ=18)

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

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑C‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                                                        >  1:2664858/2‑100 (MQ=18)

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑A‑AT‑T‑‑‑C‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                                                       >  1:1821118/1‑98 (MQ=18)

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:2633705/1‑99 (MQ=17)

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑A‑AT‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:3137849/1‑100 (MQ=21)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:2412234/1‑99 (MQ=17)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:1724615/1‑99 (MQ=18)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑C‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                                                                                       >  1:3089461/1‑99 (MQ=18)

               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                     >  1:3198703/9‑100 (MQ=18)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑AT‑C‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                   >  2:1129459/7‑98 (MQ=18)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:367496/7‑100 (MQ=18)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:2919671/7‑100 (MQ=18)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:2153548/7‑100 (MQ=17)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑C‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                    >  1:2827292/7‑96 (MQ=21)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑C‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                    >  1:2731792/7‑96 (MQ=18)

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

                 ggataaggcgTTCACGCCGCATCCGGCATACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:1138195/7‑100 (MQ=18)

                 ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                    >  1:1326323/7‑100 (MQ=18)

                  gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                                                                                   >  1:889198/6‑100 (MQ=18)

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

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

                    taaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑G‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactacgt                                                                                 >  1:3191400/4‑93 (MQ=25)

                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  2:618983/1‑100 (MQ=255)

                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               <  1:53518/100‑1 (MQ=255)

                      aggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑T‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCATTtatgcgtt                                                                              >  1:1497246/2‑97 (MQ=25)

                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:212371/100‑1 (MQ=255)

                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:1448762/100‑1 (MQ=255)

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

                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:47441/100‑1 (MQ=255)

                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  2:812876/1‑100 (MQ=255)

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

                                              aC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTCGCCATCAGTCTTGTATGGCttt                                                       <  2:1896388/100‑1 (MQ=255)

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

                                                 aaaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     <  1:2301231/100‑1 (MQ=255)

                                                 aaaaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     >  2:738477/1‑100 (MQ=255)

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

                                                      cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                <  2:3179306/100‑1 (MQ=255)

                                                      cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                >  2:3123350/1‑100 (MQ=255)

                                                       gTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  2:1507270/1‑100 (MQ=255)

                                                        tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:3046452/1‑100 (MQ=255)

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

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

                                                          cAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:2022323/1‑100 (MQ=255)

                                                          cAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:1573585/1‑100 (MQ=255)

                                                            aaTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  2:1202125/100‑1 (MQ=255)

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

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

                                                                     cATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 >  1:1396724/1‑100 (MQ=255)

                                                                         gC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTGTGGCTTTTAATTATTCGATAGAGCCGTTTTc                               >  2:2832350/1‑98 (MQ=255)

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

                                                                                a‑aT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         >  2:1590772/1‑100 (MQ=255)

                                                                                  aT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        <  2:725961/100‑1 (MQ=255)

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

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

                                                                                                tAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                 >  2:2094993/1‑100 (MQ=255)

                                                                                                 aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  2:2766287/100‑1 (MQ=255)

                                                                                                  ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  2:2452811/100‑1 (MQ=255)

                                                                                                  ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:2963702/100‑1 (MQ=255)

                                                                                                   gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              <  2:569575/100‑1 (MQ=255)

                                                                                                    catGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGGTTTCCTa              <  2:2720194/97‑1 (MQ=255)

                                                                                                      tGATAAGCGGAGCGCATGAGGCATTTTTGCCATGGCCATCAGTCTGGTATGGCTTTTAATCATGCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:1119886/100‑1 (MQ=255)

                                                                                                      tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:1678994/100‑1 (MQ=255)

                                                                                                      tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:1396710/100‑1 (MQ=255)

                                                                                                      tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  1:3251954/1‑100 (MQ=255)

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

                                                                                                          aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  1:738485/100‑1 (MQ=255)

                                                                                                              ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   >  2:2550115/1‑100 (MQ=255)

                                                                                                               gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  >  1:1974523/1‑100 (MQ=255)

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑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: `‑`