Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,390,754 G→T A252A (GCC→GCA rsgA ← ribosome small subunit‑dependent GTPase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,7540GT95.5% 72.0 / NA 22A252A (GCC→GCArsgAribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (14/7);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTC  >  NC_000913/4390693‑4390813
                                                              |                                                           
aCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAgctgcg                                                          >  1:2893389/1‑65 (MQ=255)
  tCACGCACTCCTGGGGAGTCAATCACTTCCCCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGt                                                        <  1:847715/65‑1 (MQ=255)
  tCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGt                                                        <  1:1041365/65‑1 (MQ=255)
  tCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGt                                                        <  1:886054/65‑1 (MQ=255)
  tCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCCGGAAGTGATACAGCCGAGCTGCGGt                                                        >  1:3184098/1‑65 (MQ=255)
          tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGgtgt                                                <  1:2771027/65‑1 (MQ=255)
              ggggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTg                                            >  1:1453432/1‑65 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGAccc                                        >  1:1487032/1‑64 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGAccc                                        >  1:2688652/1‑64 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCg                                        >  1:1796862/1‑65 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCg                                        >  1:817074/1‑65 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCg                                        >  1:1988799/1‑65 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCg                                        >  1:291383/1‑65 (MQ=255)
                  aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGAGTGCTGACCg                                        >  1:1633186/1‑65 (MQ=255)
                                  gccgTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCag                        >  1:138826/1‑65 (MQ=255)
                                   ccgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCaga                       <  1:1696188/65‑1 (MQ=255)
                                         gggAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                 <  1:2032840/65‑1 (MQ=255)
                                          ggAAGTTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                 <  1:2660549/65‑1 (MQ=255)
                                              gTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGCTgt            >  1:2997655/1‑64 (MQ=255)
                                                  aTACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCa         <  1:491889/65‑1 (MQ=255)
                                                   tACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg        <  1:2270926/65‑1 (MQ=255)
                                                    agtgCCGAGCTGCGGTCGTGGTGGGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAgg       <  1:3056225/62‑1 (MQ=255)
                                                    aCAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAgg       >  1:2862421/1‑65 (MQ=255)
                                                    aCAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAgg       >  1:1973862/1‑65 (MQ=255)
                                                    aCAGCCGAGCTGCGGTCGTGGGGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGCCAgg       >  1:722130/1‑65 (MQ=255)
                                                       gCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtc    >  1:2958795/1‑65 (MQ=255)
                                                         ggAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtctc  <  1:914625/64‑1 (MQ=255)
                                                         cGAGCTGCGGTCGTGGTGTGCTGGCCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtctc  <  1:1824874/65‑1 (MQ=255)
                                                         cGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtctc  >  1:945671/1‑65 (MQ=255)
                                                              |                                                           
ACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTC  >  NC_000913/4390693‑4390813

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