Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,701,585 Δ1 bp coding (629/1173 nt) malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,5820G.97.8% 195.6 / ‑2.4 46coding (626/1173 nt)malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (22/22);  minor base T (1/0);  total (24/22)
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

TGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGA  >  NC_000913/1701489‑1701677
                                                                                             |                                                                                               
tgtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcgag                                                                                           >  2:511114/1‑96 (MQ=255)
tgtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcgag                                                                                           <  1:906078/100‑5 (MQ=255)
 gtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGagc                                                                                         <  1:426008/100‑1 (MQ=255)
       gTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGc                                                                                   <  2:1003496/100‑1 (MQ=255)
        tGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCa                                                                                  >  2:237222/1‑100 (MQ=255)
         gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCat                                                                                 >  1:60004/1‑100 (MQ=255)
         gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCat                                                                                 >  1:959552/1‑100 (MQ=255)
          cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGGTTGGGCGAGCCGACGCat                                                                                 >  2:1045029/1‑99 (MQ=255)
               aGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTccc                                                                           <  2:609305/100‑1 (MQ=255)
                gCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCACATATTCCCt                                                                          >  1:615175/1‑100 (MQ=255)
                  tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTgg                                                                        <  1:1273197/100‑1 (MQ=255)
                  tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTgg                                                                        <  2:520487/100‑1 (MQ=255)
                    gagaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAg                                                                      <  1:237224/100‑1 (MQ=255)
                     agaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGt                                                                     >  2:1041722/1‑100 (MQ=255)
                      gaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAg                                                                      <  2:25062/98‑1 (MQ=255)
                       ntCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTaa                                                                   >  1:1007151/2‑100 (MQ=255)
                         cATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAAtg                                                                 <  2:60003/100‑1 (MQ=255)
                           tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAAtgtg                                                               <  2:216012/100‑1 (MQ=255)
                           tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAAtgtg                                                               <  2:1272471/100‑1 (MQ=255)
                            ggCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                              >  1:511958/1‑100 (MQ=255)
                                     gtgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgag                                                     <  2:185127/100‑1 (MQ=255)
                                     gtgCGAGCGTCATGGTGTGCGGGTTATTGCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgag                                                     <  2:1086511/100‑1 (MQ=255)
                                        ccaGCGTCAAGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGAc                                                   >  2:1126066/3‑99 (MQ=255)
                                           gCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATCTGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGCGAGTgg                                                >  1:270112/1‑99 (MQ=255)
                                                 tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTATTGTGGCTCGCGGAGACTGGGCGTTGc                                         <  1:794500/100‑1 (MQ=255)
                                                 tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                         >  2:319175/1‑100 (MQ=255)
                                                   gtgtGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCtn                                       <  2:1007142/100‑2 (MQ=255)
                                                    tgtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTaa                                      >  1:98505/1‑100 (MQ=255)
                                                     gtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAAc                                     <  1:755762/100‑1 (MQ=255)
                                                       gCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:391331/1‑100 (MQ=255)
                                                       gCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:693462/1‑100 (MQ=255)
                                                        cGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:186426/1‑100 (MQ=255)
                                                        cGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTc                                  >  1:706871/1‑100 (MQ=255)
                                                          ggTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCTTTGCTAACGTCgg                                <  1:876580/100‑1 (MQ=255)
                                                          ggTTATTTCCGATGAAATCCATATGGAGATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCgg                                <  2:108773/100‑1 (MQ=255)
                                                               tttCCGATGAAATCAATATGGATATGGTTTGGGCGAGGAGCCGCATATTACCTGGAGTAATGTGGCGCGAGGAGACTGGGCGTTGCTAACGTCGGGCttg                           >  2:241661/1‑98 (MQ=255)
                                                                ttCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCg                           <  2:628831/99‑1 (MQ=255)
                                                                ttCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                          <  1:960343/100‑1 (MQ=255)
                                                                tgccGATGAAATCCTTATGGATATGTTTTGGTCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                          >  2:386028/3‑100 (MQ=255)
                                                                  ccGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                         >  2:19230/1‑99 (MQ=255)
                                                                     aTGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGt                     >  2:657996/1‑100 (MQ=255)
                                                                      tGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGtt                    <  1:672323/100‑1 (MQ=255)
                                                                        aaaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGATTc                  <  2:197030/100‑1 (MQ=255)
                                                                          aTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                >  2:741081/1‑100 (MQ=255)
                                                                             cATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAatat             >  1:65380/1‑100 (MQ=255)
                                                                                atggatatggTTTGGGCGAGCAGCCGCATATTCCCTGTAGTTATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTc           <  2:916167/99‑1 (MQ=255)
                                                                                   gatatggTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGc       >  2:505026/1‑100 (MQ=255)
                                                                                        ggTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGa  >  1:158560/1‑100 (MQ=255)
                                                                                             |                                                                                               
TGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGA  >  NC_000913/1701489‑1701677

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