BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,028	0	T	C	21.1%	429.1 / 112.6	218	intergenic (‑238/+70)	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 (97/67); new base C (44/0); total (151/67)
Fisher's exact test for biased strand distribution p-value = 3.12e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 220 total aligned reads displayed.

TATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGT  >  NC_000913/3705929‑3706123

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

   aCAGATTGTTCACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑ACAACCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTc                                                                                                <  2:1763744/100‑1 (MQ=255)

    cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                                                                             <  2:851077/100‑1 (MQ=255)

     aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑Gt                                                                                            >  2:1390612/1‑100 (MQ=255)

     aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑Gt                                                                                            <  2:1661503/100‑1 (MQ=255)

      gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTa                                                                                           >  1:2167492/1‑100 (MQ=255)

         cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGc                                                                                        >  1:1967971/6‑100 (MQ=18)

          gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGcc                                                                                       >  2:2100254/1‑100 (MQ=255)

          ataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGcc                                                                                       >  2:1094990/5‑100 (MQ=18)

            tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTg                                                                                     <  1:1113944/100‑1 (MQ=255)

             gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                                                                    >  2:31022/1‑100 (MQ=255)

                  aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                                                               >  1:2023699/1‑100 (MQ=255)

                   gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                                                              >  1:1612443/1‑100 (MQ=255)

                    tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                             >  2:2026162/1‑100 (MQ=255)

                     gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                            <  2:333444/100‑1 (MQ=255)

                       gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                          <  1:1467010/100‑1 (MQ=255)

                        cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                         >  1:387968/1‑100 (MQ=255)

                             tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCa                                                                    >  1:526289/6‑100 (MQ=21)

                             tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                    >  1:868749/6‑100 (MQ=18)

                             tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  1:1755779/6‑100 (MQ=21)

                              gtgggtcGGATAAGGCGTTCACGCCGCATCCGGCCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  1:121162/7‑100 (MQ=18)

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

                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  1:2374983/5‑100 (MQ=18)

                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  2:2076139/5‑100 (MQ=18)

                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  1:2683839/5‑100 (MQ=18)

                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  1:2496886/5‑100 (MQ=18)

                              gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  2:822404/5‑100 (MQ=18)

                               taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                  >  2:2386912/4‑100 (MQ=18)

                               taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                  >  1:1320791/4‑100 (MQ=17)

                               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AGAATCGTGCAAATTCAATACTTTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                  <  2:1469505/100‑1 (MQ=255)

                               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                  >  1:1071859/1‑100 (MQ=255)

                               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                  <  1:2439469/100‑1 (MQ=255)

                                aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                 >  2:576354/3‑100 (MQ=17)

                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                 >  2:688363/3‑100 (MQ=25)

                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                 >  2:2572293/3‑100 (MQ=21)

                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                 >  1:1101506/3‑100 (MQ=17)

                                aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                 >  1:2510372/3‑100 (MQ=17)

                                aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GC‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                 >  2:490034/3‑100 (MQ=17)

                                 ggtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                >  1:2277393/11‑100 (MQ=21)

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

                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                >  2:1438086/2‑100 (MQ=18)

                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                >  1:1668125/2‑100 (MQ=21)

                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                >  1:921400/2‑100 (MQ=18)

                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                                >  1:1782546/2‑100 (MQ=18)

                                  gCCGGATGCGGCGTTAACGCGTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                               >  2:1645094/1‑100 (MQ=255)

                                  gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                                                               >  1:1504541/1‑99 (MQ=18)

                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                               >  2:1994576/1‑99 (MQ=17)

                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                               >  1:347873/1‑99 (MQ=21)

                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgactgt                                                              >  1:2604508/1‑95 (MQ=21)

                                   ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                             >  2:1208603/9‑100 (MQ=18)

                                   ccggataaggcgTTCACGCCGCATCCGGCCTAC‑TAAATCGTGCTAATTCAATATATTGT‑‑G‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                              >  2:1750837/9‑98 (MQ=18)

                                   ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                              >  2:738360/9‑98 (MQ=21)

                                   ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                              >  1:358023/9‑98 (MQ=17)

                                   ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                              >  1:2233213/9‑98 (MQ=17)

                                    cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgc                                                             >  2:1039662/8‑97 (MQ=17)

                                    cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                             >  1:187123/1‑100 (MQ=255)

                                     ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                           >  1:1159360/7‑98 (MQ=21)

                                     ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                            >  2:1547137/7‑100 (MQ=18)

                                     ggataaggcgTTCACGCCTTATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                             >  2:37474/7‑96 (MQ=18)

                                     ggataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GA‑GA‑TGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                                                            >  1:2557327/7‑100 (MQ=14)

                                     ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                            >  1:1002661/7‑100 (MQ=21)

                                      gataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGTCCTGATAAGCGTAACGCATAAGGCaaa                                                              >  1:49881/6‑95 (MQ=17)

                                      gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                           <  2:434665/100‑1 (MQ=255)

                                       aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                          >  2:1138951/1‑100 (MQ=255)

                                        taaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcgcc                                                        >  1:226536/4‑95 (MQ=21)

                                        taaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GA‑GA‑TGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCATTtttgc                                                         >  1:559345/4‑100 (MQ=14)

                                          cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                       <  1:1046350/100‑1 (MQ=255)

                                          cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                       >  2:185546/1‑100 (MQ=255)

                                           ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                      >  2:2454022/1‑100 (MQ=255)

                                            gcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgttt                                                     >  1:469707/1‑100 (MQ=17)

                                            gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                     >  2:345926/1‑100 (MQ=255)

                                              gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                   <  2:2259897/100‑1 (MQ=255)

                                              gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                   <  1:2008538/100‑1 (MQ=255)

                                              gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                   <  2:882843/100‑1 (MQ=255)

                                              gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                   <  2:2556884/100‑1 (MQ=255)

                                               ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                  <  1:190784/100‑1 (MQ=255)

                                               ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                  <  2:2222702/100‑1 (MQ=255)

                                                 aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAATTTTGCCTTTGCCAt                                                >  1:2655377/1‑100 (MQ=255)

                                                  aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                               <  2:2166999/100‑1 (MQ=255)

                                                  aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                               >  2:2233977/1‑100 (MQ=255)

                                                    gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                             >  2:2477066/1‑100 (MQ=255)

                                                       ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                          <  1:1930475/100‑1 (MQ=255)

                                                           ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                      >  2:2495364/1‑100 (MQ=255)

                                                            cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                     >  2:858496/1‑100 (MQ=255)

                                                             ggCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                    <  2:1129491/100‑1 (MQ=255)

                                                             ggCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTGC‑‑GTAGGCCTGATAAGCAGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                    >  2:1885296/1‑100 (MQ=255)

                                                                        aTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                          <  1:732565/100‑1 (MQ=255)

                                                                        aTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                          >  2:2061328/1‑100 (MQ=255)

                                                                          cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                        >  1:2100191/1‑100 (MQ=255)

                                                                            tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGTCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                      <  2:96074/100‑1 (MQ=255)

                                                                             gCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                     >  1:583433/1‑100 (MQ=255)

                                                                                 aatCAATACATTGC‑‑A‑GA‑CT‑TTTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAGCAGTCTTGTATGGCTTTTAATTATTCGATac                 <  2:2528288/98‑2 (MQ=37)

                                                                                    cAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga              <  2:1071851/100‑1 (MQ=255)

                                                                                    cAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga              >  1:1439733/1‑100 (MQ=255)

                                                                                    cAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga              <  1:1185751/100‑1 (MQ=255)

                                                                                             gC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa     >  1:1344104/1‑100 (MQ=255)

                                                                                              c‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa    <  1:21213/100‑1 (MQ=255)

                                                                                                 a‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg   <  2:1339450/100‑1 (MQ=255)

                                                                                                 a‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg   >  2:1375052/1‑100 (MQ=255)

                                                                                                   gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt  >  1:700684/1‑100 (MQ=255)

TATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGT  >  NC_000913/3705929‑3706123

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