Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,833,262 T→G V127G (GTC→GGC)  ascF → PTS cellobiose/arbutin/salicin transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,833,2620TG95.0% 58.5 / ‑3.3 20V127G (GTC→GGC) ascFPTS cellobiose/arbutin/salicin transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (9/10);  total (9/11)
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.27e-01

GTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGAC  >  NZ_CP009273/2833181‑2833347
                                                                                 |                                                                                     
gTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAActgc                                                                               >  2:485204/1‑90 (MQ=255)
 tCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAActgct                                                                              <  2:111650/90‑1 (MQ=255)
      ttGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCaa                                                                         >  1:448071/1‑90 (MQ=255)
        ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                       >  1:158818/1‑90 (MQ=255)
        ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                       >  1:340930/1‑90 (MQ=255)
        ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                       <  1:485204/90‑1 (MQ=255)
          tGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGAt                                                                     <  1:29141/90‑1 (MQ=255)
               ggATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCg                                                                <  2:254242/90‑1 (MQ=255)
                          gcgcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAgcggc                                                     >  2:471957/1‑90 (MQ=255)
                                  cGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGAc                                             >  2:326933/1‑90 (MQ=255)
                                         aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAgcgg                                                      <  1:475948/74‑1 (MQ=255)
                                         aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAgcgg                                                      >  2:475948/1‑74 (MQ=255)
                                         aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                      >  2:192202/1‑90 (MQ=255)
                                                        ccGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                       <  1:266824/90‑1 (MQ=255)
                                                        ccGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                       <  1:313828/90‑1 (MQ=255)
                                                        ccGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAAGATCGCCGACCTTAAcc                       >  1:253494/1‑90 (MQ=255)
                                                            cGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTc                   <  1:94806/90‑1 (MQ=255)
                                                                  tCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtg             <  2:274133/90‑1 (MQ=255)
                                                                   cggcggATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgt            <  1:312910/90‑1 (MQ=255)
                                                                             aTGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAGGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGAc  <  1:547830/90‑1 (MQ=255)
                                                                                 |                                                                                     
GTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGAC  >  NZ_CP009273/2833181‑2833347

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