Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,507,458 T→C N146S (AAC→AGC)  fryC ← putative PTS enzyme IIC component FryC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,507,4580TC98.2% 184.5 / ‑4.7 55N146S (AAC→AGC) fryCputative PTS enzyme IIC component FryC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (51/3);  total (52/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ACAGGAAATCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGG  >  NC_000913/2507321‑2507593
                                                                                                                                         |                                                                                                                                        
aCAGGAAATCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTc                                                                                                                                        >  2:1247821/1‑140 (MQ=255)
        tCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTg                                                                                                                                >  1:1471622/1‑140 (MQ=255)
           ggTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATCGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGt                                                                                                                             <  2:1471622/140‑1 (MQ=255)
                  ggATTTTCTGGTTCATCCACCTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTgg                                                                                                                      <  2:556384/140‑1 (MQ=255)
                    aTTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGt                                                                                                                    <  2:768977/140‑1 (MQ=255)
                           ggTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATc                                                                                                             >  1:590013/1‑140 (MQ=255)
                              tCATCCACTTCACAAGGTAACCGGTAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAAc                                                                                                          >  1:992898/1‑140 (MQ=255)
                               cATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAAcc                                                                                                         >  1:1490953/1‑140 (MQ=255)
                                  ccACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACccgc                                                                                                      >  1:1325227/1‑140 (MQ=255)
                                       tCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATa                                                                                                 >  1:217758/1‑140 (MQ=255)
                                             ggTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCa                                                                                           >  1:1057908/1‑140 (MQ=255)
                                                aaCCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGcc                                                                                        >  2:1430598/1‑140 (MQ=255)
                                                  ccGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCgg                                                                                      >  1:577940/1‑140 (MQ=255)
                                                       aaCGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACg                                                                                 >  1:403210/1‑140 (MQ=255)
                                                         cGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGcc                                                                               >  2:968449/1‑140 (MQ=255)
                                                                aaaaatCAGCGCACAAAGGAGGGGGGCGGGCCCGGGAGAGGGGGCCCCCCTTTGAATGGTGGTGGGATCGAACCTCAGCGGTTGCGGGGGCGGGGAAACCACCACCCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGc                                                                        <  2:1117774/137‑1 (MQ=255)
                                                                aaaaatCAGCCCACAAGGGAGGGGGGAGGGCCCGGGGGAAGGGGCCGCCTTTTGTATGGTGGGGGAATCGAAGCTCGGCGTTTGCGGGGGGGGGGAAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGc                                                                        <  2:577940/137‑1 (MQ=255)
                                                                    atCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCAcc                                                                    >  2:597777/1‑140 (MQ=255)
                                                                    atCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCAcc                                                                    >  1:791092/1‑140 (MQ=255)
                                                                    atCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCAcc                                                                    >  2:182074/1‑140 (MQ=255)
                                                                      cAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGa                                                                  >  1:103784/1‑140 (MQ=255)
                                                                         cgcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATag                                                               >  1:1117774/1‑140 (MQ=255)
                                                                         cgcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATag                                                               >  1:1529902/1‑140 (MQ=255)
                                                                         cgcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATag                                                               >  2:954482/1‑140 (MQ=255)
                                                                         cgcACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATag                                                               >  1:903605/1‑140 (MQ=255)
                                                                         cgcACCAATGAAGGTGGACGGAACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATag                                                               >  2:1179118/1‑140 (MQ=255)
                                                                           cACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATagag                                                             >  1:334045/1‑140 (MQ=255)
                                                                              cAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTc                                                          >  1:1126416/1‑140 (MQ=255)
                                                                                aTGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGc                                                        >  1:4624/1‑140 (MQ=255)
                                                                                     ggTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACaa                                                   >  2:1047467/1‑140 (MQ=255)
                                                                                        ggACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAag                                                >  1:869494/1‑140 (MQ=255)
                                                                                        ggACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAag                                                >  1:818816/1‑140 (MQ=255)
                                                                                        ggACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAag                                                >  1:1578410/1‑140 (MQ=255)
                                                                                         gACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAagc                                               >  1:743279/1‑140 (MQ=255)
                                                                                          aCGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAagca                                              >  2:975506/1‑140 (MQ=255)
                                                                                          aCGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAagca                                              >  2:756807/1‑140 (MQ=255)
                                                                                              cACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCaa                                          >  1:1354026/1‑140 (MQ=255)
                                                                                              cACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCaa                                          >  1:1369377/1‑140 (MQ=255)
                                                                                               aCCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCaaa                                         >  2:1447379/1‑140 (MQ=255)
                                                                                                             cGCCCATTGCATAGCGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCaa                          >  2:1579201/1‑140 (MQ=255)
                                                                                                                          gTGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAacac             >  1:1146051/1‑140 (MQ=255)
                                                                                                                           tGCTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcg            >  1:1063428/1‑140 (MQ=255)
                                                                                                                             cTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcgcc          >  1:280000/1‑140 (MQ=255)
                                                                                                                             cTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAACAACCAACATCCGGAATGGAAAACCCAAAAAAAcgccc         >  2:1214180/1‑139 (MQ=255)
                                                                                                                             cTGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcgccg         >  1:55187/1‑140 (MQ=255)
                                                                                                                              tGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcgccg         >  2:1454981/1‑140 (MQ=255)
                                                                                                                              tGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcgccg         >  2:160294/1‑140 (MQ=255)
                                                                                                                              tGGGATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACAcgccg         >  1:1445523/1‑140 (MQ=255)
                                                                                                                                ggATCGAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACgccgcc       >  2:558327/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  2:27263/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  2:789512/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  1:1505303/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  1:36743/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  2:1216477/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  2:1629483/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  1:157928/1‑140 (MQ=255)
                                                                                                                                     gAAGCTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAgg  >  2:1376322/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
ACAGGAAATCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGTTCAGCAGTTGCGTCGGCTGGGTAGACATCAACCCGCCGATAAAACCAGCCGGGAACGCCAGTTTGCCACCGATAGAGTTCGCCACAAAAGCAGCAAACATCGGAATGGCAAAACCAAACAACACGCCGCCGAAGG  >  NC_000913/2507321‑2507593

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