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	24.4%	450.4 / 147.6	253	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 (98/79); new base C (55/2); total (171/82)
Fisher's exact test for biased strand distribution p-value = 5.89e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 262 total aligned reads displayed.

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGT  >  NC_000913/3705931‑3706128

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

taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At‑t‑‑‑tt                                                                                                      >  1:1665000/1‑100 (MQ=255)

taCAGATTGCTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At‑t‑‑‑tt                                                                                                      >  1:1395319/1‑100 (MQ=255)

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

  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑g                                                                                                  <  2:2575269/100‑1 (MQ=255)

  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑g                                                                                                  <  1:1514521/100‑1 (MQ=255)

     aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAg                                                                                               <  1:1022406/100‑1 (MQ=255)

        gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGcc                                                                                            <  1:2213964/100‑1 (MQ=255)

        gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGcc                                                                                            >  2:812691/1‑100 (MQ=255)

        gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGcc                                                                                            >  1:1525502/1‑100 (MQ=255)

                 gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                                                                   >  2:1306295/1‑100 (MQ=255)

                  tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                  <  1:225885/100‑1 (MQ=255)

                  tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                  <  1:771746/100‑1 (MQ=255)

                    cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                <  2:730320/100‑1 (MQ=255)

                      cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                              <  2:2717074/100‑1 (MQ=255)

                      cTTTGTTTATGCCGGAGGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                              <  2:626922/100‑1 (MQ=255)

                        ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCa                                                                            >  2:1747108/9‑100 (MQ=18)

                         tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGG‑GA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCAt                                                                           >  2:543247/8‑100 (MQ=18)

                         cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                           >  2:382255/8‑100 (MQ=17)

                          ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                         >  2:2105223/7‑100 (MQ=18)

                           tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGC‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                         >  1:1220173/1‑100 (MQ=255)

                            gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                         >  2:2697775/5‑99 (MQ=17)

                             tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                       >  1:730472/1‑100 (MQ=255)

                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:2140386/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                      >  1:2302770/3‑100 (MQ=18)

                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTGCGGAGA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                     >  2:2051672/3‑99 (MQ=18)

                              aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                      >  2:2460252/3‑100 (MQ=18)

                              aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                      >  2:512177/1‑100 (MQ=255)

                               tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                     <  2:2130947/100‑1 (MQ=255)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  2:1443462/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCa                                                                     >  2:167696/2‑100 (MQ=17)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AAT‑CT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCANCAGGCa                                                                     >  2:2587315/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑AGG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                     >  1:894326/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑AGG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                     >  1:787710/2‑100 (MQ=18)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑AGAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:2500614/2‑100 (MQ=17)

                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑AGG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                    >  2:2547557/2‑100 (MQ=21)

                               ggCCGGATAAGGCGATTACGCCGCATCCTGCCTAC‑AAAAGCATGCAAATTCCATATATTAC‑‑AGA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  2:1247841/2‑100 (MQ=14)

                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                   >  1:1083590/1‑100 (MQ=18)

                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑A‑A‑TT‑C‑‑‑AACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                     >  2:1663604/1‑100 (MQ=21)

                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                    >  1:1802046/1‑99 (MQ=17)

                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                    >  1:1270798/1‑99 (MQ=17)

                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑AGG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                    >  2:2205404/1‑99 (MQ=18)

                                 ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                  >  1:1707809/9‑100 (MQ=18)

                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑TCG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                   >  1:1428842/9‑98 (MQ=21)

                                 ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                   >  1:2412885/1‑100 (MQ=255)

                                 ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                   >  2:1078306/1‑100 (MQ=255)

                                 ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                   >  2:2343011/1‑98 (MQ=18)

                                  cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑AGA‑GT‑T‑‑‑CAC‑‑GTAGGACTGATAAGCGTAGCGCATCAGGCAgtt                                                                  >  1:5191/8‑100 (MQ=17)

                                  cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                                                                  >  2:153092/8‑97 (MQ=18)

                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                >  2:1472089/7‑98 (MQ=21)

                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                >  2:2236084/7‑98 (MQ=21)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                 >  1:2088474/7‑100 (MQ=17)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑AGAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                 >  2:2415926/7‑100 (MQ=17)

                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑AGG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGTGCATCAGGCAtttat                                                                >  2:1862266/7‑98 (MQ=21)

                                    gataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑ATG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                >  1:1142236/6‑100 (MQ=17)

                                    gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                                                                >  1:483435/6‑99 (MQ=18)

                                    gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                <  2:369744/100‑1 (MQ=255)

                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                               <  1:1103607/100‑1 (MQ=255)

                                      taaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                              >  1:2337810/4‑100 (MQ=18)

                                      tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                              >  2:2099407/1‑100 (MQ=255)

                                         ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                           >  2:455464/1‑100 (MQ=255)

                                           cgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatttggcgttg                                                          >  2:1392888/1‑93 (MQ=25)

                                           cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                         <  1:2667617/100‑1 (MQ=255)

                                            gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                        <  2:291204/100‑1 (MQ=255)

                                              tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                      >  1:683596/1‑100 (MQ=255)

                                               aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                     >  1:1721843/1‑100 (MQ=255)

                                                  gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                  <  2:1594643/100‑1 (MQ=255)

                                                   ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                 >  2:1094451/1‑100 (MQ=255)

                                                    cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                >  1:1997671/1‑100 (MQ=255)

                                                     ttATCCGTCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                               <  1:1563710/100‑1 (MQ=255)

                                                        tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                            <  1:2875642/100‑1 (MQ=255)

                                                        tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                            >  2:2697195/1‑100 (MQ=255)

                                                          cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTAt                                          >  1:1883816/1‑100 (MQ=18)

                                                              cTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                      >  2:1993701/1‑100 (MQ=255)

                                                              cTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                      >  1:436188/1‑100 (MQ=255)

                                                               tAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                     >  1:1703074/1‑100 (MQ=255)

                                                                 c‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                   <  2:2784966/100‑1 (MQ=255)

                                                                   aaaaTCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                  <  1:500813/100‑1 (MQ=255)

                                                                    aaaTCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                 >  1:341174/1‑100 (MQ=255)

                                                                        cGTGCAAATTCAATACATTGC‑‑AGA‑AC‑T‑‑‑TTC‑‑GTAGGCCTGCTAAGCGGAGAGCATAAGGCATTTTTGCCTTTGCCATCAGCCTTGTATGTTTTTTAattat                             <  1:963589/100‑1 (MQ=25)

                                                                         gTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                            <  2:1213455/100‑1 (MQ=255)

                                                                          tGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                           <  1:1757680/100‑1 (MQ=255)

                                                                          tGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                           <  1:1840272/100‑1 (MQ=255)

                                                                           gCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                          >  2:887178/1‑100 (MQ=255)

                                                                             aaaTTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                        <  1:819096/100‑1 (MQ=255)

                                                                             aaaTTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                        <  1:384412/100‑1 (MQ=255)

                                                                              aaTTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                       >  2:642843/1‑100 (MQ=255)

                                                                                 tCAATACATTGC‑‑AGA‑TT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                    <  2:2857606/100‑1 (MQ=255)

                                                                                 tCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                    >  2:2286430/1‑100 (MQ=255)

                                                                                  cAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                   <  1:1463185/100‑1 (MQ=255)

                                                                                     tACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGt                <  1:1119223/100‑1 (MQ=255)

                                                                                           gC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa          >  2:1159336/1‑100 (MQ=255)

                                                                                            c‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa         >  1:2368988/1‑100 (MQ=255)

                                                                                            c‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATATTTGCCTTTGCCACCAGACTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTTGaa         >  2:2760226/1‑100 (MQ=38)

                                                                                               aGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg        <  1:375629/100‑1 (MQ=255)

                                                                                                 a‑aT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg      >  1:2603784/1‑100 (MQ=255)

                                                                                                 a‑aT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg      <  1:455459/100‑1 (MQ=255)

                                                                                                   aT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa     >  2:19558/1‑100 (MQ=255)

                                                                                                      t‑‑‑ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTGTGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg   <  2:858903/100‑1 (MQ=255)

                                                                                                          ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt  >  2:1833383/1‑100 (MQ=255)

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGT  >  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: `‑`