Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,938,397 G→A R250H (CGT→CAT)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,938,3970GA97.8% 170.9 / ‑4.4 46R250H (CGT→CAT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGT  >  NC_000913/1938249‑1938542
                                                                                                                                                    |                                                                                                                                                 
ggCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATg                                                                                                                                                 <  1:508720/151‑1 (MQ=255)
    aaGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGa                                                                                                                                             <  1:405529/151‑1 (MQ=255)
           gAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCgg                                                                                                                                       <  2:226306/150‑1 (MQ=255)
           gAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGgt                                                                                                                                      <  1:344641/151‑1 (MQ=255)
            aaCTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGgtg                                                                                                                                     >  2:626762/1‑151 (MQ=255)
                 tGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGa                                                                                                                                 >  1:659344/1‑150 (MQ=255)
                 tGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGa                                                                                                                                 >  1:630286/1‑150 (MQ=255)
                    ccGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAAtt                                                                                                                             <  1:139829/151‑1 (MQ=255)
                    ccGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAAtt                                                                                                                             <  2:690191/151‑1 (MQ=255)
                     cGTCTCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTg                                                                                                                            <  2:164815/151‑1 (MQ=255)
                             ggCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCg                                                                                                                    >  2:642949/1‑151 (MQ=255)
                                      tGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTc                                                                                                           >  2:295801/1‑151 (MQ=255)
                                         aGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGc                                                                                                        >  2:151006/1‑151 (MQ=255)
                                         aGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGc                                                                                                        >  1:647026/1‑151 (MQ=255)
                                            aTGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAtt                                                                                                     <  1:145855/151‑1 (MQ=255)
                                             tgtgATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTc                                                                                                    >  2:328562/1‑151 (MQ=255)
                                             tgtgATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTc                                                                                                    >  2:36923/1‑151 (MQ=255)
                                                   gCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                                                                >  1:300742/1‑149 (MQ=255)
                                                       aaaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGtt                                                                                          <  2:667296/151‑1 (MQ=255)
                                                            gttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATcc                                                                                     >  1:308873/1‑151 (MQ=255)
                                                              tgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgt                                                                                   >  2:267560/1‑151 (MQ=255)
                                                              tgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                                                                    <  1:197797/150‑1 (MQ=255)
                                                                ccAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcg                                                                                 <  2:758930/151‑1 (MQ=255)
                                                                         aaCGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCa                                                                        <  2:537675/151‑1 (MQ=255)
                                                                          aCGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAg                                                                       <  1:179895/151‑1 (MQ=255)
                                                                                ggAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAAc                                                                 >  1:181696/1‑151 (MQ=255)
                                                                                   aGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGa                                                              <  2:263417/151‑1 (MQ=255)
                                                                                          tGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTa                                                        <  2:75312/150‑1 (MQ=255)
                                                                                           gCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAAt                                                      >  1:719785/1‑151 (MQ=255)
                                                                                           gCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAAt                                                      <  1:204814/151‑1 (MQ=255)
                                                                                             aGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCa                                                    <  1:275895/151‑1 (MQ=255)
                                                                                              gCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCAc                                                   <  1:601666/151‑1 (MQ=255)
                                                                                                 aGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGc                                                >  1:696705/1‑151 (MQ=255)
                                                                                                 aGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGc                                                >  1:32909/1‑151 (MQ=255)
                                                                                                             aTGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAgatgat                                    >  2:266877/1‑151 (MQ=255)
                                                                                                               gACATCTTCATCGCCTCTGACGTGGTAATGGTTGCTCATGGCTACCTCGGTGTGGAAATTGGCGATCCGGAACTGGTCGGCATTCAGAAAGCGTTTATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATgg                                  <  1:719640/151‑1 (MQ=255)
                                                                                                                  atcatcCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGt                               <  2:412133/151‑1 (MQ=255)
                                                                                                                    catcCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTc                              <  2:144193/150‑1 (MQ=255)
                                                                                                                           gCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGAtt                       >  2:492647/1‑150 (MQ=255)
                                                                                                                                tGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAAc                 >  1:382391/1‑151 (MQ=255)
                                                                                                                                tGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAAc                 >  2:711851/1‑151 (MQ=255)
                                                                                                                                 gACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAAcc                <  1:719601/151‑1 (MQ=255)
                                                                                                                                      ggTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATg           >  2:423569/1‑151 (MQ=255)
                                                                                                                                            ggTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcg     >  2:749376/1‑151 (MQ=255)
                                                                                                                                               tGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcgcg   <  1:80262/150‑1 (MQ=255)
                                                                                                                                                gCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGt  <  1:470446/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGT  >  NC_000913/1938249‑1938542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑