Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,947,467 A→C 16.7% intergenic (+4/‑30) metZ → / → metW tRNA‑Met/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC16.7% 62.5 / 9.2 30intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (9/16);  new base C (2/3);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTT  >  NC_000913/2947319‑2947552
                                                                                                                                                    |                                                                                     
aaaTCGGTGATCGGTGTTAAAAAAAGGTTGCATGAAAACGCGAGCGTAGTATAGTGCGCATCGACTGAGGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTGGTCGGTTCAATTCCGGCCCCCGCAACCAATTa                                                                                       <  2:347323/149‑1 (MQ=21)
  aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                     <  6:328501/149‑1 (MQ=255)
     cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGCGCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                  <  1:75240/149‑1 (MQ=255)
       tGAGCGGGTTAAAAAAAAGGGTGCATGAAAACGCGAGCGGAGTAGAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                               <  6:172042/149‑1 (MQ=18)
        gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                               <  4:42870/149‑1 (MQ=255)
         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                               <  6:164290/148‑1 (MQ=255)
            ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                           <  5:43757/149‑1 (MQ=255)
                  aaaaaaaGGTGGCAGGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                      <  8:204437/148‑1 (MQ=255)
                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                               >  1:51100/1‑137 (MQ=255)
                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                               <  2:51100/137‑1 (MQ=255)
                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                   <  4:149022/149‑1 (MQ=255)
                                       gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                <  6:135239/149‑1 (MQ=255)
                                       gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                              <  1:320485/118‑1 (MQ=255)
                                       gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                              >  2:320485/1‑118 (MQ=255)
                                        cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                            <  2:396388/136‑1 (MQ=255)
                                        cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                            >  1:396388/1‑136 (MQ=255)
                                           gtggAGTTTAGTGAGCATCCCCTGACGCTGGTTGCAGCAGCCTGTTAGCTCGTCGGGATCATAACCCGAATGTCTTATGTTCAAATCATTCCCCCGCAACCAATTAAACTTTGATGAAGTAACGCAGTCCGGTGACGCGGGGTGGagca                                            <  4:372068/147‑1 (MQ=9)
                                                      tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                              >  8:275063/1‑103 (MQ=255)
                                                      tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                              <  7:275063/103‑1 (MQ=255)
                                                                                               tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                            <  7:310194/97‑1 (MQ=255)
                                                                                               tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                            >  8:310194/1‑97 (MQ=255)
                                                                                                     tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                 >  8:180872/1‑86 (MQ=255)
                                                                                                       aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                             <  5:158205/56‑1 (MQ=255)
                                                                                                       aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                             >  6:158205/1‑56 (MQ=255)
                                                                                                          aCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAAGGAAGAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                     >  8:5924/1‑93 (MQ=14)
                                                                                                                                   ggCCCCCGCAACCAATTAAAATTTGATGAATTAAAGCAGTACGGTGACGCGGGGTGGa                                               >  4:271047/1‑58 (MQ=255)
                                                                                                                                       cccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGaa             >  4:119705/1‑88 (MQ=255)
                                                                                                                                       cccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGaa             <  3:119705/88‑1 (MQ=255)
                                                                                                                                                 aTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAt                     <  4:365267/70‑1 (MQ=255)
                                                                                                                                                 aTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGtt  <  2:136649/89‑1 (MQ=255)
                                                                                                                                                   tAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg            <  6:257792/77‑1 (MQ=40)
                                                                                                                                                   tAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg            >  5:257792/1‑77 (MQ=40)
                                                                                                                                                    |                                                                                     
AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTT  >  NC_000913/2947319‑2947552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: