Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	129,360	G→T	100%	intergenic (+24/+47)	lpd → / ← yacH	dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	129,360	0	G	T	99.0%	324.4 / ‑4.0	100	intergenic (+24/+47)	lpd/yacH	dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (52/47);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATACTG  >  NC_000913/129263‑129456                                                                                                  |                                                                                                 cTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCa                                                                                                >  1:1148125/1‑100 (MQ=255)  tgtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               >  1:1516926/1‑100 (MQ=255)  tgtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               >  2:1026169/1‑100 (MQ=255)  tgtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAAAAGAAGTAAGTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               >  2:1229203/1‑100 (MQ=255)   gtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               <  2:2421741/99‑1 (MQ=255)    tgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAAtt                                                                                             >  2:59991/1‑100 (MQ=255)     gggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTa                                                                                            >  1:1246662/1‑100 (MQ=255)     gggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTa                                                                                            <  1:2219821/100‑1 (MQ=255)      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                           >  1:265166/1‑100 (MQ=255)      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                           >  2:2002186/1‑100 (MQ=255)      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                           <  1:2444437/100‑1 (MQ=255)      ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaa                                                                                           >  2:2068788/1‑100 (MQ=255)        ccTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaaaa                                                                                         <  1:963238/100‑1 (MQ=255)         cTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaaaaa                                                                                        >  1:927541/1‑100 (MQ=255)         cTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTaaaaa                                                                                        <  2:918305/100‑1 (MQ=255)            gcggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGc                                                                                     <  1:1260161/100‑1 (MQ=255)            gcggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACAGCCGTCAATTAAAAAAGc                                                                                     <  1:488794/100‑1 (MQ=255)             cggcAGAAGGGTTCGAAGGAAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGc                                                                                     <  2:2532202/99‑1 (MQ=255)              ggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCgg                                                                                   >  1:45959/1‑100 (MQ=255)               gcATAAGTGTTCGAAGGTAGCATTACCGACCTGCCTAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  2:627330/100‑1 (MQ=255)               gcAGAAGTGTTCGAAGGTAGCTTTACCGACCTGACGAACCCGAAAGCTAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  2:1400975/100‑1 (MQ=255)               gcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  1:1859463/100‑1 (MQ=255)               gcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  1:1528173/100‑1 (MQ=255)                cAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCt                                                                                 >  2:2300168/1‑100 (MQ=255)                cAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCt                                                                                 >  1:387838/1‑100 (MQ=255)                      tgtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAc                                                                           >  1:636265/1‑100 (MQ=255)                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  1:1089660/1‑100 (MQ=255)                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  1:485371/1‑100 (MQ=255)                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  2:1893934/1‑100 (MQ=255)                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  1:1156281/1‑100 (MQ=255)                       gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAc                                                                           >  2:1060688/1‑99 (MQ=255)                        ttagAAGGTAGCATGACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgc                                                                         <  2:451213/97‑1 (MQ=255)                         tCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgcc                                                                        >  2:1099106/1‑100 (MQ=255)                               gTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtttt                                                                  >  2:2492459/1‑100 (MQ=255)                                tAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGTGGCTAACCACGCCGCttttt                                                                 >  2:833726/1‑100 (MQ=255)                                tAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttt                                                                 <  2:86275/100‑1 (MQ=255)                                  gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                               >  2:2284411/1‑100 (MQ=255)                                  gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                               <  1:1631763/100‑1 (MQ=255)                                  gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                               >  2:2346805/1‑100 (MQ=255)                                    aTTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTTTTTTTTTTGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTAc                                                             >  2:2413365/1‑100 (MQ=255)                                    aTTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTAc                                                             >  1:2413558/1‑100 (MQ=255)                                     ttACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTAc                                                             >  2:1725008/1‑99 (MQ=255)                                      tACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTTTTTTTTTTTTACATCCGTCAATTAAAAAAGCGGCTACCCACGCCGCTTTTTTTACGt                                                           >  2:2403312/1‑100 (MQ=255)                                       aCCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTc                                                          >  1:2211968/1‑100 (MQ=255)                                        ccGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCt                                                         >  2:1330330/1‑100 (MQ=255)                                         cGACCTGCCGAACCCGAAAGTGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTg                                                        <  1:51384/100‑1 (MQ=255)                                         cGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTg                                                        <  1:285265/100‑1 (MQ=255)                                         cGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTg                                                        <  2:1180311/100‑1 (MQ=255)                                          gACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGc                                                       >  1:2116157/1‑100 (MQ=255)                                              tGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAAtt                                                   <  2:1834785/100‑1 (MQ=255)                                               gCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAAtt                                                   <  2:1541153/99‑1 (MQ=255)                                                ccGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTa                                                 <  2:1089652/100‑1 (MQ=255)                                                ccGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTa                                                 <  2:1156273/100‑1 (MQ=255)                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                                <  2:2387680/100‑1 (MQ=255)                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                                <  1:2321261/100‑1 (MQ=255)                                                 cGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                                >  2:1308585/1‑100 (MQ=255)                                                      ccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCCCGCCGCTTTTTTTACGTCTGCACTTTACCTTTc                                           >  1:1491767/1‑100 (MQ=255)                                                      ccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTc                                           <  2:701580/100‑1 (MQ=255)                                                      ccGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTc                                           >  1:2364266/1‑100 (MQ=255)                                                       cagaaGCGAAGGAGAATTAATTTTTCGTCTGCCGGAACATCCGTGAAGTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          <  2:2398007/97‑1 (MQ=255)                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          >  2:869577/1‑100 (MQ=255)                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          >  1:1110287/1‑100 (MQ=255)                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          <  2:503160/100‑1 (MQ=255)                                                       cGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          <  1:2337642/100‑1 (MQ=255)                                                        gAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCa                                         >  2:1699083/1‑100 (MQ=255)                                                         aaaGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAg                                        <  2:2046635/100‑1 (MQ=255)                                                            gCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtct                                     >  2:462396/1‑100 (MQ=255)                                                             cGAATAAGAAGTAATTTTACATTTGCAAGAACATCCATCAATTAAATAAGCTGATAAATACGCCGCTTTTTTTACATCTACAATTTACCTTTGCattat                                     >  2:268852/1‑95 (MQ=255)                                                              gaagaagaaGTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTGTCCAGtctgc                                   <  2:2323797/100‑3 (MQ=255)                                                               aagaagaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  >  1:354997/1‑100 (MQ=255)                                                               aagaagaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  >  1:354926/1‑100 (MQ=255)                                                               aagaagaagTAATTTGTCGTTTGCCGCGACATCCGTCAATTAACAAAGCGGCTAACCACGCCGCTGTTGTTACGTCTGCAATTTACCTTTCCAGtcttct                                  <  2:1161569/100‑1 (MQ=255)                                                                agaagaagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTcttctt                                 <  1:492591/100‑1 (MQ=255)                                                                   agaagTAATTTTTTGTTTTCTTGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCt                              >  2:2060610/1‑100 (MQ=255)                                                                     aagTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTcc                            <  2:881879/100‑1 (MQ=255)                                                                      agTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCa                           <  1:940193/100‑1 (MQ=255)                                                                      agTAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCa                           >  1:1384648/1‑100 (MQ=255)                                                                        tAATTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACg                         <  1:737414/100‑1 (MQ=255)                                                                        tAATTTTTCGGTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACg                         <  1:180274/100‑1 (MQ=255)                                                                          aTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGtt                       >  2:2200885/1‑100 (MQ=255)                                                                          aTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGtt                       >  2:1560481/1‑100 (MQ=255)                                                                            ttttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCa                     <  1:1767115/100‑1 (MQ=255)                                                                            ttgtCGTTCGCCTGACCATCCTTCAATTACAAACGCGCTTTCCCACTCCGCTTTTTTTACGCCTGCCATTTTCCTTCCCACTCTTCTTGCTCCACCTTCa                     <  2:1585972/97‑1 (MQ=255)                                                                              ttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTACCCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCaga                   <  2:2472253/100‑1 (MQ=255)                                                                                cGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCCCGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCgga                   <  2:2293351/98‑1 (MQ=255)                                                                                cGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagaga                 <  1:2200900/100‑1 (MQ=255)                                                                                cGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagaga                 <  1:1980453/100‑1 (MQ=255)                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag                <  1:35891/100‑1 (MQ=255)                                                                                   ttGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGAc              <  1:29940/100‑1 (MQ=255)                                                                                     gCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGt            >  2:2069368/1‑100 (MQ=255)                                                                                      ccGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGtt           >  1:175324/1‑100 (MQ=255)                                                                                      ccGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGtt           >  2:2255916/1‑100 (MQ=255)                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg         >  2:1266068/1‑100 (MQ=255)                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg         >  1:1329065/1‑100 (MQ=255)                                                                                          aaCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCa       <  2:163954/100‑1 (MQ=255)                                                                                          aaCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCCTCTTGCTCCACGTTCAGAGAGACGTTCGCa       <  2:1713777/100‑1 (MQ=255)                                                                                           aCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCAt      <  1:2068802/100‑1 (MQ=255)                                                                                           aCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCAt      >  1:1644567/1‑100 (MQ=255)                                                                                              tCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAct   >  2:351028/1‑100 (MQ=255)                                                                                               ccGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATActg  >  2:552460/1‑100 (MQ=255)                                                                                                  |                                                                                                 CTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATACTG  >  NC_000913/129263‑129456

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