Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,062,210 T→C G390G (GGT→GGC agp → bifunctional glucose‑1‑phosphatase/inositol phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,062,2100TC88.2% 37.5 / ‑1.0 17G390G (GGT→GGCagpbifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (11/4);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAT  >  NZ_CP009273/1062132‑1062280
                                                                              |                                                                      
gagTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGa                                                             >  1:388186/1‑90 (MQ=255)
  gTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACg                                                           >  2:353322/1‑90 (MQ=255)
    gCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGc                                                                     >  1:456444/1‑78 (MQ=255)
    gCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGc                                                                     <  2:456444/78‑1 (MQ=255)
    gCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCt                                                         >  2:68338/1‑90 (MQ=255)
           aGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGtt                                                  <  2:815524/90‑1 (MQ=255)
                    aTGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGa                                         >  1:302140/1‑90 (MQ=255)
                      gCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATg                                       >  2:40838/1‑90 (MQ=255)
                      gCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATg                                       >  1:162830/1‑90 (MQ=255)
                                         aGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCgtgt                    <  2:156410/90‑1 (MQ=255)
                                               cTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaatg              >  2:361228/1‑90 (MQ=255)
                                               cTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaatg              >  2:620838/1‑90 (MQ=255)
                                                      gCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGt       <  2:565442/90‑1 (MQ=255)
                                                        gtgtGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                   <  2:543588/60‑1 (MQ=255)
                                                        gtgtGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                   >  1:543588/1‑60 (MQ=255)
                                                        gtgtGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGa     >  1:668707/1‑90 (MQ=255)
                                                           tgACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAt  >  2:668681/1‑90 (MQ=255)
                                                                              |                                                                      
GAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAT  >  NZ_CP009273/1062132‑1062280

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: