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,2100TC94.1% 49.6 / ‑5.9 17G390G (GGT→GGCagpbifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/5);  minor base A (0/1);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAG  >  NZ_CP009273/1062134‑1062285
                                                                            |                                                                           
gTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACg                                                                >  1:245991/1‑90 (MQ=255)
gTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACg                                                                >  2:721641/1‑90 (MQ=255)
          gTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGttt                                                      >  1:511939/1‑90 (MQ=255)
          gTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGttt                                                      >  2:264811/1‑90 (MQ=255)
                  aTGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGa                                              <  1:275152/90‑1 (MQ=255)
                  aTGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGa                                              >  1:346807/1‑90 (MQ=255)
                    gCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATg                                            >  1:373356/1‑90 (MQ=255)
                        aTGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                        <  2:245991/90‑1 (MQ=255)
                               aaCCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGa                                 <  2:244409/90‑1 (MQ=255)
                                  ccTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAg                              >  2:842297/1‑90 (MQ=255)
                                    tGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCg                            >  2:686749/1‑90 (MQ=255)
                                     gCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCgt                           >  2:167906/1‑90 (MQ=255)
                                                cggaGCGTGTCACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCTGATGAATACGTTTGATAGCGTGTTGAATGAAg                <  2:25927/87‑1 (MQ=255)
                                                      gtgtGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAg                              >  2:156550/1‑70 (MQ=255)
                                                      gtgtGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAg                              <  1:156550/70‑1 (MQ=255)
                                                           ccGCTGGAATTATGGGGATGCCCGATAGCCGCTGATGGTTTCTGCCCCATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATaa     <  2:772456/89‑1 (MQ=255)
                                                              cTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAg  >  1:306725/1‑90 (MQ=255)
                                                                            |                                                                           
GTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAG  >  NZ_CP009273/1062134‑1062285

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