Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,253,701 C→G 9.5% L600F (TTG→TTC)  metG ← methionine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,253,7010CG 9.5% 117.5 / 2.3 41L600F (TTG→TTC) metGmethionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base G (2/2);  total (22/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 = 1.00e+00

GGACACGCCAAAGCGCATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCTTGGCTTCG  >  NC_002947/1253554‑1253842
                                                                                                                                                   |                                                                                                                                             
ggACACGCCAAAGCGCATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTgcagc                                                                                                                                               >  4:193102/1‑148 (MQ=255)
         aaaGCGCATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGGAACTTGTCgg                                                                                                                                     >  6:195488/1‑149 (MQ=255)
              gcATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGGAACTTGTCGGCAccc                                                                                                                                >  2:212317/1‑149 (MQ=255)
              gcATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCAccc                                                                                                                                >  2:75815/1‑149 (MQ=255)
                aTCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGc                                                                                                                              >  5:21323/1‑149 (MQ=255)
                aTCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGc                                                                                                                              >  5:26058/1‑149 (MQ=255)
                  cTTACGGCGTTTGAGGTTGTCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCACCCTTGAACCCGGAGAACCCGTTGCGGCGATCCTAACCAATCTCCAGGCTCAGTTGCAGCAACTTGTCGGCACCCGCCa                                                                                                                            >  2:197254/1‑149 (MQ=255)
                              gAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCt                                                                                                                >  2:212050/1‑149 (MQ=255)
                               aGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCtt                                                                                                               <  6:298323/149‑1 (MQ=255)
                                          catcatGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAgcg                                                                                                    <  3:193102/149‑1 (MQ=255)
                                            tcatGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAgcgc                                                                                                   >  5:24421/1‑148 (MQ=255)
                                              atGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATCTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCcac                                                                                                <  1:325543/149‑1 (MQ=255)
                                                  tCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGAGGATCAGCGCCACACGc                                                                                            <  1:212050/149‑1 (MQ=255)
                                                              ccAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACgg                                                                                  >  7:165231/1‑147 (MQ=255)
                                                               cAGCTTGGCCGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCTTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCCGCGCCACACGCAGGTCGACGGCAg                                                                               >  8:200238/1‑149 (MQ=255)
                                                                    tGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACg                                                                                   >  7:99595/1‑140 (MQ=255)
                                                                    tGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACg                                                                                   <  8:99595/140‑1 (MQ=255)
                                                                         gggTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCAGCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATc                                                                     <  6:24421/149‑1 (MQ=255)
                                                                           gTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTc                                                                                                                     >  1:114340/1‑99 (MQ=255)
                                                                           gTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTc                                                                                                                     <  2:114340/99‑1 (MQ=255)
                                                                           gTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCg                                                                    >  5:42093/1‑148 (MQ=255)
                                                                            tCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGaa                                                                  >  4:64565/1‑149 (MQ=255)
                                                                                 gTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGGAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGa                                                             <  5:195488/149‑1 (MQ=255)
                                                                                     gCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGAtctc                                                         >  2:321645/1‑149 (MQ=255)
                                                                                       tGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTg                                                       >  7:146520/1‑149 (MQ=255)
                                                                                         aCTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGGAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGcc                                                     <  1:212317/149‑1 (MQ=255)
                                                                                         aCTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTg                                                       <  5:322291/147‑1 (MQ=255)
                                                                                            tGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGa                                                   >  5:130655/1‑148 (MQ=255)
                                                                                             gATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGa                                                   >  5:185951/1‑147 (MQ=255)
                                                                                                 ccGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCgg                                             <  3:64565/149‑1 (MQ=255)
                                                                                                  cGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCggg                                            <  1:75815/149‑1 (MQ=255)
                                                                                                       aaCACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCtt                                       <  6:21323/149‑1 (MQ=255)
                                                                                                       aaCACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCtt                                       <  6:26058/149‑1 (MQ=255)
                                                                                                       aaCACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCtt                                       >  2:241909/1‑149 (MQ=255)
                                                                                                       aaCACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGACGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCtt                                       >  1:248256/1‑149 (MQ=255)
                                                                                                        acacGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCtt                                       <  8:146520/148‑1 (MQ=255)
                                                                                                        acacGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGAACTCTGCCGACAGCGGGTCCtt                                       <  6:42093/148‑1 (MQ=255)
                                                                                                                             cAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCccgcc                 >  1:68856/1‑149 (MQ=255)
                                                                                                                                     aGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCtt         <  1:241909/149‑1 (MQ=255)
                                                                                                                                     aGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCtt         <  2:68856/149‑1 (MQ=255)
                                                                                                                                           aGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCTTGGCTTc   >  7:78371/1‑149 (MQ=255)
                                                                                                                                            gTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCTTGGCTTCg  <  8:165231/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GGACACGCCAAAGCGCATCTTCCGGGGTTTGAGGTTGGCTACCATCATGGTCAGGCGGCCTTCCAGCTTGGACGGGTCTGGGTAGGCTGACTTGATGCCGGAGAACACGTTGCGGCGCTCGTCACCAATGTCCAGGGTCAGTTGCAGCAACTTGTCGGCACCCGCCACGGCTTCAGCCTTGACGATCAGCGCCACACGCAGGTCGACGGCAGCGAAGGTATCGAACTCGATCTCTGCCGACAGCGGGTCCTTGGCCAGCTCGCCATTGCCCGCCGGGGCCTTGGCTTCG  >  NC_002947/1253554‑1253842

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