Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	433,649	T→G	27.9%	G65G (GGT→GGG)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	27.9%	23.9 / 2.5	31	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (6/16);  new base G (6/3);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 5.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCT  >  NC_000913/433419‑433870                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGAACg                                                                                                                                                                                                                        >  1:370558/1‑238 (MQ=255)        aaGAGAACGCGCGCATGAAGGACTAAGCCGTGCGGGACGGGTATTACATGGCGGGGGCGCTAAAGCTGGCGCCAAGAGGACGTTTTACCACGAATCCCAACCCGAGTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAGGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                 <  2:70904/238‑1 (MQ=255)                        ggaggaCAAAGCCGTGCAGGACGAGTATGCCAGGGCGCGGGCTGTAAAGCGGGCGCAACGAGGCCGTTTTCCCACGCTACCCACCCCGAATGTCGGGGGCGTCATTGTGAAAGATGGCGAAATTGTCGGTGAAGTTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGt                                                                                                                                                                                                 <  2:350628/238‑1 (MQ=255)                              cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGGGGACGCCCCAGTCACCCCCGCACCCTTTAGCcat                                                                                                                                                                                           >  2:287053/1‑238 (MQ=255)                                    cggtgaGGACGAGTATTGCAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                    <  2:166748/234‑1 (MQ=255)                                        gaGGACGAGTTTTACATGGCGGGGGCGCTAAAGCGGGCACAACGAGGACGTTTTACCACGCACCCCAACCCGAATGCCGGGTGCGTCATTGTCAAAGATGCCAAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                 <  2:255401/237‑1 (MQ=255)                                           gACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCTATGGTACAATCGAACCCTGTAGCCATCATGGTCGTACGc                                                                                                                                                                              >  1:337784/1‑238 (MQ=255)                                             cgggtcttACAGGGCGCGGGCCGTAAAGCGGGCGAACAGAGGGCGTTTCACCCCGCATCCCACCCGGAATGGCGGGGGCGTCATTGTCAAAGATGGCAAAATGGTCGGTGAAGGTTTCCACAAACGTGGGGGTGACCCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                            <  2:391009/232‑1 (MQ=255)                                                gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGTGCGACCGCCGATGTCACACTCGAACCCCGTAGCCATCATGGGCGTGAGCCACCgg                                                                                                                                                                        >  1:211973/1‑238 (MQ=255)                                                 tATTACATGGGGCGGGCGTTAAAGCGGGCGCAACGAGGACGTTTTACCACGCAGCCCAACCCGAATGTCGGGGGCGCAATTGACAAAGAGGGCAAAATTGTGGGGGAAGGTTACCCCCAACGTGCGGGTGAACCACAGGCCGAGGTACACGGGTGGCGTATGGCGGGGGAAAAACCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                       <  2:262926/239‑1 (MQ=255)                                                     aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGCGGAAAAAACCAAAGGGGCGGCCGCCTATGTTCCACTCGGACCCTGTAGGCACCAACGGCGCACGCCACCGTGCtgt                                                                                                                                                                   >  1:68769/1‑239 (MQ=255)                                                           cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCCCACTCGAACCCTGTAGCCATCCTTGTCGTACGCCCCCGTGCTGTGACGc                                                                                                                                                              >  2:178563/1‑238 (MQ=255)                                                                  ccTAAAGCTGGCGAAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTAACCCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATc                                                                                                                                                       <  2:383641/237‑1 (MQ=255)                                                                                             tACCAAGCATCCTACCCCGAATGTCGGGTGCGTCATTGTCAAGGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                           <  2:333928/239‑1 (MQ=255)                                                                                               ccACGCACCCCAACCCGAATGGCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                                                                                                         <  2:245655/239‑1 (MQ=255)                                                                                                    cATCCCAACCCGAAGGTGGGGGGCGGCATTGTCAAAGAGGGCGAAATTGGCGTGGAAGGTTCCCCCCAGCGTGCGGGTGAACCCAATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                    <  2:337784/239‑1 (MQ=255)                                                                                                     aacccagcgggAATGTGGGGTGCGTAATGGTCAAAGATGGCGAAGTTGTGGGTTAAGTTTACCACAAAGGTGCGGGTGAACCACATGCGGAAGTACACGCGTTGCGTATGGCGGGTGAAAAGGCCAATGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                    <  2:68769/228‑1 (MQ=255)                                                                                                      tCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCACAAGGGCGCCCCCCCCCGGCACCCTCTGACCCTCTACCCCTCATGGGCGTCCGCCCCCGCCCTGTGACCCCCTCCGCGCCCCCGCCGTGGAGCGCGCGGTTACCTcccctgacg                                                                                                                  >  2:254920/1‑231 (MQ=255)                                                                                                                   gtggggTGCGGAATTGTAAAAGGTGGCAAAATTGGCGGTGAGGTTTACCGCCACCGTGCGGGTGAACCACATGCGGAAGTCCCCGCGTGGCGTAGGGGGGGTGAAAAAGCCAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc                                                                                                      <  2:132066/235‑1 (MQ=255)                                                                                                                                             cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCCCCTATGTCACACCCGAAACCTGTAGCCCTCATGGGCGTACGCCACCGTGCTGTGACGCACTCATTCCCGCCCGCGTCACGCGCGTGGTTGCCTCGATGCCAGAACCTAACCCCCGGGTGCCTGGGCGTGGaacttaacct                                                                            >  2:2056/1‑229 (MQ=255)                                                                                                                                                 aTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCAAATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCa                                                                       <  1:287053/239‑1 (MQ=255)                                                                                                                                                        gTGAAGGTTACACCCAACGTGCGGGGGACCGACATGCGGAAGTACACGGGTTGCGTATGGCGGGGGAAAAAGCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                                 <  1:34361/238‑1 (MQ=255)                                                                                                                                                                        aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCCAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCCTCCTGGTCGTACGCCACCGTGCTGTGACGCACTCCTCCCCGCCGGCGTAGCGCGCGGGGTTGCCCCGATGCAAGATCCCAAACCGCCGGGCGCTGGGCGGGGGATTTCCCGTCCGCACCCGGCCGGCCCTGACCTCCGc                                                 >  2:292427/1‑238 (MQ=255)                                                                                                                                                                          ggggggcgaaaCCCCAAGCCGAAGTACGCGCGTGGCGTATGGCGGGGAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                               <  1:178563/229‑1 (MQ=255)                                                                                                                                                                           tGCGGGGGAACCACAGGCGGAAGTACACGCGTTGCGTATGGCGGGTGAAAAGGCCAAAGTGGCACCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                               <  1:16329/237‑1 (MQ=255)                                                                                                                                                                            gggggTGAAACACATGCCGAAGTACACGCGGTGCGTAGGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTCCGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                             <  2:471326/236‑1 (MQ=255)                                                                                                                                                                                                aGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGCCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCCCGATGCAAGAGCCTAACCCCCGCGTCCCTGGGCGTGGCCTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCCGATGGTGAGTGAAGcca                        >  2:50828/1‑238 (MQ=255)                                                                                                                                                                                                     acggggtGCGTAGGGGGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                   <  2:370558/233‑1 (MQ=255)                                                                                                                                                                                                      ggctttGCGTATGGCGGGTAAAAATCCAAAAGGTTCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt                  <  2:165296/235‑1 (MQ=255)                                                                                                                                                                                                                       gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGGACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGGGGACCTTACCCGCTGCAAAAGGCTGGCATTGACCTCAGCCCCGGCCCGAGTGTGAGTGAAGCCGAGCAATTTAATAAAAGCGTtct  >  2:322576/1‑238 (MQ=255)                                                                                                                                                                                                                       gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtct  >  2:460751/1‑238 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCT  >  NC_000913/433419‑433870

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