Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,404,215	0	T	G	44.3%	92.9 / 13.1	62	G244G (GGT→GGG)	acuI	putative acryloyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (7/27);  new base G (27/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 3.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGG  >  NC_000913/3404071‑3404363                                                                                                                                                  |                                                                                                                                                      cccTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGtttt                                                                                                                                                 <  2:1235166/151‑1 (MQ=255)  ccTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCCGTGGTCTGGC‑GGGGGGgttt                                                                                                                                                 >  1:1054082/1‑150 (MQ=255)  ccTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGGCTGGC‑GGGGGGgttt                                                                                                                                                 >  2:942474/1‑150 (MQ=255)     cGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGgttt                                                                                                                                                 >  2:997552/1‑147 (MQ=255)     cGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTct                                                                                                                                              >  2:1153865/1‑150 (MQ=255)         atgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGcc                                                                                                                                        <  1:912131/151‑1 (MQ=255)         atgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑TGcc                                                                                                                                        >  2:33012/1‑151 (MQ=255)          tgggTTTTCCTAATCCCGTCCTCTGGAAAAACAAGTATGGGCTGGGGCAAATGACACCGTTGGCGACAAAGCGCTGGCAAAAGTGCTGGCGCGAAGGAAGTGCGGCGGCGGCGTGGAGGAAAAAGGTCTCGA‑GGGTGGtttttttct                                                                                                                                            <  2:209487/144‑6 (MQ=255)                 gCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACgg                                                                                                                                <  1:894882/151‑1 (MQ=255)                 gCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACg                                                                                                                                 <  2:1092518/150‑1 (MQ=255)                 gCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACg                                                                                                                                 <  2:574213/150‑1 (MQ=255)                    gAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCa                                                                                                                             <  2:873724/151‑1 (MQ=255)                    gAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATGGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACgggc                                                                                                                              <  2:85556/150‑3 (MQ=255)                       tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGTTTTTCTC‑TGCCCAACACCGGCATGc                                                                                                                          >  2:1215270/1‑151 (MQ=255)                                 tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGGTTTTCTC‑TTCCCACCACGGGCATGGCATTTAtttt                                                                                                                >  1:15870/1‑149 (MQ=255)                                  ggAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTAtt                                                                                                                  <  1:1103767/148‑1 (MQ=255)                                      aaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGGTGCGTGGCCGGCTGTAGTGTGGC‑GGGGGGGTTTAATC‑TCCCAACCAAGGACATGCTATTTATTCTGCg                                                                                                             >  1:537901/1‑149 (MQ=255)                                        aCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAAt                                                                                                         <  2:116805/151‑1 (MQ=255)                                        aCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAACCACGGGCATGCCATTTATTTTGCGTaa                                                                                                          >  2:992501/1‑150 (MQ=255)                                          aaGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAAACACGGGCATGCCATTTTTTCTGCggta                                                                                                          >  2:390672/1‑145 (MQ=255)                                                gggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGttt                                                                                                 <  2:344205/151‑1 (MQ=255)                                                 ggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGttt                                                                                                 >  1:1200698/1‑150 (MQ=255)                                                     ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCa                                                                                              <  2:671874/149‑1 (MQ=255)                                                        gCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGG‑GGGGGGGTTTTCTCTTGCCAA‑CACGGGCATGCCCTTTTTTTTGCGGAAATGCCGGTTGCAAgggg                                                                                         >  1:342079/1‑151 (MQ=255)                                                          aaTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑TGCCAACCACGGTCATGCCCTTTTTTTTGCGTAATTTCCGTTTTCCAgggggg                                                                                       >  2:168966/1‑149 (MQ=255)                                                              gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGAtt                                                                                   <  2:185885/151‑1 (MQ=255)                                                              gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGAt                                                                                    <  2:15870/150‑1 (MQ=255)                                                                     ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGC‑GGGGGGTTTTTCCC‑TTCCCACCAACGTAATAGCACTTATTTTGCGGAATTTCCGGTTGGAAAGGGGGGGTTCAGCaa                                                                             >  1:1125734/1‑150 (MQ=255)                                                                      tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑TGCCAACCACGGTCATGCCATTTTTTCTTCGTTATGTCCGTTTGCAAGGGGGGGGTTCAGTaa                                                                             >  2:618727/1‑149 (MQ=255)                                                                       ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑CGCCAACCACGGGCATGCCAATTTTTTTGCGGAATGTCCGGTTGCAAGgggggg                                                                                      >  2:784641/1‑139 (MQ=255)                                                                             aaaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATTT‑TGCCAA‑CACGGGCATGCCCTTTTTTTTGCGTAATGTCCGTTTTCAAGGGGGGGGTTTAGTAAAGatggc                                                                      >  2:657288/1‑145 (MQ=255)                                                                               aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGccacca                                                                  <  2:526860/151‑1 (MQ=255)                                                                                    tGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAACCACGGTCATGCCCTTTATTCTGCGTAATGTTCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGaa                                                               >  1:331722/1‑149 (MQ=255)                                                                                      gCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTTTC‑TGCCCAACACGGGGAAGGCATTTATTTTTCGGAATTTTCGGTTTCCAGGGGGGGGGTCAGGAAAAACGCCCCCCGaa                                                               >  2:368113/1‑147 (MQ=255)                                                                                       cAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAACCACGGGCATGCCCTTTTTTCTGCGGAAAGTCCGTTTGCAAGGGGGGGGTTCAGTAATGAAGCCACCAGAAcgcc                                                           >  2:437408/1‑150 (MQ=255)                                                                                          aaGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑TGCCAACCACGGGCATGCCATTTTTTCTGCGGAATGTCCGGTTGCAAGGGGGGGATTTAGTAATGGAGCCCCCAGAAcgccgc                                                         >  2:1195442/1‑149 (MQ=255)                                                                                            gTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTC‑TGCCCACCACCGGCATGCCATTTTTTTTGCGTAATGTCCCTTTTCAAGGGGGGGGTTCAGTAATGACGCCACCAGAACGCCGCGca                                                      >  2:864435/1‑150 (MQ=255)                                                                                              gCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGcaca                                                    <  1:844107/150‑1 (MQ=255)                                                                                              gCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTATC‑TGCCAACCCCGGTCATGCCCTTTTTTCTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAAAGCCCCGcac                                                     >  2:1244758/1‑149 (MQ=255)                                                                                               cTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTATC‑TGCCCACCACGGTCATGCCCTTTATTTTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAACGCCGCGcac                                                     >  2:613912/1‑148 (MQ=255)                                                                                                 ggCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTATTC‑TTCCCAACCCGGGCATGCCCTTTATTCTTCGTAAAGTCCGTTTGCCAGGGGGGGGTTTAGTAATTACGCCCCCAGAACGCCCCGCCCAAGcc                                                >  2:1098974/1‑151 (MQ=255)                                                                                                  gcgcAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCt                                               <  1:470245/151‑1 (MQ=255)                                                                                                    gcAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTgg                                             <  2:221579/151‑1 (MQ=255)                                                                                                    gcAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGTTTTTCTC‑TGCCAAACACGGGCATGCCATTTTTTTTTCGTAATGTCCGTTTGCAAGGGGGGGGTTTAGTAATGACGCCACCCGAAAGCCGCGCACAAGcc                                                >  1:988459/1‑148 (MQ=255)                                                                                                     cAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTC‑TGCCAACCACGGTCATGCCATTTATTTTGCGTAATGTCCGTTTGCAAAGGGGGGATTCAGTAATTACGCCACCAGAACGCCGCGCACAACCCCgg                                             >  1:502169/1‑150 (MQ=255)                                                                                                      aaaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGc                                            <  2:1253330/150‑1 (MQ=255)                                                                                                               aCGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGGTTTACTC‑TGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTgg                                   >  1:213030/1‑150 (MQ=255)                                                                                                                cggcggCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTc                                 <  1:1082497/151‑1 (MQ=255)                                                                                                                  gcggcTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGc                               >  1:1156126/1‑151 (MQ=255)                                                                                                                  gcggcTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGc                               <  1:833654/151‑1 (MQ=255)                                                                                                                       tGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAtt                          >  1:327260/1‑151 (MQ=255)                                                                                                                         cGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAttt                         >  2:65393/1‑150 (MQ=255)                                                                                                                         cGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTa                        <  1:289821/151‑1 (MQ=255)                                                                                                                          gTGGCGGCCTGTGGACTGGC‑GGGGGGTTTTACTC‑TGCCAACCACGGTCCTGCCATTTTTTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGGAATGACGCCAACAGAACGCCGCGCACAAGCCTGGCAGCGACTGGGCGCCGAttt                         >  2:1128589/1‑149 (MQ=255)                                                                                                                               ggcCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa                  <  2:933903/151‑1 (MQ=255)                                                                                                                                gcCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt                 <  1:40167/151‑1 (MQ=255)                                                                                                                                        tCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCta          >  1:1126215/1‑150 (MQ=255)                                                                                                                                           ggC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCtata        <  1:221640/149‑1 (MQ=255)                                                                                                                                           ggC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCtata        <  1:743785/149‑1 (MQ=255)                                                                                                                                           ggC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCtat         <  1:280653/148‑1 (MQ=255)                                                                                                                                            ac‑GGGAGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATAccc     >  2:920510/2‑151 (MQ=255)                                                                                                                                               gggTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATAccc     >  2:46205/1‑149 (MQ=255)                                                                                                                                                ggtggtTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAgg  <  1:121671/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTC‑TGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGG  >  NC_000913/3404071‑3404363

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