Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,441,796 A→C W42G (TGG→GGG)  prmB ← 50S ribosomal protein L3 N(5)‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,441,7960AC96.9% 103.7 / ‑4.2 32W42G (TGG→GGG) prmB50S ribosomal protein L3 N(5)‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (15/16);  total (15/17)
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.51e-01

CGGCGGATCACGCGTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCC  >  NZ_CP009273/2441672‑2441929
                                                                                                                            |                                                                                                                                     
cggcggGTCACGCGTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGCAGGCCACCCCAGCTCTACGGCTTCATCCCCCCGGCTATCGGTAc                                                                                                                         >  2:1107609/1‑139 (MQ=255)
 ggcggATCACGCGTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTAcc                                                                                                                        <  2:1166680/139‑1 (MQ=255)
     gATCACGCGTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGa                                                                                                                    <  2:945790/139‑1 (MQ=255)
           gcgTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTCCATCCCCCGGGTTATCGGTACCGTGACCGTAc                                                                                                              <  1:1107609/139‑1 (MQ=255)
                caacaaTACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCgg                                                                                                                                    <  2:479413/112‑1 (MQ=255)
                caacaaTACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCgg                                                                                                                                    >  1:479413/1‑112 (MQ=255)
                caacaaTACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGa                                                                                                          >  2:459811/1‑138 (MQ=255)
                         ggTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCcgc                                                                                                <  1:93996/139‑1 (MQ=255)
                             tttttCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCt                                                                                             >  1:851003/1‑138 (MQ=255)
                                  cGCTGGAGGTCAGACGCGCGGTGCGCATATCTGCCGGAATATCCAGCGGCAGGTAGAGCGAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCg                                                                                       <  1:939799/138‑1 (MQ=255)
                                     tGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCt                                                                                    >  1:244196/1‑139 (MQ=255)
                                             aGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGAc                                                                             <  1:363080/138‑1 (MQ=255)
                                             aGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGAc                                                                             <  2:435611/138‑1 (MQ=255)
                                             aGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTGATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGAc                                                                             <  2:1059161/138‑1 (MQ=255)
                                                                 ttCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATg                                                        >  1:714329/1‑139 (MQ=255)
                                                                 ttCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATg                                                        >  2:386840/1‑139 (MQ=255)
                                                                   ccGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGt                                                      >  2:503525/1‑139 (MQ=255)
                                                                       aaTATCCAGCGGCAGGTACAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATAGTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCCGAATGGTTTGc                                                  <  2:966795/139‑1 (MQ=255)
                                                                         tatCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAg                                                <  2:564769/139‑1 (MQ=255)
                                                                                 ggCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTa                                        >  1:192625/1‑139 (MQ=255)
                                                                                  gCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTAc                                       <  2:851003/139‑1 (MQ=255)
                                                                                               aGGCAACACCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAAc                           >  1:463319/1‑138 (MQ=255)
                                                                                                    acacCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAAt                     >  1:838485/1‑139 (MQ=255)
                                                                                                     cacCAGCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAAtt                    >  1:461867/1‑139 (MQ=255)
                                                                                                          gCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATc               >  1:898629/1‑139 (MQ=255)
                                                                                                          gCTGTACGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATc               <  2:463319/139‑1 (MQ=255)
                                                                                                              tACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCa             <  2:635236/137‑1 (MQ=255)
                                                                                                               aCGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTa                 <  1:626308/132‑1 (MQ=255)
                                                                                                               aCGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTa                 >  2:626308/1‑132 (MQ=255)
                                                                                                                cGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTa         >  2:112377/1‑139 (MQ=255)
                                                                                                                cGGCTTCATCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTa         <  2:192625/139‑1 (MQ=255)
                                                                                                                       aTCCCCCGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATtcctcc  <  1:459811/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                     
CGGCGGATCACGCGTTCAACAATACGGTGTTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCGCGCTGAAGCGGCTCACCGACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCC  >  NZ_CP009273/2441672‑2441929

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