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	96.2%	80.5 / ‑4.3	26	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (17/8);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815660‑3816132                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          cGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                  >  2:43876/1‑250 (MQ=255)         aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTACTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGTTGAGTAGCTCTTc                                                                                                                                                                                                                          >  2:226802/1‑250 (MQ=255)             aGGGGCTTTTGCGCCCGCATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAACGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACAATGGCTTCGAT‑TCCCCCTCGGCCCAGAGCCACCAAGATGAGTAGCTCTTCATggg                                                                                                                                                                                                                     >  2:294152/1‑251 (MQ=255)                    tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                               <  2:172732/250‑1 (MQ=255)                       gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                                                                                                                                            <  2:104533/250‑1 (MQ=255)                                 ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 >  1:207538/1‑251 (MQ=255)                                                      tatAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                            >  2:195908/3‑251 (MQ=255)                                                           cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGccc                                                                                                                                                                       >  2:323288/1‑251 (MQ=255)                                                                                       atcaggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGACGAGTAGCTCTTCATGGGTGAACGGCGCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                           <  1:88814/247‑1 (MQ=255)                                                                                                 gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                                                                                                                                                      >  1:23001/1‑166 (MQ=255)                                                                                                 gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                                                                                                                                                      <  2:23001/166‑1 (MQ=255)                                                                                                 gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                                                                                 <  1:43876/251‑1 (MQ=255)                                                                                                         ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                         >  2:299446/1‑251 (MQ=255)                                                                                                             tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAACTCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTCCAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTAAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                                                                                                                      >  2:201309/1‑250 (MQ=255)                                                                                                               ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                                                                                   <  1:18869/251‑1 (MQ=255)                                                                                                                  tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGt                                                                                                                  <  1:205943/249‑1 (MQ=255)                                                                                                                              taaCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGTCTCGCCTTCTCCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGTCCATCACTACGTTCATGTCGGTCTCTGCAGCAGAGTCTTCAACGTATTCCAGacagca                                                                                                     >  2:237178/3‑245 (MQ=255)                                                                                                                                            aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                                                                             >  1:282844/1‑244 (MQ=255)                                                                                                                                            aCAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCGACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                                                                             <  2:282844/244‑1 (MQ=255)                                                                                                                                                        cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  2:47886/1‑251 (MQ=255)                                                                                                                                                              atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAAGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAAc                                                                    >  1:31004/1‑251 (MQ=255)                                                                                                                                                               ccaGTCGCTTTAAAAATCACTGGGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAGGCCCCC‑TCGGGCCAGTGCCAACCAGATGAGTAGCTCTTCATGGGTGAACGGCGCGCCTTCTGCCTTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCCCTGCGGCAGATTCTTCAACGTATGCCAGATGTCAAACCGCTTCGCCTTTCACAAGTCCGCCAGAAACt                                                                   <  2:331660/250‑1 (MQ=255)                                                                                                                                                                 aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTGATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc                                                                 >  2:116796/1‑251 (MQ=255)                                                                                                                                                                                      gCCAGCGACGCCTTCTGCGTCGCTACAGTGGATTCGAT‑TCCCCCTCGGGACAGAGCCAACAAGATGAGTAGAGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGTGCCTCACTACGTTCATGTCGGTCCCTGCGGCAGAGTCTTCAACGTCTTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGCCAGAACATGCGGCTACCATTCCTTTCAACgg                                            >  2:179542/1‑251 (MQ=255)                                                                                                                                                                                           cgccgcCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCAGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                       >  1:8310/1‑251 (MQ=255)                                                                                                                                                                                                           gCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtt                       >  1:347408/1‑251 (MQ=255)                                                                                                                                                                                                                                 c‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  1:250209/1‑250 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          CGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815660‑3816132

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