Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,326,882 A→G intergenic (‑72/‑187) rplU ← / → ispB 50S ribosomal protein L21/octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,326,8820AG97.6% 141.7 / ‑5.6 42intergenic (‑72/‑187)rplU/ispB50S ribosomal protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/20);  minor base C (0/1);  total (21/21)
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.90e-01

GCTTACTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAT  >  NZ_CP009273/3326766‑3327008
                                                                                                                    |                                                                                                                              
gCTTACTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCt                                                                                                          >  1:68224/1‑139 (MQ=255)
gCTTACTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCt                                                                                                          >  1:432099/1‑139 (MQ=255)
     cTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGAc                                                                                                     <  1:132352/139‑1 (MQ=255)
     cTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGAc                                                                                                     <  1:422798/139‑1 (MQ=255)
            ttgttTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCt                                                                                              <  2:221946/139‑1 (MQ=255)
               tttACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTAc                                                                                            <  2:67719/138‑1 (MQ=255)
                         cTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATacac                                                                                 >  2:323393/1‑139 (MQ=255)
                            tGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACNNg                                                                                >  2:133638/1‑137 (MQ=255)
                                      cgTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGAc                                                                                                     >  1:125517/1‑106 (MQ=255)
                                      cgTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGAc                                                                                                     <  2:125517/106‑1 (MQ=255)
                                      cgTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaaaaa                                                                     >  1:254691/1‑138 (MQ=255)
                                        tACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAacac                                                                  <  1:285815/139‑1 (MQ=255)
                                                     ccGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGt                                                      >  1:404566/1‑138 (MQ=255)
                                                      cGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGt                                                           <  1:33547/132‑1 (MQ=255)
                                                      cGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGCCAATACACGAAGAAAAAAAACACAACTTGt                                                           >  2:33547/1‑132 (MQ=255)
                                                       gCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAAc                                                   >  1:206757/1‑139 (MQ=255)
                                                            cgcgCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTa                                              >  2:24232/1‑139 (MQ=255)
                                                                   acCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGc                                       >  2:267271/1‑139 (MQ=255)
                                                                        ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                    <  2:432099/137‑1 (MQ=255)
                                                                        ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                    <  2:254691/137‑1 (MQ=255)
                                                                        ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                    <  2:68224/137‑1 (MQ=255)
                                                                               aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAAc                                                               >  2:165067/1‑103 (MQ=255)
                                                                               aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAAc                                                               <  1:165067/103‑1 (MQ=255)
                                                                                     agTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                     <  1:399646/139‑1 (MQ=255)
                                                                                     agTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                     <  2:404566/139‑1 (MQ=255)
                                                                                     agTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                     <  1:24232/139‑1 (MQ=255)
                                                                                      gTGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCa                    <  1:323393/139‑1 (MQ=255)
                                                                                       tGCGCTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc                   >  2:257071/1‑139 (MQ=255)
                                                                                         cgcTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                 >  1:88148/1‑139 (MQ=255)
                                                                                         cgcTATAAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                 >  2:301031/1‑139 (MQ=255)
                                                                                            tataAATATTCACAATAGGCCGCGCATATTCCACAAAACGCACCCCTTTGACAAGTCCTACAGTCATTACTCAAAAAAAAAAAACACAACTTGTACGTTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                <  1:133638/137‑1 (MQ=255)
                                                                                              taAATATTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat              <  2:206757/137‑1 (MQ=255)
                                                                                                   atacacaATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                <  1:285434/127‑1 (MQ=255)
                                                                                                   atTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAACAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                >  2:285434/1‑130 (MQ=255)
                                                                                                   atTCACAATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt        >  1:17950/1‑138 (MQ=255)
                                                                                                      cacaATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTa                                                          >  1:311746/1‑85 (MQ=255)
                                                                                                      cacaATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTa                                                          <  2:311746/85‑1 (MQ=255)
                                                                                                      cacaATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc     >  1:355613/1‑138 (MQ=255)
                                                                                                      cacaATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc     >  1:375694/1‑138 (MQ=255)
                                                                                                        caATAGGGCGCGGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAt  <  2:336518/139‑1 (MQ=255)
                                                                                                                cgcgGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                     >  1:194512/1‑112 (MQ=255)
                                                                                                                cgcgGATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                     <  2:194512/112‑1 (MQ=255)
                                                                                                                    |                                                                                                                              
GCTTACTCGGTGTTGTTTACCACCACTTTGGAAAACCGCGTACATATAAAACTCCGCTTCCGCGCACACCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAT  >  NZ_CP009273/3326766‑3327008

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