Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,568,321 G→A G336G (GGC→GGT glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,568,3210GA97.8% 149.7 / ‑4.2 45G336G (GGC→GGTglgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/22);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

ATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGA  >  NC_000913/3568244‑3568395
                                                                             |                                                                          
aTGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAAccaccg                                                                         <  1:687566/81‑1 (MQ=255)
 tGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCgg                                                                        <  2:687567/81‑1 (MQ=255)
 tGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCgg                                                                        <  2:687568/81‑1 (MQ=255)
   ttGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGa                                                                       <  2:687569/80‑1 (MQ=255)
    tGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGaaa                                                                     >  2:687627/1‑81 (MQ=255)
      cAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAAcc                                                                   <  2:687570/81‑1 (MQ=255)
      cAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAAc                                                                    >  1:687621/1‑80 (MQ=255)
        gaatgaatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAg                                                                 <  2:687571/81‑1 (MQ=255)
            gaatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAg                                                             >  2:687641/1‑81 (MQ=255)
             aatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAg                                                             <  2:687572/80‑1 (MQ=255)
             aatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGt                                                            >  2:687656/1‑81 (MQ=255)
              atTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGt                                                            >  1:687619/1‑80 (MQ=255)
                   cgcgAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAggg                                                      <  1:687573/81‑1 (MQ=255)
                     cgAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTc                                                    <  2:687574/81‑1 (MQ=255)
                       aaCGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCAt                                                  <  2:687575/81‑1 (MQ=255)
                       aaCGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCAt                                                  <  1:687576/81‑1 (MQ=255)
                        aCGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATc                                                 <  1:687577/81‑1 (MQ=255)
                         cgcgcgAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATcc                                                >  1:687646/1‑81 (MQ=255)
                                gaacagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGc                                          >  2:687635/1‑80 (MQ=255)
                                gaacagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGc                                          >  1:687634/1‑80 (MQ=255)
                                 aacagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTa                                        <  1:687578/81‑1 (MQ=255)
                                  acagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTAc                                       >  1:687660/1‑81 (MQ=255)
                                     gaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCgg                                    <  2:687579/81‑1 (MQ=255)
                                     gaacGGACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCg                                     >  1:687655/1‑80 (MQ=255)
                                         ggACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGc                                 >  1:687676/1‑80 (MQ=255)
                                         ggACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGc                                 >  1:687643/1‑80 (MQ=255)
                                         ggACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGc                                 >  1:687633/1‑80 (MQ=255)
                                         ggACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCg                                <  1:687580/81‑1 (MQ=255)
                                         ggACTGCACCACCACCGAACCGGAGATCACACAACCACCGAAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCg                                <  1:687581/81‑1 (MQ=255)
                                          gACTGCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGa                               >  1:687644/1‑81 (MQ=255)
                                           aCTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGAt                              <  1:687582/81‑1 (MQ=255)
                                              gCACCACCACCGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCt                           >  1:687628/1‑81 (MQ=255)
                                                accaccaccGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGc                         <  2:687583/81‑1 (MQ=255)
                                                   accaccGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCAc                       >  2:687647/1‑80 (MQ=255)
                                                    ccaccGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGa                     <  1:687584/81‑1 (MQ=255)
                                                     caccGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGaa                    >  2:687636/1‑81 (MQ=255)
                                                      accGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAAt                   >  2:687650/1‑81 (MQ=255)
                                                        cGAACCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAAttt                 >  2:687683/1‑81 (MQ=255)
                                                           aCCGGAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGc              >  2:687661/1‑81 (MQ=255)
                                                              ggAGATCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTgg           >  1:687637/1‑81 (MQ=255)
                                                               gagaTCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGc          <  1:687585/81‑1 (MQ=255)
                                                                agaTCACACAACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCg         <  1:687586/81‑1 (MQ=255)
                                                                     acacaACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAAt    <  2:687588/81‑1 (MQ=255)
                                                                     acacaACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAAt    <  2:687587/81‑1 (MQ=255)
                                                                      cacaACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATg   <  1:687589/81‑1 (MQ=255)
                                                                       acaACCACCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGa  >  1:687645/1‑81 (MQ=255)
                                                                             |                                                                          
ATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGA  >  NC_000913/3568244‑3568395

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