Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,940,705 G→A A208V (GCG→GTG)  yieP ← putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,940,7050GA98.2% 224.2 / NA 57A208V (GCG→GTG) yiePputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/25);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTCATGC  >  NC_000913/3940558‑3940848
                                                                                                                                                   |                                                                                                                                               
ccATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACAt                                                                                                                                              >  2:288007/1‑151 (MQ=255)
ccATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACAt                                                                                                                                              >  2:290067/1‑151 (MQ=255)
 cATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATc                                                                                                                                             >  2:454347/1‑151 (MQ=255)
     gccgcGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGAc                                                                                                                                          >  2:591872/1‑150 (MQ=255)
       cgcgATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAAt                                                                                                                                       >  2:440688/1‑151 (MQ=255)
            cccACGCCAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCt                                                                                                                                  <  1:300072/150‑1 (MQ=255)
                 gacAACCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGc                                                                                                                             <  2:279149/149‑1 (MQ=255)
                    aaccccGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGAt                                                                                                                               <  2:88781/143‑1 (MQ=255)
                    aaaCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATg                                                                                                                              <  1:465929/147‑1 (MQ=255)
                    aaaCCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGc                                                                                                                             >  1:483913/1‑148 (MQ=255)
                       cccGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATc                                                                                                                        <  1:113395/150‑1 (MQ=255)
                                     cttcttTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACtg                                                                                                          <  2:189313/150‑1 (MQ=255)
                                       tcttTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACtgt                                                                                                         <  1:48560/149‑1 (MQ=255)
                                         tttttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCg                                                                                                     >  2:108888/1‑151 (MQ=255)
                                            ttATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTg                                                                                                  >  2:609774/1‑151 (MQ=255)
                                              aTCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGt                                                                                                >  2:155096/1‑151 (MQ=255)
                                                cGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTa                                                                                               >  2:588787/1‑150 (MQ=255)
                                                cGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTa                                                                                               >  2:20486/1‑150 (MQ=255)
                                                cGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGAT‑ATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAAt                                                                                             <  1:403119/151‑1 (MQ=255)
                                                      ttCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACg                                                                                        >  2:132764/1‑151 (MQ=255)
                                                        cAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTg                                                                                      >  2:136152/1‑151 (MQ=255)
                                                        cAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTg                                                                                      >  1:201717/1‑151 (MQ=255)
                                                                  ttATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAgtgt                                                                            >  2:414722/1‑151 (MQ=255)
                                                                    aTCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGtgtg                                                                           >  2:431945/1‑150 (MQ=255)
                                                                    aTCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGa                                                                          <  1:643508/151‑1 (MQ=255)
                                                                      cGGGGTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGAt                                                                         >  2:175324/1‑150 (MQ=255)
                                                                       gggTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATa                                                                        >  1:307145/1‑150 (MQ=255)
                                                                                  tATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGa                                                             <  1:529222/150‑1 (MQ=255)
                                                                                   aTCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGaa                                                            >  2:560699/1‑150 (MQ=255)
                                                                                     cAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATa                                                          <  2:372134/150‑1 (MQ=255)
                                                                                     cAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAg                                                         <  1:40228/151‑1 (MQ=255)
                                                                                       ggTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCg                                                       <  2:389930/151‑1 (MQ=255)
                                                                                          gAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAgg                                                    >  1:306064/1‑151 (MQ=255)
                                                                                             cGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCg                                                  >  2:504521/1‑150 (MQ=255)
                                                                                              gTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGa                                                 >  1:198849/1‑150 (MQ=255)
                                                                                                agcagcGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATg                                              >  1:1792/1‑151 (MQ=255)
                                                                                                   agcGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAg                                            <  1:669762/150‑1 (MQ=255)
                                                                                                     cgcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTc                                         <  1:559893/151‑1 (MQ=255)
                                                                                                           aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAAc                                   <  2:73393/151‑1 (MQ=255)
                                                                                                           aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAAc                                   <  1:97466/151‑1 (MQ=255)
                                                                                                           aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAAc                                   <  1:667886/151‑1 (MQ=255)
                                                                                                                 tttAAATGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTg                             >  2:645367/1‑151 (MQ=255)
                                                                                                                    aaaTGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgcg                          >  1:286397/1‑151 (MQ=255)
                                                                                                                      aTGCCGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAGTGGCGCt                        <  1:400474/151‑1 (MQ=255)
                                                                                                                         ccGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt                     <  2:299416/151‑1 (MQ=255)
                                                                                                                         ccGCGTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt                     >  1:44656/1‑151 (MQ=255)
                                                                                                                            cgTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTc                  >  1:355304/1‑151 (MQ=255)
                                                                                                                             gTCGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTc                  >  2:244644/1‑150 (MQ=255)
                                                                                                                               cGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa               <  2:683184/151‑1 (MQ=255)
                                                                                                                               cGCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa               <  2:310859/151‑1 (MQ=255)
                                                                                                                                gCCATCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGCCGCTCATTTCAt                >  2:48056/1‑149 (MQ=255)
                                                                                                                                   aTCGCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGata           <  1:83806/151‑1 (MQ=255)
                                                                                                                                      gCTTTGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGc        >  2:44896/1‑151 (MQ=255)
                                                                                                                                        tttGGATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTc      >  2:149236/1‑151 (MQ=255)
                                                                                                                                           ggATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTCAt    <  2:424927/150‑1 (MQ=255)
                                                                                                                                           ggATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTCATg   <  1:431931/151‑1 (MQ=255)
                                                                                                                                            gATAATCACATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTCATGc  >  1:138672/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCATCGCCGCGATCCACGGCAAACCCGCCATACTGCGCTTCTTTTTATCGCTCATTCAATGTCCTGTTATCGGGTTATCACTTATCAGGTGAGCGTAGCAGCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATATGCTCATGC  >  NC_000913/3940558‑3940848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑