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)
 gtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGagc                                                                                         <  1:426008/100‑1 (MQ=255)
       gTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGc                                                                                   <  2:1003496/100‑1 (MQ=255)
        tGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCa                                                                                  >  2:237222/1‑100 (MQ=255)
         gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCat                                                                                 >  1:60004/1‑100 (MQ=255)
         gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCat                                                                                 >  1:959552/1‑100 (MQ=255)
          cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGGTT‑GGGCGAGCCGACGCat                                                                                 >  2:1045029/1‑99 (MQ=255)
               aGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTccc                                                                           <  2:609305/100‑1 (MQ=255)
                gCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCACATATTCCCt                                                                          >  1:615175/1‑100 (MQ=255)
                  tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTgg                                                                        <  1:1273197/100‑1 (MQ=255)
                  tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTgg                                                                        <  2:520487/100‑1 (MQ=255)
                    gagaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAg                                                                      <  1:237224/100‑1 (MQ=255)
                     agaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGt                                                                     >  2:1041722/1‑100 (MQ=255)
                      gaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAg                                                                      <  2:25062/98‑1 (MQ=255)
                       ntCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTaa                                                                   >  1:1007151/2‑100 (MQ=255)
                         cATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAAtg                                                                 <  2:60003/100‑1 (MQ=255)
                           tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAAtgtg                                                               <  2:216012/100‑1 (MQ=255)
                           tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAAtgtg                                                               <  2:1272471/100‑1 (MQ=255)
                            ggCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                              >  1:511958/1‑100 (MQ=255)
                                     gtgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgag                                                     <  2:185127/100‑1 (MQ=255)
                                     gtgCGAGCGTCATGGTGTGCGGGTTATTGCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgag                                                     <  2:1086511/100‑1 (MQ=255)
                                        ccaGCGTCAAGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGAc                                                   >  2:1126066/3‑99 (MQ=255)
                                           gCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATCTGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGCGAGTgg                                                >  1:270112/1‑99 (MQ=255)
                                                 tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTATTGTGGCTCGCGGAGACTGGGCGTTGc                                         <  1:794500/100‑1 (MQ=255)
                                                 tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                         >  2:319175/1‑100 (MQ=255)
                                                   gtgtGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCtn                                       <  2:1007142/100‑2 (MQ=255)
                                                    tgtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTaa                                      >  1:98505/1‑100 (MQ=255)
                                                     gtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAAc                                     <  1:755762/100‑1 (MQ=255)
                                                       gCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:391331/1‑100 (MQ=255)
                                                       gCGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:693462/1‑100 (MQ=255)
                                                        cGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                   >  1:186426/1‑100 (MQ=255)
                                                        cGGGTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTc                                  >  1:706871/1‑100 (MQ=255)
                                                          ggTTATTTCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCTTTGCTAACGTCgg                                <  1:876580/100‑1 (MQ=255)
                                                          ggTTATTTCCGATGAAATCCATATGGAGATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCgg                                <  2:108773/100‑1 (MQ=255)
                                                               tttCCGATGAAATCAATATGGATATGGTTT‑GGGCGAGGAGCCGCATATTACCTGGAGTAATGTGGCGCGAGGAGACTGGGCGTTGCTAACGTCGGGCttg                           >  2:241661/1‑98 (MQ=255)
                                                                ttCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCg                           <  2:628831/99‑1 (MQ=255)
                                                                ttCCGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                          <  1:960343/100‑1 (MQ=255)
                                                                tgccGATGAAATCCTTATGGATAT‑GTTTTGGTCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                          >  2:386028/3‑100 (MQ=255)
                                                                  ccGATGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                         >  2:19230/1‑99 (MQ=255)
                                                                     aTGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGt                     >  2:657996/1‑100 (MQ=255)
                                                                      tGAAATCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGtt                    <  1:672323/100‑1 (MQ=255)
                                                                        aaaTCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGATTc                  <  2:197030/100‑1 (MQ=255)
                                                                          aTCCATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                >  2:741081/1‑100 (MQ=255)
                                                                             cATATGGATATGGTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAatat             >  1:65380/1‑100 (MQ=255)
                                                                                atggatatggTTT‑GGGCGAGCAGCCGCATATTCCCTGTAGTTATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTc           <  2:916167/99‑1 (MQ=255)
                                                                                   gatatggTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGc       >  2:505026/1‑100 (MQ=255)
                                                                                        ggTTT‑GGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGa  >  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: ‑