Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 218,835 (A)8→9 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,8271.A97.5% 137.9 / ‑3.4 40intergenic (‑52/+60)proS/trmOproline‑‑tRNA ligase/tRNA m(6)t(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (26/13);  total (27/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

CATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑‑AAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCA  >  NC_000913/218698‑218961
                                                                                                                                  |                                                                                                                                       
caTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑‑aaaaaaaaag                                                                                                                              >  1:910481/1‑138 (MQ=255)
             cTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTAc                                                                                                                  <  2:333448/140‑1 (MQ=255)
                aTCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAg                                                                                                               >  2:146775/1‑140 (MQ=255)
                               gCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTc                                                                                                >  1:141430/1‑140 (MQ=255)
                                                         cAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAAcgc                                                                      >  1:849684/1‑140 (MQ=255)
                                                            gTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGt                                                                   >  2:726005/1‑140 (MQ=255)
                                                              aTTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTc                                                                 >  1:748249/1‑140 (MQ=255)
                                                                 ggCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAg                                                              >  1:349929/1‑140 (MQ=255)
                                                                     aGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCa                                                          <  1:726005/140‑1 (MQ=255)
                                                                      gTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaa                                                         >  1:991680/1‑140 (MQ=255)
                                                                      gTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaa                                                         >  1:957167/1‑140 (MQ=255)
                                                                         cGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGa                                                      >  1:136127/1‑140 (MQ=255)
                                                                            aTGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACtt                                                   >  1:859530/1‑140 (MQ=255)
                                                                                 gttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCaaaa                                              >  2:736857/1‑140 (MQ=255)
                                                                                 gttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCaaaa                                              >  1:188477/1‑140 (MQ=255)
                                                                                 gttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCaaaa                                              >  1:302008/1‑140 (MQ=255)
                                                                                      ggTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGc                                         >  2:334167/1‑140 (MQ=255)
                                                                                        ttCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGt                                       >  1:92794/1‑140 (MQ=255)
                                                                                        ttCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑‑AAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGTTTCCAGCGCAAAGACTTCAAAACCTGCGTc                                      >  2:1068960/1‑140 (MQ=255)
                                                                                          ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCg                                     >  2:1051430/1‑140 (MQ=255)
                                                                                          ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCg                                     >  2:891307/1‑140 (MQ=255)
                                                                                              ttGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGa                                 <  1:549686/140‑1 (MQ=255)
                                                                                                ggAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGAcg                               <  2:748249/140‑1 (MQ=255)
                                                                                                  aaGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGAcgc                              >  2:661579/1‑140 (MQ=255)
                                                                                                  aaGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGAcgcg                             >  1:591523/1‑140 (MQ=255)
                                                                                                        gAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                       <  2:143392/140‑1 (MQ=255)
                                                                                                        gAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                       <  2:991680/140‑1 (MQ=255)
                                                                                                             ggCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGt                  <  2:349929/140‑1 (MQ=255)
                                                                                                              gCTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGtt                 >  2:160857/1‑140 (MQ=255)
                                                                                                               cTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTa                >  1:362761/1‑140 (MQ=255)
                                                                                                               cTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTa                >  2:524537/1‑140 (MQ=255)
                                                                                                               cTCAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTa                >  1:993580/1‑140 (MQ=255)
                                                                                                                 cAAGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaaa              <  2:141430/140‑1 (MQ=255)
                                                                                                                   aGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAAt            >  2:765182/1‑140 (MQ=255)
                                                                                                                   aGGCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAAt            <  2:957167/140‑1 (MQ=255)
                                                                                                                    ggCGAGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATc           <  2:993580/140‑1 (MQ=255)
                                                                                                                        aGCCTGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGt       >  2:669428/1‑140 (MQ=255)
                                                                                                                            tGGGACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGcc   >  1:942730/1‑140 (MQ=255)
                                                                                                                             gggACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCGAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa  <  2:880432/140‑1 (MQ=255)
                                                                                                                             gggACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa  <  2:362761/140‑1 (MQ=255)
                                                                                                                             gggACAAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa  <  2:942730/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑‑AAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCA  >  NC_000913/218698‑218961

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