Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 568,752 A→G G71G (GGA→GGG ybcO → DUF1364 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273568,7520AG95.2% 59.8 / ‑3.4 21G71G (GGA→GGGybcODUF1364 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (5/15);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568680‑568838
                                                                        |                                                                                      
ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:529411/90‑1 (MQ=255)
ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:523618/90‑1 (MQ=255)
ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  1:335800/90‑1 (MQ=255)
ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  1:439145/90‑1 (MQ=255)
ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  1:607271/90‑1 (MQ=255)
      tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  2:352983/90‑1 (MQ=255)
                gCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  1:638199/90‑1 (MQ=255)
                   tGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCgaga                                                    <  2:112793/90‑1 (MQ=255)
                       ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                >  1:188459/1‑90 (MQ=255)
                       ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                >  1:188460/1‑90 (MQ=255)
                       ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                <  1:127779/90‑1 (MQ=255)
                                 cgacgaAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTg                                      >  2:163030/1‑90 (MQ=255)
                                       aaTCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGaa                                >  1:79320/1‑90 (MQ=255)
                                           gACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAagag                            <  2:643751/90‑1 (MQ=255)
                                             ccgccgCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAggg                          >  2:572686/1‑90 (MQ=255)
                                                        aTTTTGTTGACGCTGGGTATGCAAAAGAATGTGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAgg               <  2:148731/90‑1 (MQ=255)
                                                           ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  2:687856/90‑1 (MQ=255)
                                                             gttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGa          <  1:532446/90‑1 (MQ=255)
                                                              ttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGaa         <  2:305606/90‑1 (MQ=255)
                                                                 aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGGTATTAAGGCGGGAATAc      >  1:361948/1‑90 (MQ=255)
                                                                     tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  2:315691/90‑1 (MQ=255)
                                                                        |                                                                                      
CCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568680‑568838

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