Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,528,938 T→A L505H (CTT→CAT)  ptsI → PEP‑protein phosphotransferase of PTS system (enzyme I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,528,9380TA97.0% 99.5 / NA 33L505H (CTT→CAT) ptsIPEP‑protein phosphotransferase of PTS system (enzyme I)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (11/21);  minor base G (1/0);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACACTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAG  >  CP009273/2528798‑2529074
                                                                                                                                            |                                                                                                                                        
gTACCCTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATg                                                                                                                                        <  1:362329/143‑1 (MQ=255)
gTACACTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGCTCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTgggt                                                                                                                                  >  2:29302/1‑147 (MQ=255)
               tGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTg                                                                                                                         <  1:317851/143‑1 (MQ=255)
                     tGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTacac                                                                                                                   >  1:560696/1‑143 (MQ=255)
                               atgatTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCt                                                                                                         >  1:1810/1‑143 (MQ=255)
                                tgatTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTg                                                                                                        <  1:209534/143‑1 (MQ=255)
                                         acCTTTACCAGCCAATGTCACCGCCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGCCTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGt                                                                                               <  1:462545/143‑1 (MQ=255)
                                           cTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGAGGGGTCt                                                                                             >  1:50362/1‑143 (MQ=255)
                                                 ccgCCAATGTCACCGTCCGTGCTGAACTTGCTCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGAAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCCGTTGCTGGGGATGGGTCTGGACGa                                                                                       <  1:34870/141‑1 (MQ=255)
                                                 cAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATtctc                                                                                 <  2:331160/149‑1 (MQ=255)
                                                    ccAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAAtt                                                                                    <  1:559310/143‑1 (MQ=255)
                                                      aaTGTCGCCGTCCGTGCTGAACTTGATCAAGGAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACCCTTCTGTTGCTGGGGATGGGTCTGGACGAATtct                                                                                  <  1:55787/143‑1 (MQ=255)
                                                          tCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCGTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTAt                                                                              >  1:1021/1‑143 (MQ=255)
                                                                  cgtGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGTGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCtc                                                                >  2:700290/1‑148 (MQ=255)
                                                                      gtGAACTTGATCAAGCAAGTTATTGATGCTTCTCTTGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAACTCTCTATGAGCGCCATTTCTATccc                                                            <  2:354379/148‑1 (MQ=255)
                                                                      cTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGACAAATGGACTGGCATGTGTGGTGAGCATGGTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATccc                                                            >  2:665766/1‑149 (MQ=255)
                                                                       tGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCcg                                                           >  2:256346/1‑149 (MQ=255)
                                                                                     caagTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTaa                                                   <  1:489428/143‑1 (MQ=255)
                                                                                        gTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTAt                                          <  2:2225/149‑1 (MQ=255)
                                                                                           aTTGCTGCTTCTCATGCTGAAGGCAAATGGACTGGCCTGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGAt                                             <  1:246361/143‑1 (MQ=255)
                                                                                           aTTGATGGTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGGGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGAt                                             <  1:246317/143‑1 (MQ=255)
                                                                                           aTTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCg                                       <  2:453500/149‑1 (MQ=255)
                                                                                                    tctcATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTGTGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGaa                              >  2:386399/1‑149 (MQ=255)
                                                                                                      tcATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAaca                                  >  1:487603/1‑143 (MQ=255)
                                                                                                       caggCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAacac                                 <  1:108382/140‑1 (MQ=255)
                                                                                                                 gCAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGaa                       <  1:176610/143‑1 (MQ=255)
                                                                                                                  cAAATGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGa                <  2:274085/149‑1 (MQ=255)
                                                                                                                   aaaTGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGGTGCGaa               >  2:511364/1‑149 (MQ=255)
                                                                                                                   aaaTGGACTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGaa               >  2:18831/1‑149 (MQ=255)
                                                                                                                      tGGACTGGCATGTGTGGTGGGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCCATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGgt            <  2:55969/149‑1 (MQ=255)
                                                                                                                          cTGGGATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGGCTGGACGAATTCGCTATGAGCTCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAg              <  1:572218/143‑1 (MQ=255)
                                                                                                                          cTGGCATGTGTGGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAg              <  1:341770/143‑1 (MQ=255)
                                                                                                                                tgtgtgGTGAGCATGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTagcagag  <  2:591676/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
GTACACTCTGGCAGTTGACCGTGGTAATGATATGATTTCACACCTTTACCAGCCAATGTCACCGTCCGTGCTGAACTTGATCAAGCAAGTTATTGATGCTTCTCATGCTGAAGGCAAATGGACTGGCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAG  >  CP009273/2528798‑2529074

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