Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,108 A→T 7.4% N51I (AAT→ATT)  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,1080AT 7.4% 184.0 / 4.2 54N51I (AAT→ATT) malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base T (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGT  >  NC_000913/1700960‑1701255
                                                                                                                                                    |                                                                                                                                                   
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                     >  3:151478/1‑149 (MQ=255)
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                     >  2:323119/1‑149 (MQ=255)
   gttttttCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCa                                                                                                                                                  >  1:127333/3‑149 (MQ=255)
   gATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCa                                                                                                                                                  <  2:127333/149‑1 (MQ=255)
         tCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCt                                                                                                                                            <  1:295973/149‑1 (MQ=255)
         tCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCt                                                                                                                                            <  4:119979/149‑1 (MQ=255)
            aaGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGAt                                                                                                                                         <  4:472587/149‑1 (MQ=255)
                   tGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                                                                  <  1:293229/149‑1 (MQ=255)
                   tGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                                                                  >  1:414854/1‑149 (MQ=255)
                   tGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                                                                  <  2:118886/149‑1 (MQ=255)
                    ggATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCg                                                                                                                                 >  1:88982/1‑149 (MQ=255)
                       tCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAt                                                                                                                              >  1:121805/1‑149 (MQ=255)
                       tCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAt                                                                                                                              <  3:417487/149‑1 (MQ=255)
                        cGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAtt                                                                                                                             <  2:195878/149‑1 (MQ=255)
                         gTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAttt                                                                                                                            <  4:149446/149‑1 (MQ=255)
                             tGGCACATGTTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCt                                                                                                                        >  4:260292/1‑149 (MQ=255)
                                         tACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGa                                                                                                            >  1:476982/1‑149 (MQ=255)
                                           cacaGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                                                                          >  3:14751/1‑149 (MQ=255)
                                              aGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                                       >  1:100754/1‑149 (MQ=255)
                                              aGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                                       >  2:219990/1‑149 (MQ=255)
                                                tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                                                                     <  2:414794/149‑1 (MQ=255)
                                                 gggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                                    <  2:414854/149‑1 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGt                                                                                                >  3:288557/1‑149 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGt                                                                                                <  1:219990/149‑1 (MQ=255)
                                                       aTGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGttt                                                                                              <  4:151478/149‑1 (MQ=255)
                                                            gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcgc                                                                                         <  4:25220/149‑1 (MQ=255)
                                                             cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcgcg                                                                                        <  2:100754/149‑1 (MQ=255)
                                                              tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCgg                                                                                       <  2:420827/149‑1 (MQ=255)
                                                                    tttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTg                                                                                 >  2:384417/1‑149 (MQ=255)
                                                                     ttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                                                                                <  2:121805/149‑1 (MQ=255)
                                                                         gCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCa                                                                                                                                                  >  2:405756/1‑79 (MQ=255)
                                                                         gCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCa                                                                                                                                                  <  1:405756/79‑1 (MQ=255)
                                                                                    cTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAc                                                                 >  4:175057/1‑149 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                                    <  2:365842/111‑1 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                                    >  1:365842/1‑111 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa                                                              <  1:411630/149‑1 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa                                                              <  2:88982/149‑1 (MQ=255)
                                                                                               cACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                                                                    <  4:276039/87‑1 (MQ=255)
                                                                                               cACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAg                                                                                                                    >  3:276039/1‑87 (MQ=255)
                                                                                               cACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTaca                                                      <  3:175057/149‑1 (MQ=255)
                                                                                                aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTacac                                                     >  1:284895/1‑149 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTTGCTACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                                                                          <  2:210156/94‑1 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaa                                                                                                          >  1:210156/1‑94 (MQ=255)
                                                                                                   aTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGc                                                  >  1:368034/1‑149 (MQ=255)
                                                                                                            aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATtc                                         >  1:98934/1‑149 (MQ=255)
                                                                                                             tGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTGTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATtct                                        >  3:58466/1‑149 (MQ=255)
                                                                                                               gATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCa                                      <  1:323119/149‑1 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                               <  1:384417/149‑1 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                               <  2:284895/149‑1 (MQ=255)
                                                                                                                              cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGCCGGTGGTGTATgg                       >  1:473142/1‑149 (MQ=255)
                                                                                                                                   cATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc                                                               >  1:48841/1‑104 (MQ=255)
                                                                                                                                   cATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc                                                               <  2:48841/104‑1 (MQ=255)
                                                                                                                                        tCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCATTGGTTTTCCACCCAGCATTACACCGCCATCGATTATCAGACGGTGGTGTATGGCCCTTCTGTc             <  3:407821/149‑1 (MQ=255)
                                                                                                                                        cgtaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTc             <  4:119144/146‑1 (MQ=255)
                                                                                                                                        atcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTc             >  2:426596/4‑149 (MQ=255)
                                                                                                                                          caGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCAt           >  2:363596/2‑149 (MQ=255)
                                                                                                                                                   aaTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTAAGACGGTGGTGTATGGCCCTTCTGTGATCTATATGGt  >  3:309426/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGT  >  NC_000913/1700960‑1701255

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