Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,777,745 A→T 39.7% I312F (ATT→TTT)  Y75_RS18735 → PTS transporter subunit EIIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,777,7450AT39.7% 18.3 / 9.2 32I312F (ATT→TTT) Y75_RS18735PTS transporter subunit EIIC
Reads supporting (aligned to +/- strand):  ref base A (2/14);  new base T (4/5);  total (10/22)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCAACCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAT  >  NC_007779/3777508‑3777980
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
cacTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGAt                                                                                                                                                                                                                                                <  1:567165/239‑1 (MQ=255)
  cTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGGTGCTGttt                                                                                                                                                                                                                                               >  1:866511/1‑238 (MQ=255)
    tCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTccc                                                                                                                                                                                                                                            <  1:1786165/239‑1 (MQ=255)
                       cTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGCTTCTGTTTTCGCCCACCCTAAACCGCCaag                                                                                                                                                                                                                       >  1:1878958/1‑237 (MQ=255)
                       cTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTTGTGGTTTCCTCAACCCTAACCGcccagc                                                                                                                                                                                                                      >  1:1781728/1‑235 (MQ=255)
                          tGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTTCTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCGTTCTGTTTCCCGCCACCCCAACCCCCCTGCtgg                                                                                                                                                                                                                   >  2:686099/1‑239 (MQ=255)
                          tGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGTTGTTTTCCGCCACCCTAACCCCCCAGCtgt                                                                                                                                                                                                                   >  1:750591/1‑238 (MQ=255)
                          tGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCGGGGTCCCCCCCCCCTTACCGCCacacggg                                                                                                                                                                                                                   >  1:1981074/1‑233 (MQ=255)
                            gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTGTGATTTCCCCCACCCTAACACCCATggcgggg                                                                                                                                                                                                                 >  1:1179162/1‑233 (MQ=255)
                            gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTGGTTTCCTCCACCCTAACCGCCAAGGtgggg                                                                                                                                                                                                                 >  1:506075/1‑237 (MQ=255)
                                     cTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCg                                                                                                                                                                                                                                           <  1:1150275/207‑1 (MQ=255)
                                                  cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGGTGCTGATTTCCGCCACCTTAACCGCCATGCTGGTGGGGAATACCGAAACGCCGGGAt                                                                                                                                                                                          >  1:209787/1‑239 (MQ=255)
                                                                   cgcAGCATCTGAAAGAGTTCAGTTTGAGCGCCGAGCCGCTAAAATCGTTTTTCCCGGAGGGCGTTTTTGCCCTGCACGGTAACTAAAAAATCTTGGGGGCGGTGGGCTTTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGGGGGGGGTTTGCTGATTCCCGCAACCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCAt                                                                                                                                                                          <  1:1301239/238‑1 (MQ=255)
                                                                      agcaTCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGGTGCTTGTTTCCCCCACCCTAAACCCCCCACCGGGGGGGAGTTACCAACCGCTGGCATTTACCCTCCTTTTCAtttt                                                                                                                                                                       >  2:774417/1‑237 (MQ=255)
                                                                                    ttCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTccc                                                                                                                                                                                                                                            <  1:1688227/159‑1 (MQ=255)
                                                                                              gcgcCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTTGGTTCCCCCACCCTAAACCCCCTGCTCGTGGGGAGGATCGCAACGGCGGGTGTTACCCTTCCTTTTCTTTTCACCGTGCTTTTTTCgggtgtcgc                                                                                                                                               >  2:398282/1‑231 (MQ=255)
                                                                                                       cgcgaaaaTCGTTGTTCCCGGAGGGGGTTTTTGCCCTGCGCGGTACCTAAAAAATCTTGGGGGCGGTGGGCATTTCTTTGGCGATGTCCTTCACTCCCGCCCCGAAAAATCGGTTAAAAGGGGCGGCTTGGCTGATTCCCGCACCCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTggc                                                                                                                                      <  1:1290151/233‑1 (MQ=255)
                                                                                                                 gttgtcCCGGGAGGGGTGTTTCGCCCTGCACGGTAACAAAAAATTTTTGGGTCCGGTGGTTTTTTTTTTAGCGATTTCCTCCCCTGCCCACCCGGAAATGGGGGAAAAGGGGGGGGGTTGGTTGATCCCCGAACCTTAACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGt                                                                                                                           <  2:506075/239‑1 (MQ=255)
                                                                                                                  tttttCCCGGGAGGGGGTTTTTCCCCTGCGCGGTACATAAAAAATTTGGGGGCCGGTGGTTTTTTTTTTGGCGTTTTCCTCCCCTCCCCCCCGGAAAACCGGGTAAAAGGGGGGGGGTTTGTTGTTCCCCGCACCTTAAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGt                                                                                                                           <  2:1298424/236‑1 (MQ=255)
                                                                                                                         tggAGGGCGGTTTTGCCGTGCAGGTTAACTAAAAATTTGTGGTGGCGGTGGGTTTTTCTTAGGGGATTTCCTTCCCTCCGGCACGGAAAAGTCGGGAAAGGGGGGGGGGTTGGCTGATCCCCCCACCCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTa                                                                                                                   <  2:322003/238‑1 (MQ=255)
                                                                                                                                 ggTTTTGCCCGTCGCGGTAACCAAAAATTTTTTGGTGGCGGGGGCTTTTTGTTGGGGGGGTACTTTCACGGCCGCACGGAAAAGCGGGGTAAAGGGGGGGGGTTGGTGGTTTCCGGAACCTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTAt                                                                                                             <  2:1298396/238‑1 (MQ=255)
                                                                                                                                                                             gggCTTTTCTTTAGGGTTGTCCTTCCCTCCCGCCCCGGAAAATCGGGTAAATGGGGGGGGTTTCTTATTTCCCGCACCTTAACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTCCCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGAc                                                                   <  1:774417/238‑1 (MQ=255)
                                                                                                                                                                                               gTACTTAACGCCCCCCCCGGAAAATCGGTAAAAAGGGGGGGCTTTGCTGATCCCAGCACCTTTAACCGCCATGCTGGTGGGAATTACCGAGCCGCTGGAATTTGCCTTCCTGTTCATTTCCCCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACt                                                <  1:1318119/239‑1 (MQ=255)
                                                                                                                                                                                                        ccgCCGCACCAGAAATCCGGGTAAAGGGGGGGGCTTGGTTGCTTCCCACATCCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACCGGATCCCg                                        <  2:548133/237‑1 (MQ=255)
                                                                                                                                                                                                         tGCCGCACCGAAAATTGGGGAAAAAGGGGGGGGTTTGTTGCTCCCACAACCTTAATCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTGCCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGa                                       <  1:1054601/239‑1 (MQ=255)
                                                                                                                                                                                                                       aatcggTAAAAGGGGGGGGCTTTGTTATTCCCCCAAACTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCa                          <  2:1179162/235‑1 (MQ=255)
                                                                                                                                                                                                                           ggTAAAAGGGGGGGCGTTTCTGTATCCCCGCACCTTTACCCGCCATGCTGGTGGGAATTCCCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcgc                       <  1:1063401/238‑1 (MQ=255)
                                                                                                                                                                                                                              aaatGTGGGGGGCTTGCTGTTCCCCGCACCTTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGata                  <  1:25485/239‑1 (MQ=255)
                                                                                                                                                                                                                              aaagGTGGGGGGTTTGCTGTTTCCCGAACCCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGata                  <  1:231183/239‑1 (MQ=255)
                                                                                                                                                                                                                                      ggcgtgttagCTTATCCCCGCAACCTTACCCGCCATGCTGGTGGGAATTACCGACCCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATatgatg             <  1:1079359/229‑1 (MQ=255)
                                                                                                                                                                                                                                          ctttGCTGTTCCCCCCAACCTCACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCTTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGAc        <  1:980106/236‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGATCCCCGCATCCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAt  <  1:203341/238‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGATCCCCCCACCTTAAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAt  <  1:398282/238‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
CACTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCAACCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAT  >  NC_007779/3777508‑3777980

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