Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 687,482 C→T A106T (GCG→ACG)  ybeY ← rRNA maturation RNase YbeY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,4820CT80.7% 72.6 / 10.2 31A106T (GCG→ACG) ybeYrRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base C (6/0);  new base T (17/8);  total (23/8)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGTCGGCTTATCCT  >  NZ_CP009273/687354‑687616
                                                                                                                                |                                                                                                                                      
aTCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGt                                                                                                                              >  1:763963/1‑139 (MQ=255)
     cATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGcc                                                                                                                         >  1:43702/1‑139 (MQ=255)
     cATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGcc                                                                                                                         >  2:762594/1‑139 (MQ=255)
               agcaagcaTAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCt                                                                                                                 >  1:78877/1‑137 (MQ=255)
                  aagcaTAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGct                                                                                                            <  1:824245/139‑1 (MQ=255)
                      aTAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGcttcct                                                                                                        <  2:165930/139‑1 (MQ=255)
                        aaTCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGcttccttc                                                                                                      >  2:826820/1‑139 (MQ=255)
                               gTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAacca                                                                                               >  2:167860/1‑139 (MQ=255)
                                tttCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAaccac                                                                                              <  1:497955/139‑1 (MQ=255)
                                tttCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAaccac                                                                                              >  2:823254/1‑139 (MQ=255)
                                     aGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGt                                                                                                                              >  1:181347/1‑102 (MQ=255)
                                     aGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGt                                                                                                                              <  2:181347/102‑1 (MQ=255)
                                     aGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGAc                                                                                         >  1:664659/1‑139 (MQ=255)
                                     aGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGAc                                                                                         >  2:19257/1‑139 (MQ=255)
                                     aGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGAc                                                                                         >  2:192778/1‑139 (MQ=255)
                                       ggCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACgg                                                                                       <  1:259626/139‑1 (MQ=255)
                                         cTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCa                                                                                     <  1:826820/139‑1 (MQ=255)
                                                                       tgtgATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTc                                                       >  2:731177/1‑139 (MQ=255)
                                                                          gATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCAt                                                    >  1:143210/1‑139 (MQ=255)
                                                                                     aaCAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCa                                         >  1:702559/1‑139 (MQ=255)
                                                                                           tGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCaa                                   <  2:43702/139‑1 (MQ=255)
                                                                                                   gCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCCTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAg                           >  2:715843/1‑139 (MQ=255)
                                                                                                    ccGGGCACCACCATATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGGAg                           >  2:855128/1‑138 (MQ=255)
                                                                                                            caccaTATGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACg                  >  2:450081/1‑139 (MQ=255)
                                                                                                              ccaTATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGtt                >  1:759722/1‑139 (MQ=255)
                                                                                                                atatGCGCCCAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTgg              <  2:232277/139‑1 (MQ=255)
                                                                                                                 tatGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGt             >  1:862927/1‑139 (MQ=255)
                                                                                                                      gcCCAGTGCGTCTCCAGTGGTTTGCCATGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGTCGGCt        >  1:830471/1‑139 (MQ=255)
                                                                                                                         cAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGTCGGCTTAt     >  2:193538/1‑139 (MQ=255)
                                                                                                                         cAGTGCGTCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGTCGGCTTAt     >  1:562159/1‑139 (MQ=255)
                                                                                                                            tGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGCTCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGCAGGAGAGCACGTTGGTCGGCTTATCCt  >  1:713141/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                      
ATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGGGAAGGAGAGCACGTTGGTCGGCTTATCCT  >  NZ_CP009273/687354‑687616

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