Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,843 A→G 9.6% F406L (TTC→CTC)  yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,8430AG 9.6% 137.8 / 5.7 52F406L (TTC→CTC) yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (3/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

ACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTC  >  NC_002947/3893695‑3893985
                                                                                                                                                     |                                                                                                                                              
aCCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATAGACTAAAGGATGTCTGTCATGATCAAC‑AGCCCGGTGAAGTCACCGaa                                                                                                                                                >  8:233016/1‑149 (MQ=255)
         aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGa                                                                                                                                        <  5:148323/148‑1 (MQ=255)
         aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGAc                                                                                                                                       <  7:110468/149‑1 (MQ=255)
          acctCAGCGCACCGCAGACCACAAAGCCTTCCCCCTCCTGAACCATGCCCGGCCCTTCAACCTCACTGGCAACCGGCACCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGAGAAGTCACCGAACTCGTTGAc                                                                                                                                       >  6:280515/4‑148 (MQ=255)
                              acaAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACATCAATGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCTTGGTAAAC‑AGCCCGGTGAATTCACAGAACTCGTTGACCACGAAGGCAATATGGGTCgg                                                                                                                  >  6:17703/1‑148 (MQ=255)
                               caAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGATCGAc                                                                                                                 >  7:76728/1‑149 (MQ=255)
                                aaaGCCTTCCCCCTCCTGAACGATGCCCGGGCCATCAACCTCACTGGCATCCGGCAGCTAACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACt                                                                                                                >  6:322103/1‑149 (MQ=255)
                                     cTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcg                                                                                                           <  5:122843/149‑1 (MQ=255)
                                          cccTCCTGAACGATGCCCGGGCCTTCAACATCACTGGCATCCGGCAGCTCACAGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAt                                                                                                      >  4:10620/1‑149 (MQ=255)
                                           cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                     >  2:128930/1‑149 (MQ=255)
                                                tGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                >  3:140388/1‑149 (MQ=255)
                                                  aaCGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCa                                                                                              <  1:128930/149‑1 (MQ=255)
                                                    cGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGAcc                                                                                                                                      <  5:234414/107‑1 (MQ=255)
                                                    cGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGAcc                                                                                                                                      >  6:234414/1‑107 (MQ=255)
                                                    cGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTGCAgg                                                                                            >  7:233484/1‑149 (MQ=255)
                                                          ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATgg                                                                                                                       <  2:304354/116‑1 (MQ=255)
                                                          ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATgg                                                                                                                       >  1:304354/1‑116 (MQ=255)
                                                           ggggCCTTCATCCTCACTGGGATCCGGCAGCTCACCGTCAATCGACTCAAGTATGTCGGTCTTGGTCAAC‑AGCCGGGTGAAGTGACCGAACTCGTTGTGCACGATGGCTTTTTGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGa                                                                                     <  5:280515/148‑1 (MQ=255)
                                                            gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCgcgc                                                                                                        <  8:228256/129‑1 (MQ=255)
                                                            gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCgcgc                                                                                                        >  7:228256/1‑129 (MQ=255)
                                                              gCCGTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGCTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGa                                                                                  <  7:233016/149‑1 (MQ=255)
                                                                cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCCCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGGATTTGCTCCAGGGCATTGAGGATg                                                                                >  8:166159/1‑149 (MQ=255)
                                                                 ttCAACCTCAATGGCATCCGGCAGCTCACCCGCAATCGACTCAAGGATGTCGGTCATCCTCAAC‑AACCCGGTGAAGACACCGAACTCGTATACCACGAAGGCAATATGGGTCGACTGACCCCCCAATAGCTCCAGGGCATTGACGAtca                                                                               >  7:7463/1‑147 (MQ=255)
                                                                  tCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCt                                                                              >  1:132449/1‑149 (MQ=255)
                                                                   gaaCCTCACTGGCATCCGGGAGCTCACCGGCAATCGACTCAAGGATTTCGGTCATGGTCAAG‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGATCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCt                                                                              <  8:76728/147‑1 (MQ=255)
                                                                      ccTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAg                                                                          >  6:8633/1‑149 (MQ=255)
                                                                      ccTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAg                                                                          <  3:307024/149‑1 (MQ=255)
                                                                          aCTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCTCATTTGCTCCAGGGCATTGAGGATGCTGAAGctc                                                                       <  1:42353/148‑1 (MQ=255)
                                                                                          tCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCca                                                                   >  5:2083/1‑136 (MQ=255)
                                                                                          tCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCca                                                                   <  6:2083/136‑1 (MQ=255)
                                                                                                   aTCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGc                                             <  5:322103/149‑1 (MQ=255)
                                                                                                   aTCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCACGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGTTCGACTGACGGAGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGc                                             <  5:17703/149‑1 (MQ=255)
                                                                                                         tCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAg                                                                          >  4:301065/1‑114 (MQ=255)
                                                                                                         tCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAg                                                                          <  3:301065/114‑1 (MQ=255)
                                                                                                         tCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACAGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCt                                       >  4:307345/1‑149 (MQ=255)
                                                                                                               aTGTCGGTCAGCGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCa                                  <  8:284834/148‑1 (MQ=255)
                                                                                                                   cGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTc                              <  3:10620/148‑1 (MQ=255)
                                                                                                                   cGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCg                             >  6:279249/1‑149 (MQ=255)
                                                                                                                     ttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATg                                         >  7:22458/5‑135 (MQ=25)
                                                                                                                     ttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATg                                         <  8:22458/131‑1 (MQ=25)
                                                                                                                     gTCATTGTCAAC‑AGCCCTGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGAGCATTTGGTGCAGGGGTTTGTGGATGGTGAAGCTCTCCATCATGTTCAACGGCTCCCGGGCCATGCTTTCCTGGTCggg                           <  8:7463/149‑1 (MQ=255)
                                                                                                                     gTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                           <  2:132449/149‑1 (MQ=255)
                                                                                                                      tCATGGTCAAC‑AGCACGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggc                          >  2:117645/1‑149 (MQ=255)
                                                                                                                         tGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                       >  3:104686/1‑149 (MQ=38)
                                                                                                                         atttCAAC‑ATCCCGGTGAAGTAACCGAACTCGTTGACAACGAAGGCAAACTGGGTCGTCTGAACGCGCATTTGCTACAGGCCATTGACGATTCTGAAGCTCTCCATCATGATCAACGGCGCCCGGGACATGCTTTCCAGGTCGggctgg                       <  4:196325/146‑1 (MQ=255)
                                                                                                                               aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                 >  5:133377/1‑149 (MQ=255)
                                                                                                                                  aGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                 >  2:205344/1‑147 (MQ=255)
                                                                                                                                    cccGGTGAAGTCTCCGAACTCGGAGACCACGAAGGCAATATTGGTCGACTGACCGCGAATTTGGTCCAGGGCATTGATGATTCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga             <  3:206088/149‑1 (MQ=255)
                                                                                                                                       ttcgtggTCAACGAGTTCGGTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCGATGCTTTCCAGGTCGGGCTGGCTGCCCgaga             >  2:206278/7‑146 (MQ=255)
                                                                                                                                       ttcgtggTCAACGAGTTCGGTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga             <  1:206278/140‑1 (MQ=255)
                                                                                                                                       ggTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGCGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagca          <  1:205344/149‑1 (MQ=255)
                                                                                                                                           aaGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttc      <  5:8633/149‑1 (MQ=255)
                                                                                                                                                aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttc  <  5:279249/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
ACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTC  >  NC_002947/3893695‑3893985

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