Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,497,025 C→A 9.1% intergenic (+33/+32) 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,0250CA 9.1% 137.7 / 6.9 55intergenic (+33/+32)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base A (2/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

GAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGC  >  NC_000913/2496878‑2497173
                                                                                                                                                   |                                                                                                                                                    
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)
                                                                                                                                                  tCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGc  <  2:797205/150‑1 (MQ=255)
                                                                                                                                                  tCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGc  <  1:688706/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGC  >  NC_000913/2496878‑2497173

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: