Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,968,501	C→A	Q159H (CAG→CAT)	rppH ←	RNA pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,968,501	0	C	A	96.3%	88.7 / NA	27	Q159H (CAG→CAT)	rppH	RNA pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (6/20);  minor base G (0/1);  total (6/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAACGCACCAAACG  >  NC_000913/2968264‑2968734                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          cccccGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAAACAGAACAGACCTCGGTATCCATCGCCAGACAGATGTCGGTCACCACAATCTTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTAAGTCGCTAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGAcctc                                                                                                                                                                                                                               >  2:219185/1‑247 (MQ=255)       ggTCGCCATCAGGTAGTAACAACGTCGATCGTGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCacac                                                                                                                                                                                                                        <  1:188994/251‑1 (MQ=255)             cATCAGGTAGTAACAACGTCGATCATGATCGGGTAGGTAGACCGGACAGACCTCGGTATCGATCTGAAGGCAGATGTCGGTAACCAGAAGATTGAACGCCGCAGGCAGGCGTTGTGCGCTGGCTACCGTTTCGACTATTTCGCGTAGGCGAGTGAGCATAATTGGCGTGTCTTAATCGCTTTTACTTGGATAGGCGGGTGAGTGGTGTGGTTTTGGTGTATGTTCGTGGTGTGACAGCAGCACACGCGCGa                                                                                                                                                                                                                  <  2:39971/251‑1 (MQ=255)                aagggagtaACAACGTCGATCAAGATCGGGCCGGTAGACAGAACAGACCTCTGTATACATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCGCATTCAGGCGTGGTGCGCTGGCTACCTTTTCTAATTGTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTATCCTCTTTTGCGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGGAGTTACATCACCACACTCGCGacc                                                                                                                                                                                                                <  2:142644/245‑3 (MQ=255)                             cGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCt                                                                                                                                                                                                  <  2:170903/251‑1 (MQ=255)                              gTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTa                                                                                                                                                                                                 <  1:94116/251‑1 (MQ=255)                                       gATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGAc                                                                                                                                                                                         >  1:112143/1‑250 (MQ=255)                                                   aGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAAt                                                                                                                                                                            <  1:43440/251‑1 (MQ=255)                                                                                   cAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACt                                                                                                                                            >  1:102045/1‑251 (MQ=255)                                                                                      aTGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTGCAGTACCCTACGGTAGACATCACGTTTAAATGACACCCCCTGTCTTACCGGATCCCAGTAACTTAc                                                                                                                                         <  2:144490/251‑1 (MQ=255)                                                                                                              cccTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCAGAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTgg                                                                                                                 <  2:34433/250‑1 (MQ=255)                                                                                                                  cATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTa                                                                                                             <  2:69121/251‑1 (MQ=255)                                                                                                                                gcgcTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCa                                                                                                <  1:136056/250‑1 (MQ=255)                                                                                                                                     ggCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCCCGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATAtt                                                                                            <  2:219876/249‑1 (MQ=255)                                                                                                                                           cTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTACCCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATAACAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACGGCTGGTTTGCATATTGATTTCTg                                                                                    >  1:57536/1‑251 (MQ=255)                                                                                                                                                     aTTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATcgccg                                                                            <  2:179669/249‑1 (MQ=255)                                                                                                                                                              aGGCGAGGGAGCATAAGTGGCGTGACTTAACCTCTGTTACGTCGATAAGCAGATGCGTTTTGTGGTGTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTGGCAGATTGATTTCTGCATCGCCGCTCACCAGCTg                                                                 <  2:62258/251‑1 (MQ=255)                                                                                                                                                               ggCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGc                                                                <  1:23636/251‑1 (MQ=255)                                                                                                                                                                 cGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCa                                                               >  1:2155/1‑250 (MQ=255)                                                                                                                                                                    gTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAagag                                                           <  1:14411/251‑1 (MQ=255)                                                                                                                                                                             aTTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCAttt                                                  <  1:173808/251‑1 (MQ=255)                                                                                                                                                                                            gctctTTTAAGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGACAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACAAGCTGCAAGAGAAACCATTTTTGTTGTTGGCCGAt                                   >  2:78985/2‑251 (MQ=255)                                                                                                                                                                                               ctTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATa                                  <  2:65233/249‑1 (MQ=255)                                                                                                                                                                                                               cAGATGCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTGTTGTTTTTGGCCGATACAAACCGGCTTCGTGTcc                <  2:112143/251‑1 (MQ=255)                                                                                                                                                                                                                    gCGTTTTGTGGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGTATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGTCTTCGTGTCCCAAc            <  2:11947/250‑1 (MQ=255)                                                                                                                                                                                                                           gtgGTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTGTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAGCGCACCaa     <  2:178907/250‑1 (MQ=255)                                                                                                                                                                                                                             ggTTTTGGCGTATTTTCATGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAACGCACCAAACg  <  1:233/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          CCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCGAACTCTTTCATTACCCTACGGTAGACATCACGTTTAAATGACACCACCTGTCTGACCGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAACGCACCAAACG  >  NC_000913/2968264‑2968734

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