Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,065,206 A→T L269Q (CTG→CAG)  yihT ← sulfofructosephosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,065,2060AT95.0% 138.9 / ‑0.4 40L269Q (CTG→CAG) yihTsulfofructosephosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (16/22);  total (17/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 = 8.08e-01

TATCCTCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCAGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGCGCCCGCCGTCATTGCC  >  NZ_CP009273/4065123‑4065292
                                                                                   |                                                                                      
tATCCTCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCa                                                                                  <  1:472663/90‑1 (MQ=255)
     tCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCg                                                                             >  2:645608/1‑90 (MQ=255)
     tCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCg                                                                             >  2:521367/1‑90 (MQ=255)
     tCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCg                                                                             >  1:627273/1‑90 (MQ=255)
      cATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGt                                                                            >  2:872397/1‑90 (MQ=255)
        tCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTg                                                                          <  2:384992/90‑1 (MQ=255)
        tCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTg                                                                          <  2:165007/90‑1 (MQ=255)
          ttAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTc                                                                        <  1:332364/90‑1 (MQ=255)
           tAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCAGCATCAGCTCGTTGTCt                                                                       >  1:429716/1‑90 (MQ=255)
           tAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGCCACGCTGCATCAGCTCGTTGTCt                                                                       >  2:919805/1‑90 (MQ=255)
               ggcgTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTa                                                                   <  1:958349/90‑1 (MQ=255)
                gcgTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTaa                                                                  >  1:221751/1‑90 (MQ=255)
                  gTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAAc                                                                >  1:325639/1‑90 (MQ=255)
                  gTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAAc                                                                >  2:266400/1‑90 (MQ=255)
                       gCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACcgacg                                                           <  1:182106/90‑1 (MQ=255)
                       gCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACcgacg                                                           <  2:1941/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  1:363514/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  1:603005/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  1:322829/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  1:737967/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  1:919805/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  2:811970/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  2:803928/90‑1 (MQ=255)
                           tcatTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCAGCATCAGCTCGTTGTCTGGTAAACCGACGACCg                                                       <  2:347719/90‑1 (MQ=255)
                                   tcgACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAATGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAg                                               <  1:620148/90‑1 (MQ=255)
                                   tcgACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAg                                               <  2:102339/90‑1 (MQ=255)
                                           atCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACg                                       <  2:540058/90‑1 (MQ=255)
                                            tCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACgg                                      >  1:131429/1‑90 (MQ=255)
                                                   gttgttGTAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCa                               >  2:713492/1‑90 (MQ=255)
                                                     tgttgtAATTTCGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAgg                             <  1:122277/90‑1 (MQ=255)
                                                                cGGTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATgc                  >  2:59707/1‑90 (MQ=255)
                                                                 ggTGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATgcg                 <  2:5826/90‑1 (MQ=255)
                                                                   tGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATgc                  <  2:310540/87‑1 (MQ=255)
                                                                   tGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATgc                  >  1:310540/1‑87 (MQ=255)
                                                                   tGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGcgcc               >  2:422636/1‑90 (MQ=255)
                                                                   tGCGCAAACGTCACGCTGCATCAGCTCGTTGTCTAGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGcgcc               >  2:604061/1‑90 (MQ=255)
                                                                    gcgcAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGcgccc              >  1:218395/1‑90 (MQ=255)
                                                                    gcgcAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGcgccc              <  2:394556/90‑1 (MQ=255)
                                                                      gcAAACGTCACGCTGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGcgcccgc            >  2:381889/1‑90 (MQ=255)
                                                                                cGCTGCATCAGCTCGTTGTCTGGTAACCCGACGCCCGATGCCCAGACTGCCCGGCCTGCCAGGAATCCCGATGCGCCCGCCGTCATTGcc  <  1:563076/90‑1 (MQ=255)
                                                                                   |                                                                                      
TATCCTCATCTTAGCGGCGTTTAGCCATCATTTCGTCGACGATATCGCCAAGTTGTTGTAATTTCGGTGCGCAAACGTCACGCAGCATCAGCTCGTTGTCTGGTAAACCGACGACCGATGCCCAGACTGCACGGCCTGCCAGGAATCCCGATGCGCCCGCCGTCATTGCC  >  NZ_CP009273/4065123‑4065292

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