Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,941,604 A→C G10G (GGT→GGG ruvC ← crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,941,6040AC93.3% 40.7 / ‑4.8 15G10G (GGT→GGGruvCcrossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (8/6);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941517‑1941686
                                                                                       |                                                                                  
atcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAccc                                                                                  >  2:168211/1‑90 (MQ=255)
atcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAccc                                                                                  >  1:622315/1‑90 (MQ=255)
 tcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 <  2:378939/90‑1 (MQ=255)
       ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCa                                                                           >  1:300295/1‑90 (MQ=255)
       ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCa                                                                           >  2:481012/1‑90 (MQ=255)
        tGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCaa                                                                          <  2:26118/90‑1 (MQ=255)
              ggATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCga                                                                    >  2:252877/1‑90 (MQ=255)
                              cAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAc                                                    <  2:300295/90‑1 (MQ=255)
                                   aGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtc                                               <  2:152998/90‑1 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  1:807001/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  1:429911/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:487697/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:501165/1‑90 (MQ=255)
                                                acctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGt                                  <  1:579278/90‑1 (MQ=255)
                                                   tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTa                               >  1:568072/1‑90 (MQ=255)
                                                         gCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGc                         >  1:149058/1‑90 (MQ=255)
                                                                               aCGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAc   <  2:738307/90‑1 (MQ=255)
                                                                                cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  1:481012/90‑1 (MQ=255)
                                                                                       |                                                                                  
ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941517‑1941686

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