Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,739,096 T→C E61G (GAA→GGA)  rimM ← ribosome maturation factor RimM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,739,0960TC94.7% 50.7 / ‑3.2 19E61G (GAA→GGA) rimMribosome maturation factor RimM
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (11/7);  total (12/7)
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.77e-01

TCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGG  >  NZ_CP009273/2739016‑2739177
                                                                                |                                                                                 
tCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGa                                                                          >  2:167850/1‑90 (MQ=255)
tCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGa                                                                          >  2:22256/1‑90 (MQ=255)
                 gcATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGAccc                                                         >  2:157893/1‑90 (MQ=255)
                            gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                              <  1:439018/90‑1 (MQ=255)
                            gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                              <  2:42129/90‑1 (MQ=255)
                                   ccTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAg                                       >  1:120852/1‑90 (MQ=255)
                                     tttCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                              >  1:318961/1‑81 (MQ=255)
                                     tttCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                              <  2:318961/81‑1 (MQ=255)
                                           cTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATa                               <  2:29511/90‑1 (MQ=255)
                                                    cATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAAt                      >  2:452845/1‑90 (MQ=255)
                                                        tCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTc                            >  1:309853/1‑80 (MQ=255)
                                                        tCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTc                            <  2:309853/80‑1 (MQ=255)
                                                            gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGc                                                       <  2:74098/49‑1 (MQ=255)
                                                            gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGc                                                       >  1:74098/1‑49 (MQ=255)
                                                            gATTGTGGTGCTTCCAGCTTCCCAGCTGGACGTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg              >  2:446591/1‑90 (MQ=255)
                                                                gtggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTc          >  1:272185/1‑90 (MQ=255)
                                                                gtggtgCTTCCAGCTTCACAGCTGGACTTGCTGCCACTGACCAGCGTGCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTc          >  1:525601/1‑90 (MQ=255)
                                                                 tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCt         >  2:432327/1‑90 (MQ=255)
                                                                      cTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGt    <  1:452845/90‑1 (MQ=255)
                                                                        tccagcttcccagctGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTgg  >  2:59408/1‑90 (MQ=255)
                                                                                |                                                                                 
TCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGG  >  NZ_CP009273/2739016‑2739177

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