Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,108 A→T 9.9% N51I (AAT→ATT)  malY → negative regulator of MalT activity/cystathionine beta‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,1080AT 9.9% 101.0 / 2.7 40N51I (AAT→ATT) malYnegative regulator of MalT activity/cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base T (2/2);  total (17/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 = 8.19e-01

TTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCAT  >  NC_000913/1700967‑1701246
                                                                                                                                             |                                                                                                                                          
tgtCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATGATGGCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAAt                                                                                                                                           <  1:780591/141‑1 (MQ=255)
      aGGTCGTAGATCGTCAGGGCATATGGTGTACACAGTGGGAGTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTGACCTTTCACGATTTGAGACATGGATTTTGCCACTGCCCCCTGCAGTATCGAGGCGCTGAATCAgcgc                                                                                                                                     <  2:686149/143‑1 (MQ=255)
        gtcgtcGTTCGTCACGTCCCATGGCGCACACAGTGGGAGTATGTCGCTGACCGTTTCGGCGCTGCTGACCTGTTACCCTTCACGATTTCAGACAGGGATTTTGGCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCt                                                                                                                                   <  2:305523/143‑1 (MQ=255)
        gTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCt                                                                                                                                   >  1:474313/1‑143 (MQ=255)
         tCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGCGGCGCTGAATCAGCGCCTg                                                                                                                                  >  1:467246/1‑143 (MQ=255)
           tgcgATCGTCACGGCACATGCTGTACACCGTGGGACTATTTCTCTGACCGTGTCGGCACCGCTGACCTGTTACCGTTCACGAGTTCATACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcga                                                                                                                                <  2:114741/140‑4 (MQ=255)
                                     acaGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCCCGGCGTATTTGGCTACAGCCGCt                                                                                                      >  2:650723/1‑143 (MQ=255)
                                         tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGag                                                                                                  >  2:272225/1‑142 (MQ=255)
                                             gtcaggccgcgGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                              <  2:791547/136‑1 (MQ=255)
                                                aTGTCGCTGGCCGCGTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                           <  2:362316/143‑1 (MQ=255)
                                                  gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGa                                                                                         <  1:762663/143‑1 (MQ=255)
                                                  gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATGATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGa                                                                                         <  2:377883/143‑1 (MQ=255)
                                                  gTCGATGACCGTTTCGGCACTGCTGACCTGTTTCCGTTCATGATTTCAGACATGGATTTTGCCATTGACCCCTGCATTATCGAGGCGCGGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGa                                                                                         <  1:824155/143‑1 (MQ=255)
                                                      cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGttt                                                                                     >  2:717833/1‑143 (MQ=255)
                                                            gTTTCGGCAGTGTTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATGATCGGGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcgcg                                                                               <  2:88021/143‑1 (MQ=255)
                                                                     gtggttaCCTGTTCCCGCTCACGATTGCAGACATGGATTTTGCCACTGCACCCTGCGATGTTGAGGCGCTGAATCAGCGCATCATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                                                                      <  2:493274/137‑1 (MQ=255)
                                                                             cTGTTCCCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCGTTATCGATGCGCTGAATCAGCCGCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtt                                                              <  2:668586/143‑1 (MQ=255)
                                                                              tGTTACCGTTCACGATTTCAGACATGGATTTTGTCACTGCCCCCTGCATTATCGAGGCGCTGAATCATCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                             <  2:307162/143‑1 (MQ=255)
                                                                              tGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCGGATGCACGGCGTATTCGGCTCCAGCCGCCGGAAAAACGATGCATTTCGCGCGCCTATTGCCCACTGGatt                                                             >  2:285731/1‑143 (MQ=255)
                                                                                         aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAgga                                                  >  1:540630/1‑141 (MQ=255)
                                                                                          cGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGCTGAGATTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAt                                                 >  1:558712/1‑143 (MQ=255)
                                                                                          cGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAt                                                 <  1:717833/143‑1 (MQ=255)
                                                                                                cAGACATGGAGTGTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc                                           <  1:293094/143‑1 (MQ=255)
                                                                                                      tgtaTTTTGTCACTGCCCCCTGCATTATCTTGGCGCTTAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGATAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCTCCCAGCATTACACCGCCATc                                     <  2:68636/140‑1 (MQ=255)
                                                                                                      tGGATTTTGCCACTGCCCCCTGCAGTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATc                                     <  1:60334/143‑1 (MQ=255)
                                                                                                          ttttGCCACTGCCCCCTGCATTAGCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGAtt                                 <  2:152404/143‑1 (MQ=255)
                                                                                                             tGCCTCTTCCCCTTGCATTCTCGCGCCGCTGTATCAGCGCCTGGTGCACGGCGTATTTGGCTACAGCCGCTGGTAAATCTATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATtctc                              <  2:391810/143‑1 (MQ=255)
                                                                                                              gtcACTGCCCCAGTGATTAGCGAGGCGCGGAATCAGCGCCTGATGCACGGCGTATTCGGCTACAGCCGCTGGAAAAACCATGAGTCTCGCGCGGCTATGGCCCACTGTTTTTCCACCCAGCATTACACCGTCATCGATTCTCa                             <  2:601346/141‑1 (MQ=255)
                                                                                                              gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCa                             <  2:243153/143‑1 (MQ=255)
                                                                                                               ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGCTTCTCAg                            >  1:152500/1‑143 (MQ=255)
                                                                                                                  tggCCCCCTGCTTTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACg                         <  2:76359/141‑1 (MQ=255)
                                                                                                                     cccccTGCATTATCGAGGCCCTAAATAAGCACGTTATTCTCGGGGTATCTGGCTACAGACGCTGGAAACACGATGAGTTTGTCGCGGCTATTACCCGCTGGGTTTCCACTCACCCTTACGCCGACACCGAGTCTCAGACggcg                      >  2:100786/1‑141 (MQ=255)
                                                                                                                          tGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACGGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTGTCCACCCAGCATTACACCGCGATCGAGTGTGAGACGggggtggg                 >  1:19904/1‑141 (MQ=255)
                                                                                                                                 tCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCTACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCt          >  1:670654/1‑143 (MQ=255)
                                                                                                                                   gAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCACACGCCGGTGTATGGCCCTTc        >  2:119065/1‑143 (MQ=255)
                                                                                                                                     ggCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGAGGGGGGTGTATGGCCCTTCTg      >  1:270891/1‑143 (MQ=255)
                                                                                                                                     ggCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGCTCTATGGCCCTCCTg      >  2:472469/1‑143 (MQ=255)
                                                                                                                                      gcgcTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTGTCTCGCGGCTATGGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGGTGCTCAGACGGGGGTGTATGGCCCttgcgc     >  2:259202/1‑139 (MQ=255)
                                                                                                                                         cTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATGGCCCACTGGCTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATAGCCCTTCTGGCAt  >  1:482796/1‑143 (MQ=255)
                                                                                                                                         cTGAACCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGCCTATTGACCTCTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCAt  <  2:621403/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
TTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCAT  >  NC_000913/1700967‑1701246

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