Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 433,649 T→G 29.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_000913433,6490TG29.9% 22.4 / 3.1 32G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (6/16);  new base G (6/4);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-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)
                                                                                                                                                                                                          gggttaggggggggaaaaaaGCCAAAGGGGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt             <  2:342379/224‑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 < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: