Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,796 G→T E367* (GAA→TAA)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,7960GT95.0% 59.5 / NA 20E367* (GAA→TAA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (11/8);  total (12/8)
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

ACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTA  >  NC_000913/1756743‑1756858
                                                     |                                                              
aCCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCt                                                  <  1:290849/68‑1 (MQ=255)
 ccGTAGACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTc                                                 >  1:28824/1‑68 (MQ=255)
   gTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGACGCTCCg                                               <  1:689991/68‑1 (MQ=255)
       aCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGa                                           >  2:649320/1‑68 (MQ=255)
       aCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGa                                           >  2:200091/1‑68 (MQ=255)
       aCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGa                                           <  2:101682/68‑1 (MQ=255)
               ttACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGtt                                   <  2:2263/68‑1 (MQ=255)
               ttACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGtt                                   <  2:197154/68‑1 (MQ=255)
               ttACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGtt                                   <  1:542110/68‑1 (MQ=255)
                 aCCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAGGCTGGATGCTCCGCTGATCGTGGTTGc                                 >  1:179884/1‑68 (MQ=255)
                      aGCGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTc                            >  2:572070/1‑68 (MQ=255)
                        cGGTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAg                          >  1:17825/1‑68 (MQ=255)
                          gTATGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAggg                        <  1:546887/68‑1 (MQ=255)
                            aTGCCGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCg                      <  2:579018/68‑1 (MQ=255)
                                cGTGGTGCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTaa                  >  1:219934/1‑68 (MQ=255)
                                   ggtgCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGGTCGTGGTTGCTACTCAGGGCGGTAAATc               >  1:21333/1‑68 (MQ=255)
                                      gCCGTTGAAACTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGc            >  2:738398/1‑68 (MQ=255)
                                              aaCTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAg    >  2:378767/1‑68 (MQ=255)
                                               aCTGCTTAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGt   >  1:272378/1‑68 (MQ=255)
                                               aCTGCTTAAAAACTGGATGCTCCGCTGATCGTGGATGCTACTCAGGGCGGTAAATCTGCTAGCGCAGt   >  2:87294/1‑68 (MQ=255)
                                                ctgcttAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTa  <  2:79143/68‑1 (MQ=255)
                                                     |                                                              
ACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTA  >  NC_000913/1756743‑1756858

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