Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,222 A→G V106A (GTG→GCG)  rbbA ← ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,2220AG96.8% 100.9 / NA 31V106A (GTG→GCG) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (12/18);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTGCC  >  NC_000913/3629089‑3629340
                                                                                                                                     |                                                                                                                      
gccgGAGAGTCTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACa                                                                                                                  <  2:772087/140‑1 (MQ=255)
 ccgGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACa                                                                                                                  >  1:727568/1‑139 (MQ=255)
 ccgGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACa                                                                                                                  <  2:727568/139‑1 (MQ=255)
   ggAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAgg                                                                                                               <  1:663820/140‑1 (MQ=255)
    gagaGTTTCCCTGCCAGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGgt                                                                                                              <  2:770611/140‑1 (MQ=255)
                gCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCAGACACCGACAGGGTGTGGTAGAGGtt                                                                                                  >  1:233072/1‑140 (MQ=255)
                                aCGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGc                                                                                  <  1:11385/140‑1 (MQ=255)
                                aCGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATACGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGc                                                                                  <  1:325133/140‑1 (MQ=255)
                                   gTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGAGGTAGAGGTTTTTGCCCAGCCCCTGCGGc                                                                               >  1:293472/1‑140 (MQ=255)
                                              gTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGAt                                                                    >  2:89206/1‑140 (MQ=255)
                                                 cTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATgcg                                                                 <  2:464187/140‑1 (MQ=255)
                                                   ggTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGcgcg                                                               >  1:621392/1‑140 (MQ=255)
                                                     tCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCGCAGCCCCTGCGGCCTCCAGGCGGTGCGCggg                                                             >  2:375376/1‑140 (MQ=255)
                                                              cATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGTTTTTTGCCTAGCCCCTATGGCATCCTGGCGTTGCGCGGGCAGACGTcg                                                    >  2:398578/1‑140 (MQ=255)
                                                              cATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTcg                                                    >  1:730076/1‑140 (MQ=255)
                                                              cATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTcg                                                    >  2:609/1‑140 (MQ=255)
                                                                   aTTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCgcggcg                                               <  2:730076/140‑1 (MQ=255)
                                                                   aTTCGCACTTCCCGCACCGCTTAGTAGTGCCCGAACATGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCgcggcg                                               <  1:375376/140‑1 (MQ=255)
                                                                      cGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATg                                            >  1:689398/1‑140 (MQ=255)
                                                                          cTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTc                                        <  2:424840/140‑1 (MQ=255)
                                                                                     gctTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCat                             >  1:156323/1‑140 (MQ=255)
                                                                                         tGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATcg                         >  2:526101/1‑140 (MQ=255)
                                                                                              tGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCgccgcc                    <  2:186030/140‑1 (MQ=255)
                                                                                               gACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGccc                   <  2:689398/140‑1 (MQ=255)
                                                                                                 ccGAACAGGCGAGCGCAAAAATCGACGTTCTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAg                 <  1:394272/140‑1 (MQ=255)
                                                                                                      cAGGCGAGCGAAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCa            <  1:526101/140‑1 (MQ=255)
                                                                                                         gcgagcgaAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATca         >  1:538534/1‑140 (MQ=255)
                                                                                                              cgaAAAAATCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTg    <  2:621392/140‑1 (MQ=255)
                                                                                                                aaaaaaTCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTGcc  <  2:156323/140‑1 (MQ=255)
                                                                                                                aaaaaaTCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTGcc  <  1:737656/140‑1 (MQ=255)
                                                                                                                aaaaaaTCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTGcc  <  1:525248/140‑1 (MQ=255)
                                                                                                                aaaaaaGCGACGTTTTCATACGCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCACGCATATCGCCGCCCAGCACCATCACATTGcc  <  1:641712/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                      
GCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCTCATTGATTCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCCAGCACCATCACATTGCC  >  NC_000913/3629089‑3629340

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