Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,026,927 A→T V30E (GTG→GAG)  mgsA ← methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,9270AT96.3% 74.1 / NA 28V30E (GTG→GAG) mgsAmethylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (11/15);  minor base G (0/1);  total (11/17)
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

ATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAGTTAACTTTCCTACAGATTACTGTAAGCACTTA  >  NC_000913/1026787‑1027062
                                                                                                                                            |                                                                                                                                       
aTGCGCCAACCTCCTGGGCACCCCCCATTGCGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTtcctcc                                                                                                                                       <  1:448015/143‑1 (MQ=255)
      gagccTGCGGGCCACCCCCCTTTGGGGCACTCTGCATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTACCGCTTGATGCCGTTCCTCCCAGCtc                                                                                                                                 <  2:9223/140‑1 (MQ=255)
            gctgGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCCGc                                                                                                                           >  1:527515/1‑143 (MQ=255)
                  cACCCCCCAGTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTTTGCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCAtttgt                                                                                                                     <  1:76986/143‑1 (MQ=255)
                     ccccccTTGGGCCCCTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCAtttgtttg                                                                                                                  <  2:400713/143‑1 (MQ=255)
                      ccccATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGGCGTTCCTCCCCGCTCATCAGCATTTGTTTGc                                                                                                                 >  1:70277/1‑143 (MQ=255)
                       cccATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCa                                                                                                                <  2:292268/143‑1 (MQ=255)
                             gggCCGCTCAACATCGCGTGGACGTTCAGGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATc                                                                                                          <  1:127779/143‑1 (MQ=255)
                                 cACTCACCATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCgtgt                                                                                                      <  1:480372/143‑1 (MQ=255)
                                  aCTCACCAGCGCGTTGACCTGCATGCCGGTCGCGCGGGAAAGTAATTTACCGGTAGTGCCTGTTGCATACTGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGtgtg                                                                                                     <  2:6851/143‑1 (MQ=255)
                                        cATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGccacca                                                                                               >  2:474564/1‑143 (MQ=255)
                                          tCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAgc                                                                                             >  2:527576/1‑143 (MQ=255)
                                                     ttCATGCCGGGCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGttt                                                                                  <  1:415758/143‑1 (MQ=255)
                                                         tGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCcgc                                                                              >  2:152647/1‑143 (MQ=255)
                                                            cGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCa                                                                           <  1:452710/143‑1 (MQ=255)
                                                                        aaaTTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTACCGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGa                                                               <  2:214525/143‑1 (MQ=255)
                                                                                  aCCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGtt                                                     <  1:373128/143‑1 (MQ=255)
                                                                                    cGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCCCGCAGGTAACGTGCGAGTCGTCAGTTcc                                                   >  2:586363/1‑143 (MQ=255)
                                                                                          tGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATg                                             <  1:480103/143‑1 (MQ=255)
                                                                                          tGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATATGTTTGCAGTGATCGTGTGCCACCAGAGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATaacg                                             >  2:208717/1‑141 (MQ=255)
                                                                                           gCCGGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGt                                            <  1:610646/143‑1 (MQ=255)
                                                                                             cTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCGGGTAAAGTGCGAGTCGTCAGTTCCATAATGtcc                                          >  1:266357/1‑141 (MQ=255)
                                                                                                  gCATACAGTCCGTGTTGTCCCAGGAACGGTTGATGCCGTCCCTCCCAGCTCATCAGCATGTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCg                                     <  1:576808/143‑1 (MQ=255)
                                                                                                     tACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAg                                  >  1:223745/1‑143 (MQ=255)
                                                                                                                    tCCAGTAACGGTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGGCACTTCCATAATGTACACCCGTAGTTAACTTCCCTCCaa                   >  2:81740/1‑142 (MQ=255)
                                                                                                                     ccAGCAACGGTTGCTGCCGTTCCGCCCAGCCCTTCAGCCTTTGTTTGCAGTGACCGCGTGCCACCGCCGCCATATTTTGCCGCGCAGGTAAAGTGCGAGTCGGCAGTTCCATAATGTACATCCGTAGTTAACTTTCCTACAGa                  <  1:203022/143‑1 (MQ=255)
                                                                                                                              gTTGATGCCGTTCCTCCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAGTTAACTTTCCTACAGATTACTGTaa         >  1:389860/1‑143 (MQ=255)
                                                                                                                                tGATGCCGTTCCTCCCAGCTCATCAGCATGTGTTTGCAGTGATCGTGTTCCACCAGCGCAATATGTTTCCGCTCATGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAGTTAACTTTCCTACAGATTACTGTAAGc       <  2:73564/143‑1 (MQ=255)
                                                                                                                                     ccGTCCCCCCAAGCCCACCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAGTTCACTTTCCTACAGATTACTGTAAGCACTTa  <  1:440854/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
ATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCTCATCAGCATTTGTTTGCAGTGATCGTGTGCCACCAGCGCAATATGTTTCCGCGCAGGTAAAGTGCGAGTCGTCAGTTCCATAATGTACATCCGTAGTTAACTTTCCTACAGATTACTGTAAGCACTTA  >  NC_000913/1026787‑1027062

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