Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,882,139 A→G E17G (GAG→GGG)  tnaA → tryptophanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,882,1390AG95.8% 84.5 / ‑6.0 24E17G (GAG→GGG) tnaAtryptophanase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (15/8);  minor base T (1/0);  total (16/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

CAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGAT  >  NZ_CP009273/3882051‑3882224
                                                                                        |                                                                                     
caGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGc                                                                                     >  1:560719/1‑91 (MQ=255)
 aGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGcc                                                                                    <  1:169745/91‑1 (MQ=255)
 aGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGcc                                                                                    <  1:265280/91‑1 (MQ=255)
   ggATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAg                                                                                  >  1:302030/1‑91 (MQ=255)
    gATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGTGCCAGt                                                                                 >  2:725529/1‑91 (MQ=255)
        aCTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTaaaa                                                                             <  2:683742/91‑1 (MQ=255)
            tAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCAGCATTCGTGTTATTGGGCCAGTAAAACGTa                                                                         >  1:424478/1‑91 (MQ=255)
                             ggATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGt                                                        <  1:470835/91‑1 (MQ=255)
                                     tAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGc                                                <  2:584395/91‑1 (MQ=255)
                                       aTGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCaa                                              >  2:177723/1‑91 (MQ=255)
                                       aTGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCaa                                              >  2:177722/1‑91 (MQ=255)
                                                        tctcCCTGAACCGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATg                             >  1:178307/1‑91 (MQ=255)
                                                                  ccGTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTcc                   <  2:589356/91‑1 (MQ=255)
                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATcc                                   >  1:770782/1‑73 (MQ=255)
                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATcc                                   <  2:770782/73‑1 (MQ=255)
                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctg                 >  2:219688/1‑91 (MQ=255)
                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctg                 >  2:7804/1‑91 (MQ=255)
                                                                    gTTCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctg                 >  2:867594/1‑91 (MQ=255)
                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctgc                >  2:148653/1‑91 (MQ=255)
                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctgc                >  2:30564/1‑91 (MQ=255)
                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctgc                <  2:461031/91‑1 (MQ=255)
                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctgc                >  1:69008/1‑91 (MQ=255)
                                                                     ttCCGCATTCGTGTTATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCctgc                >  1:585855/1‑91 (MQ=255)
                                                                                   tATTGGGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGAt  >  2:309658/1‑91 (MQ=255)
                                                                                        |                                                                                     
CAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGAT  >  NZ_CP009273/3882051‑3882224

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: