Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 218,835 (A)8→7 intergenic (‑60/+52) proS ← / ← trmO proline‑‑tRNA ligase/tRNA m(6)t(6)A37 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913218,8280A.96.5% 113.9 / ‑1.8 29intergenic (‑53/+59)proS/trmOproline‑‑tRNA ligase/tRNA m(6)t(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (13/15);  minor base C (1/0);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAG  >  NC_000913/218693‑218958
                                                                                                                                       |                                                                                                                                  
ggcaggaTCAGCGGATGGCCGATCAACTCGGCGGAGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                <  2:402383/139‑1 (MQ=255)
cgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                >  2:756541/1‑140 (MQ=255)
cgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                >  2:749015/1‑140 (MQ=255)
cgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                >  2:695756/1‑140 (MQ=255)
cgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                >  2:609335/1‑140 (MQ=255)
     caTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaaagtg                                                                                                                           <  1:749015/140‑4 (MQ=255)
      aTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaaagtga                                                                                                                          <  2:590913/140‑5 (MQ=255)
        cAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGAtt                                                                                                                       <  1:459985/140‑1 (MQ=255)
        cAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGAtt                                                                                                                       <  1:695756/140‑1 (MQ=255)
             gATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGtt                                                                                                                  >  1:18316/1‑140 (MQ=255)
                       cACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCa                                                                                                        <  1:648257/140‑1 (MQ=255)
                            cGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGa                                                                                                   <  1:23705/140‑1 (MQ=255)
                              gCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGaaa                                                                                                 <  2:486329/140‑1 (MQ=255)
                                     cAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCaaa                                                                                          <  2:877767/140‑1 (MQ=255)
                                     cAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCaaa                                                                                          <  2:18316/140‑1 (MQ=255)
                                         tgtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGt                                                                                              <  1:779618/132‑1 (MQ=255)
                                         tgtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGt                                                                                              >  2:779618/1‑132 (MQ=255)
                                                                     tGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAgc                                                          >  2:389716/1‑140 (MQ=255)
                                                                        cTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCa                                                       <  1:569049/140‑1 (MQ=255)
                                                                         tAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaa                                                      <  1:756541/140‑1 (MQ=255)
                                                                         tAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaa                                                      <  1:609335/140‑1 (MQ=255)
                                                                                             ttCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCg                                  >  2:681229/1‑140 (MQ=255)
                                                                                               ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGc                                      >  2:49425/1‑134 (MQ=255)
                                                                                               ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGc                                      <  1:49425/134‑1 (MQ=255)
                                                                                                cAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAGCCTGCGTCGGTg                               <  2:841476/140‑1 (MQ=255)
                                                                                                 aGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGa                              >  1:180198/1‑140 (MQ=255)
                                                                                                       aaGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGcc                        >  2:464276/1‑140 (MQ=255)
                                                                                                             gAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGaa                  <  1:915045/140‑1 (MQ=255)
                                                                                                                cAGGCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGt               >  1:537546/1‑140 (MQ=255)
                                                                                                                   gCTCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaa            >  2:852041/1‑140 (MQ=255)
                                                                                                                     tCAAGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaaaa          >  2:703548/1‑140 (MQ=255)
                                                                                                                        aGGCGAGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCa       >  2:659025/1‑140 (MQ=255)
                                                                                                                             aGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg  >  1:204886/1‑140 (MQ=255)
                                                                                                                             aGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg  <  2:537546/140‑1 (MQ=255)
                                                                                                                             aGCCTGGGACAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg  >  1:116710/1‑140 (MQ=255)
                                                                                                                             aGCCTGGGCCAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg  >  2:572680/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAG  >  NC_000913/218693‑218958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: