Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,813,272 A→T 18.1% N46I (AAC→ATC)  argS → arginine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,813,2720AT18.1% 49.6 / 4.7 22N46I (AAC→ATC) argSarginine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base T (2/2);  total (12/10)
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

CCATGAAAGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACG  >  NC_002947/5813134‑5813418
                                                                                                                                          |                                                                                                                                                  
ccATGAAAGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCatga                                                                                                                                         >  2:309202/1‑150 (MQ=255)
     aaaGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTg                                                                                                                                    <  2:250956/150‑1 (MQ=255)
        gACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCatga                                                                                                                                         >  1:201761/1‑142 (MQ=255)
        gACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCatga                                                                                                                                         <  2:201761/142‑1 (MQ=255)
                     aGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCa                                                                                                                    <  1:105751/150‑1 (MQ=255)
                                     ccTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCTAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCt                                                                                                    >  1:503358/1‑150 (MQ=255)
                                     ccTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCATCATGGCGATGTTGCTg                                                                                                                                    >  2:487076/1‑118 (MQ=255)
                                     ccTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCATCATGGCGATGTTGCTg                                                                                                                                    <  1:487076/118‑1 (MQ=255)
                                       tCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCg                                                                                                  >  1:51049/1‑150 (MQ=255)
                                        cACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGc                                                                                                 <  1:407799/150‑1 (MQ=255)
                                         aCCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGGTGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGcc                                                                                                >  1:436607/1‑150 (MQ=255)
                                                   tCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGa                                                                                      <  2:51049/150‑1 (MQ=255)
                                                                       gAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGa                                                                                      >  1:341476/1‑130 (MQ=255)
                                                                       gAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGa                                                                                      <  2:341476/130‑1 (MQ=255)
                                                                                             aGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCgaga                                            >  2:550394/1‑150 (MQ=255)
                                                                                                               aCAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTcc                          <  1:133814/150‑1 (MQ=255)
                                                                                                                  aGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCa                       >  2:368100/1‑150 (MQ=255)
                                                                                                                            cGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCa             <  2:503358/150‑1 (MQ=255)
                                                                                                                                    ccAGCATCATGGCGATGTTGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCa                                       >  2:61146/1‑116 (MQ=255)
                                                                                                                                    ccAGCATCATGGCGATGTTGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCa                                       <  1:61146/116‑1 (MQ=255)
                                                                                                                                    ccAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCg     >  1:392234/1‑150 (MQ=255)
                                                                                                                                       gCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCACCGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACg  >  1:477860/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CCATGAAAGACACCATTCGCCAGCTGATCCAGCAAGCCCTCACCCAACTCGTCACCGACGGTGTGCTGCCTGAAGGGCTGTCGCCGGCGATCCAGGTGGAAAACGCCCGGGACAAGACCCACGGCGACTTCGCCAGCAACATCGCCATGATGCTGGCCAAGCCGGCCGGCATGAAGCCGCGCGACCTCGCCGAAAAACTGATCAACGCCCTGCCGGCCAGCGCCGACATCAGCAAGGTCGAGATTGCAGGCCCCGGCTTCCTCAACTTCTTCCAGAACACCGACG  >  NC_002947/5813134‑5813418

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