Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	89.4%	135.6 / 8.5	47	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (27/15);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

GTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGC‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCT  >  NC_000913/3815644‑3816125                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ttATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTACGC‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTACGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAAGGCGCATCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCGTAGATACAATGGATACG‑T‑T‑Tcc‑cccccggcccatc                                                                                                                                                                                                                                          >  2:203962/2‑240 (MQ=255)          aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGC‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAg                                                                                                                                                                                                                                 <  1:245549/251‑1 (MQ=255)                     gaagaaATTGGGGGTTTTG‑GCGCCGGATTTCCGCCTAAACCCCGCCAA‑CTTTAACCCCCTCTTGCTAAACGCAAAATCAATAAAATGGCGCTGATATTGTTTCCTGTCTTCGCTCCCTATCTTTCCTTTCCACAGAAAAAAAAAAGGCGGCCCATCCGGCGCCTTAAAAAACAATTTGCCATCCCCCCCTTCCTCGTCGATTCAATTGATTCGGT‑T‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGct                                                                                                                                                                                                                       <  1:283574/250‑1 (MQ=255)                                       gcgc‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                                                                                                                                     <  2:169239/250‑1 (MQ=255)                                        cgc‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                                                                                                                                   <  1:347039/251‑1 (MQ=255)                                          c‑ccGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                                                                                                                                 >  1:148615/1‑251 (MQ=255)                                           accGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTcccc                                                                                                                                                                                                 >  2:118025/2‑248 (MQ=255)                                             cGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                               >  1:26335/1‑251 (MQ=255)                                                               c‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                             >  1:381426/1‑251 (MQ=255)                                                                          tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGa                                                                                                                                                                                                                             >  1:352981/1‑193 (MQ=255)                                                                          tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGa                                                                                                                                                                                                                             <  2:352981/193‑1 (MQ=255)                                                                             cgccTGCTTGCTAAGCGCAAATTCAAGAAACTGGCGCTGATATGGTTTCATTCCGTCGCACCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCTCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑T‑T‑TCC‑CCCTCGGGCCAGTGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                                                                <  2:256443/249‑1 (MQ=255)                                                                              acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCGGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCg                                                                                                                                                               >  1:248467/1‑251 (MQ=255)                                                                              acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCg                                                                                                                                                               >  1:218745/1‑251 (MQ=255)                                                                                             ccacaTTCAATAAACTGGCGCTGATATGGTTTCATGCCATCGCCCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGACGCAACAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                                <  2:163618/247‑1 (MQ=255)                                                                                                             ggCGCTGAGATGTGTTCATGCCTTCTCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACGCATCAGTTGCCTTAAAAAGCAGTCTGCCAGCGCCGCCTTCTGCGTCTCTACAATGGAGTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGtct                                                                                                                                <  2:381960/251‑1 (MQ=255)                                                                                                              gcgcTGATATGGTCTCATTCCTTCGCCCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACCCTGCAGGCTCCTTAAAAACCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCG‑A‑T‑CCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCTCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                                                                                                               <  1:461037/251‑1 (MQ=255)                                                                                                                gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                                                                                                      <  1:254857/210‑1 (MQ=255)                                                                                                                gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                                                                                                      >  2:254857/1‑210 (MQ=255)                                                                                                                                gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                                                             >  1:443806/1‑251 (MQ=255)                                                                                                                                    tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                                                                                                         >  1:180594/1‑251 (MQ=255)                                                                                                                                                 tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaa                                                                                             >  1:90619/1‑250 (MQ=255)                                                                                                                                                  acgtttCTACAGACAAAAAAAATGCGACTAATCAGTCGCTTTAAAAATCAGTGTGCCAGCGCCGCCTTCTGCGTCGCTACCATGGATGCG‑A‑TGCCC‑CC‑TCGGGCCAGAGCGAACAAGATGAGTAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTTCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTGCAGATCGCAAAc                                                                                           <  2:280584/248‑1 (MQ=255)                                                                                                                                                          aCAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGACGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTATTCCAGATCGAAACCCGCTTCGc                                                                                   >  2:303213/1‑251 (MQ=255)                                                                                                                                                           cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                                          >  1:440434/1‑195 (MQ=255)                                                                                                                                                           cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                                          <  2:440434/195‑1 (MQ=255)                                                                                                                                                                      aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                                                                       >  1:111274/1‑251 (MQ=255)                                                                                                                                                                          cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGAc                                                                   >  1:360822/1‑251 (MQ=255)                                                                                                                                                                          cGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGACATCACTACGTTCCTGTCAGTCTCTGCGACAGAGGCTTCAACAGATTACAGATCGCAAACCCCTTCCCCGTTCACAATTCCGAc                                                                   >  2:305354/1‑251 (MQ=255)                                                                                                                                                                           gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGAAGAGTATTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                   >  1:248780/1‑250 (MQ=255)                                                                                                                                                                               catccGTCGCATTATAAATCATTTTGCCAGCGCCGACTTCTGCGTCGCTACAAAGGATTCG‑A‑T‑TCCACC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAATCTTCAACGTCTTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCAGAccgca                                                               >  2:446346/1‑246 (MQ=255)                                                                                                                                                                                atcaGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa                                                                                              >  1:324635/1‑218 (MQ=255)                                                                                                                                                                                atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa                                                                                              <  2:324635/218‑1 (MQ=255)                                                                                                                                                                                  ccgTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                                                           <  2:243496/249‑1 (MQ=255)                                                                                                                                                                                   aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGAAACTGc                                                          >  2:4496/1‑251 (MQ=255)                                                                                                                                                                                   aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc                                                          >  1:258460/1‑251 (MQ=255)                                                                                                                                                                                    gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGACAGAAACTGCg                                                         >  2:440118/1‑251 (MQ=255)                                                                                                                                                                                     tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                        >  1:348502/1‑251 (MQ=255)                                                                                                                                                                                       gCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                      >  1:290562/1‑251 (MQ=255)                                                                                                                                                                                          ttAAAAATCAGTTTGCCAGCGACGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCCACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAcc                                                   >  2:442593/1‑251 (MQ=255)                                                                                                                                                                                           tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAcc                                                   >  2:171217/1‑250 (MQ=255)                                                                                                                                                                                              aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc                                               >  1:260188/1‑251 (MQ=255)                                                                                                                                                                                                            gcgcCGCCTTCTGCGCAGCTACAATGGATTCG‑ACT‑CCC‑CC‑GCGGGCCAAGGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTGTCATCGGATTg                                 <  2:97161/251‑1 (MQ=255)                                                                                                                                                                                                              gccgccTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAAATGCGGCTACCATCCCTTTCATCGGATTGGtt                              >  1:195297/1‑251 (MQ=255)                                                                                                                                                                                                                     tCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCa                          <  2:348502/249‑1 (MQ=255)                                                                                                                                                                                                                        gCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtt                      <  1:8294/250‑1 (MQ=255)                                                                                                                                                                                                                         cGTCGCTACAATGGATTCG‑A‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGCAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                    >  2:13315/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaTGGATTCG‑A‑T‑TCC‑CCCGCCGTCCAGAGCAAACAAGATGAGTAGCTCTTCATGGTTGAACTGCTCGCCTTCTTCCGTCCCTCGCACTTCAATGATGCGACCGTCTTCGGTTATCACTACGTTCATGTCGGTCGCGGCGGCAGAGTCTGCAACGTATTCCAGCTCGCAAACCGCTTCGCCGTTCACAATTCCGACATAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcca           <  2:415622/251‑1 (MQ=255)                                                                                                                                                                                                                                             gattccc‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCt  >  2:479386/6‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGC‑CCGGATTTCAGCGTAAACTC‑GCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑A‑T‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCT  >  NC_000913/3815644‑3816125

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