Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,374 G→A G226D (GGC→GAC)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,3740GA97.7% 120.5 / NA 43G226D (GGC→GAC) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/22);  total (20/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 = 1.00e+00

TTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGC  >  NC_000913/1756235‑1756506
                                                                                                                                           |                                                                                                                                    
gttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCt                                                                                                                                   <  2:338535/142‑1 (MQ=255)
         gAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACtt                                                                                                                          >  2:905762/1‑143 (MQ=255)
           aCAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAACATCGAAAACCGGGAAGACCTCAACAACttgg                                                                                                                        >  1:674382/1‑141 (MQ=255)
           aCAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTcg                                                                                                                        <  2:284006/143‑1 (MQ=255)
                   tAGTCTTTGTTGCTGTTTTCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGGGAGCAGCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATc                                                                                                                <  1:750403/143‑1 (MQ=255)
                            ttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGCCGGCGAAAACATCCACATCATCTCCAAAATCGCAAACCAGGAAGACCCCAACAACTCCGACGAACTCCTCGAAGcc                                                                                                       >  2:461130/1‑143 (MQ=255)
                                 gctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGa                                                                                                  >  1:775384/1‑143 (MQ=255)
                                   ttccttcgcccGTAAGCCTTCGGCCGTTATCGAAATCCGTGTGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAACCCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACg                                                                                                <  2:459852/133‑1 (MQ=255)
                                   ttCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGACAACCAGGCAGACCTCGATAACTTCGACGAAATCCTCGAAGCCTCTCACg                                                                                                >  1:243348/1‑143 (MQ=255)
                                       catATTCGTAAGCGTTGTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcat                                                                                            >  2:197016/3‑143 (MQ=255)
                                                  gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTgcgc                                                                                 >  1:476039/1‑143 (MQ=255)
                                                              tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTgg                                                                     <  1:82813/143‑1 (MQ=255)
                                                              tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTATGACGGCATCATGGTTGCGCGTGGCGACCggg                                                                     >  2:81462/1‑143 (MQ=255)
                                                                 cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAAACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtg                                                                  >  1:657097/1‑143 (MQ=255)
                                                                      tCCGTTTGCTCCTGAAAGCGCACGGCAGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                                             <  2:439076/143‑1 (MQ=255)
                                                                           gAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATccc                                                        <  1:428884/143‑1 (MQ=255)
                                                                            cgtacgtggaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCGCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCg                                                       <  1:849674/134‑1 (MQ=255)
                                                                            aGCACCTGCAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCg                                                       <  1:332030/143‑1 (MQ=255)
                                                                             gCACCTGAAAGCGCTCGGCGGCGAAAACATCCACACCATCTCCAAACTCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCgg                                                      <  2:347297/143‑1 (MQ=255)
                                                                              cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGt                                                     <  1:531421/143‑1 (MQ=255)
                                                                               aCCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTa                                                    <  1:905762/143‑1 (MQ=255)
                                                                                 cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTaga                                                  <  1:174376/143‑1 (MQ=255)
                                                                                 cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGCAGAAATCCCGGTaga                                                  >  2:522528/1‑143 (MQ=255)
                                                                                  tGAAAGCGCACGACACCGAAAACATCCACATAATCTCCAAACTCGAAAACGAAGAAGAACTCAACAACTTCGACGAAAACCTCGAAGCCTCTGACGGCATCATGGTTGAGCGTAGCGATCTGGGTGTGGaacacacgactgac                                                 >  2:77319/1‑131 (MQ=255)
                                                                                   gAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTggacg                                                >  1:14125/1‑141 (MQ=255)
                                                                                     aaGCTCACGGCGGCGAAAACAGCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACGTCGACGAAATCCGCGAAGCCTCTGACGGCATCAGGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa                                              <  2:701012/143‑1 (MQ=255)
                                                                                     aaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa                                              <  1:465290/143‑1 (MQ=255)
                                                                                      aGCCCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGaagaag                                             <  1:577341/143‑1 (MQ=255)
                                                                                          cACGGCGGCGAAAACATCCACAGCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTCGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTAt                                         <  1:712148/143‑1 (MQ=255)
                                                                                                ggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGTAAGACCAAAACGTCTTCGATGTTATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGc                                   <  2:784678/143‑1 (MQ=255)
                                                                                                 gcgAAAACATCCACATCACCTCCATAATCGAACACCCGGAAGACCTCGACAACTTCGACGAAATCCTCGAAGCCTCGGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc                                  <  1:333267/143‑1 (MQ=255)
                                                                                                 gcgAAAACAGCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGcc                                  <  1:326693/143‑1 (MQ=255)
                                                                                                  cgAAAACATCCTCATCATCACCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATACTTGCCCGTGGCAACCTCGGTGTACAAATCCAAGTCGACGACGTTATCTCCGccc                                 >  2:177436/1‑143 (MQ=255)
                                                                                                   gAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAGTCCGGGTAGAAGAAGTTATCTTGGCGCa                                >  1:284948/1‑143 (MQ=255)
                                                                                                        cATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaaga                           <  1:146067/143‑1 (MQ=255)
                                                                                                              catcctCTCCAAAATCGAAAACCAGGAAGACCTCAACAACATCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatga                     <  1:369174/143‑1 (MQ=255)
                                                                                                                  atcTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAACCCCGGTAGCAGAAGTCATCTTCGCCCAGAACATGATGATCGa                 >  2:807467/1‑143 (MQ=255)
                                                                                                                     tcCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTTTCTTCGCCCAGAAGATGATGATAGaaaa              >  1:355700/1‑143 (MQ=255)
                                                                                                                      ccAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTGGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaa             <  1:650050/143‑1 (MQ=255)
                                                                                                                         aaaTCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGt          >  1:578668/1‑143 (MQ=255)
                                                                                                                            tCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTCGAAATCCCGGTCGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATc       >  1:36200/1‑143 (MQ=255)
                                                                                                                               aaaaCCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAATCCTCTGACGGCATGATGGTAGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGACGAAGTTATGTTCGCCGAGAAGCGGATGACCGGAAACTGTAGCCGt    >  1:142928/1‑143 (MQ=255)
                                                                                                                                 aaCCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAACATGTATCCGTGc  >  2:852/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGC  >  NC_000913/1756235‑1756506

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