Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,932,990 A→G E362G (GAA→GGA)  pyk → pyruvate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,932,9900AG100.0% 116.6 / NA 31E362G (GAA→GGA) pykpyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/19);  total (12/19)

GCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGG  >  NZ_CP009273/1932902‑1933071
                                                                                        |                                                                                 
gCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGggag                                                                                 <  2:236658/91‑1 (MQ=255)
    tGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGc                                                                             <  1:457944/91‑1 (MQ=255)
         cgcgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTg                                                                        <  2:51778/91‑1 (MQ=255)
           cgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGcc                                                                      <  2:139673/91‑1 (MQ=255)
                gCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTc                                                                 >  1:211821/1‑91 (MQ=255)
                gCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTc                                                                 <  1:295564/91‑1 (MQ=255)
                  cTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAg                                                               <  2:94571/91‑1 (MQ=255)
                     ggTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCaa                                                            >  2:390050/1‑91 (MQ=255)
                         cGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTa                                                        >  2:15709/1‑91 (MQ=255)
                          ggAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTAc                                                       <  1:99155/91‑1 (MQ=255)
                            aaaaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGc                                                     <  1:390050/91‑1 (MQ=255)
                                 aTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTa                                                <  2:468734/91‑1 (MQ=255)
                                            cAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGaaa                                     <  2:792644/91‑1 (MQ=255)
                                                      aaaCACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACgg                           <  1:521691/91‑1 (MQ=255)
                                                       aaCACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGc                          <  1:687879/91‑1 (MQ=255)
                                                        acacCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCg                         <  2:685631/91‑1 (MQ=255)
                                                           ccGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatc                      >  1:429465/1‑91 (MQ=255)
                                                            cGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatca                     <  2:364349/91‑1 (MQ=255)
                                                            cGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatca                     <  1:355462/91‑1 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACgg                           >  1:730419/1‑84 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACgg                           <  2:730419/84‑1 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcat                    >  2:162926/1‑91 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcat                    >  1:706752/1‑91 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcat                    >  2:417064/1‑91 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcat                    >  1:600691/1‑91 (MQ=255)
                                                             gTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAAAGTAAGCAGGTAAACACCACAAAGGCGGTAAGGGGatcat                    >  2:637044/1‑91 (MQ=255)
                                                              tCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcatc                   <  2:645961/91‑1 (MQ=255)
                                                                  gACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCa               <  1:519062/91‑1 (MQ=255)
                                                                      ttCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGAc           <  1:15709/91‑1 (MQ=255)
                                                                      ttCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGAc           <  1:560258/91‑1 (MQ=255)
                                                                           ttCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGTAATGTGTTACGGCGGTCCCCACCATGACCGAAt      <  2:811584/91‑1 (MQ=255)
                                                                               aCAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCggg  >  1:464394/1‑91 (MQ=255)
                                                                               aCAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCggg  >  1:135064/1‑91 (MQ=255)
                                                                                        |                                                                                 
GCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGG  >  NZ_CP009273/1932902‑1933071

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