Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,810,641 C→G 11.4% C335W (TGC→TGG thrS → threonine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,810,6410CG11.4% 103.5 / 4.0 35C335W (TGC→TGGthrSthreonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (13/18);  new base G (2/2);  total (15/20)
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.93e-01

CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGC  >  NC_002947/2810500‑2810784
                                                                                                                                             |                                                                                                                                               
cTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCa                                                                                                                                         <  1:42381/150‑1 (MQ=255)
 cactaccaaaTCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAg                                                                                                                                        <  2:258038/143‑1 (MQ=255)
   ccAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgt                                                                                                                                      <  2:732195/150‑1 (MQ=255)
   ccAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgt                                                                                                                                      >  1:60437/1‑150 (MQ=255)
     aGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgtgt                                                                                                                                    <  1:93442/150‑1 (MQ=255)
                             tCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACTGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAcc                                                                                                            >  2:740552/1‑150 (MQ=255)
                                   gCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAcc                                                                                                      >  1:736072/1‑150 (MQ=255)
                                     aTCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTg                                                                                                    >  2:1067843/1‑150 (MQ=255)
                                       ccTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc                                                                                                  <  2:1017228/150‑1 (MQ=255)
                                            gggAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGc                                                                                             <  1:398737/150‑1 (MQ=255)
                                                       ggCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAg                                                                                  >  1:102586/1‑150 (MQ=255)
                                                        gCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGt                                                                                 <  2:904734/150‑1 (MQ=255)
                                                          cACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAAATGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTc                                                                               >  1:145943/1‑150 (MQ=255)
                                                            cTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAc                                                                                                             <  2:861248/118‑1 (MQ=255)
                                                            cTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAc                                                                                                             >  1:861248/1‑118 (MQ=255)
                                                                 ccAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc                                                                                                  >  2:147359/1‑124 (MQ=255)
                                                                 ccAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGAATACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc                                                                                                  <  1:147359/124‑1 (MQ=255)
                                                                      tATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCAc                                                                   >  1:611918/1‑150 (MQ=255)
                                                                                      tCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCg                                                   <  1:513310/150‑1 (MQ=255)
                                                                                             ttCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCt                                            >  1:445672/1‑150 (MQ=255)
                                                                                                gtcgGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCa                                         >  2:1176076/1‑150 (MQ=255)
                                                                                                         ccGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGcatcat                                <  1:162366/150‑1 (MQ=255)
                                                                                                                         aaaCCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTc                <  1:1176076/150‑1 (MQ=255)
                                                                                                                         aaaCCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTc                <  2:736072/150‑1 (MQ=255)
                                                                                                                             cGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTAc            >  1:1134989/1‑150 (MQ=255)
                                                                                                                              gATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACa           <  2:256982/150‑1 (MQ=255)
                                                                                                                               aTGAACTGCCCGGGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAg          <  2:145943/150‑1 (MQ=255)
                                                                                                                                 gAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGa        <  2:1134989/150‑1 (MQ=255)
                                                                                                                                  accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGgcgc                                               >  1:491938/3‑110 (MQ=255)
                                                                                                                                  accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGgcgc                                               <  2:491938/108‑1 (MQ=255)
                                                                                                                                  accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGAc       >  2:458485/3‑150 (MQ=255)
                                                                                                                                     tGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGAt    <  2:611918/150‑1 (MQ=255)
                                                                                                                                       cccGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAc                                                                                                       >  2:658149/1‑49 (MQ=255)
                                                                                                                                       cccGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAc                                                                                                       <  1:658149/49‑1 (MQ=255)
                                                                                                                                       cacGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATgc  <  1:458485/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGC  >  NC_002947/2810500‑2810784

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