Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,497,023 A→T 9.4% intergenic (+31/+34) ypdK → / ← alaC inner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,497,0230AT 9.4% 127.7 / 4.7 53intergenic (+31/+34)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base T (2/3);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAG  >  NC_000913/2496878‑2497170
                                                                                                                                                 |                                                                                                                                                   
gAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAg                                                                                                                                                 >  2:766096/1‑150 (MQ=255)
  aGTCATTCAGGCGTATCAGTAATTACGCGGAGAGATGTATAGTGATATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                               <  2:152465/150‑1 (MQ=255)
     cATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTcc                                                                                                                                            >  2:1267500/1‑150 (MQ=255)
     cATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTcc                                                                                                                                            >  2:246158/1‑150 (MQ=255)
         cAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGtt                                                                                                                                        <  2:799400/150‑1 (MQ=255)
           ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGTGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttt                                                                                                                                      <  1:1186796/150‑2 (MQ=255)
           ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAGACAGGAGCCATCAGGCTCCTGTTTg                                                                                                                                      <  2:720555/150‑1 (MQ=255)
           ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttt                                                                                                                                      <  1:175906/150‑2 (MQ=255)
             cggaaCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGGCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCGCCTGTTTGCa                                                                                                                                    <  1:1076630/147‑1 (MQ=255)
               tAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCAtt                                                                                                                                  <  2:866168/150‑1 (MQ=255)
                aaCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCATTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttttcatt                                                                                                                                 <  1:1059572/150‑7 (MQ=255)
                    ctAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTAATGGTCATCGTTTGTGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTga                                                                                                                             <  1:1306351/149‑1 (MQ=255)
                                 gagaTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGACATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTAt                                                                                                                >  2:1236902/1‑150 (MQ=255)
                                     tGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcg                                                                                                            >  1:921472/1‑150 (MQ=255)
                                     tGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcg                                                                                                            >  1:544784/1‑150 (MQ=255)
                                       tAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcgcg                                                                                                          >  2:957854/1‑150 (MQ=255)
                                       tAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCACTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcgcg                                                                                                          >  1:732152/1‑150 (MQ=255)
                                             gAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCg                                                                                                    >  1:1267461/1‑150 (MQ=255)
                                              aaaTATTTCTTTATGGGCATTTGTTTTATGGACATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGAAATGAAAAACTGGAGCCATCAGGCTCCTGTTTACATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                   <  2:634076/150‑1 (MQ=255)
                                                  atTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAat                                                                                               >  2:982992/1‑150 (MQ=255)
                                                   tttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAata                                                                                              >  2:317116/1‑150 (MQ=255)
                                                     tCTGTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCACTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAatatt                                                                                            >  1:20686/1‑149 (MQ=255)
                                                       tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGaa                                                                                                >  1:511187/1‑144 (MQ=255)
                                                       tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGaa                                                                                                <  2:511187/144‑1 (MQ=255)
                                                          ttGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgc                                                                                       <  2:692597/149‑1 (MQ=255)
                                                                    gttttATGGTCATCTTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCggg                                                                             <  2:927796/149‑1 (MQ=255)
                                                                            gTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAAc                                                                     <  1:1042381/150‑1 (MQ=255)
                                                                              cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCg                                                                   >  2:424055/1‑150 (MQ=255)
                                                                                  gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCgg                                                               >  1:1299294/1‑150 (MQ=255)
                                                                                  gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTATCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCgg                                                               >  1:327024/1‑150 (MQ=255)
                                                                                       cgCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACAGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGGTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACgg                                                          >  1:652352/2‑150 (MQ=255)
                                                                                         ccGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTACACAGTCGGCACGGaa                                                        >  1:872737/1‑150 (MQ=255)
                                                                                                      ccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTAGGCACGGAACATCGCTTTAATc                                           >  1:246166/1‑150 (MQ=255)
                                                                                                       cTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                          <  1:50446/150‑1 (MQ=255)
                                                                                                       cTGATGATCTAAAAGCAAACCGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                          <  1:632318/150‑1 (MQ=255)
                                                                                                              tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTACTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGAt                                   >  2:298070/1‑150 (MQ=255)
                                                                                                               cTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                  <  1:694609/150‑1 (MQ=255)
                                                                                                                 aaaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGc                                >  1:437635/1‑150 (MQ=255)
                                                                                                                  aaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGGAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGcc                               <  1:618792/150‑1 (MQ=255)
                                                                                                                   aaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCt                              <  2:1267461/150‑1 (MQ=255)
                                                                                                                       aaaCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                          >  1:640826/1‑150 (MQ=255)
                                                                                                                               cAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                  <  1:1206766/150‑1 (MQ=255)
                                                                                                                               cAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                  <  1:1267500/150‑1 (MQ=255)
                                                                                                                               aaaTGAAAAACAGGAGACATCATGCTCCTGTTTGCATTTCTGACGAAGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                  <  1:246158/149‑1 (MQ=255)
                                                                                                                                    aaaaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCgg             <  2:544784/150‑1 (MQ=255)
                                                                                                                                    aaaaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGAGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCGCTCTGAGGGCCTGACGAATACGGTCGCgg             >  2:5423/1‑150 (MQ=255)
                                                                                                                                     caaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTATAATCCCTCTCATGGCCTGACGAATACGGTCGCGGt            >  2:344100/2‑150 (MQ=255)
                                                                                                                                      aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAAGATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtt           <  1:138648/150‑1 (MQ=255)
                                                                                                                                      aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGAGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTGTAATGCCTCTGATGGCCTGACGAGTTCGGTCGCGGtt           >  2:1065633/1‑150 (MQ=255)
                                                                                                                                        aCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCAGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt         >  2:754165/1‑150 (MQ=255)
                                                                                                                                        aCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt         >  1:1242319/1‑150 (MQ=255)
                                                                                                                                              gCCATCAGGCTCCTGTTTGCATTTCTGCCGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTACACCGTCGGCACGGAACATCGCTTTAAGCCCTCTGTTGGCCTGACGAATACGGTCGCGGTTTTCAATCa   >  2:176352/1‑150 (MQ=255)
                                                                                                                                               ccATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAg  >  1:446121/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAG  >  NC_000913/2496878‑2497170

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