Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 161,062 G→A 7.3% A142V (GCG→GTG)  sfsA ← sugar fermentation stimulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913161,0620GA 7.3% 163.4 / 2.8 55A142V (GCG→GTG) sfsAsugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base A (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GGGTGAAAACCGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGAAATTGATTCATTAAGGA  >  NC_000913/160917‑161209
                                                                                                                                                 |                                                                                                                                                   
gggTGAAAACCGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTa                                                                                                                                                  <  4:150068/149‑1 (MQ=255)
gggTGAAAACCGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGATACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCt                                                                                                                                                   >  6:176858/1‑148 (MQ=255)
  gTGAAAACCGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAAc                                                                                                                                                <  8:111048/149‑1 (MQ=255)
       aaCCGTGTATTGGCTGTATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCGTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGGTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTGCGCTAACGTAAc                                                                                                                                           <  4:95476/149‑1 (MQ=255)
          cGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGa                                                                                                                                        >  6:2303/1‑149 (MQ=255)
                aTGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCAc                                                                                                                                  >  6:32759/1‑149 (MQ=255)
                aTGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGATACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCa                                                                                                                                   <  5:176858/148‑1 (MQ=255)
                       aaTGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCaa                                                                                                                             <  6:223689/147‑1 (MQ=255)
                       aaTGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAata                                                                                                                           >  2:166972/1‑149 (MQ=255)
                       aaTGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAata                                                                                                                           >  4:17228/1‑149 (MQ=255)
                            agcaCGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCa                                                                                                                      <  4:146137/149‑1 (MQ=255)
                            agcaCGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGG‑AAATATCCCTGTTCGTTCTCCGCTAACGTAACCGCTTTCACTGCAATATAGCAg                                                                                                                     >  3:168818/1‑149 (MQ=255)
                            agcaCGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGGTCAGTGACCGCATC‑GGTAAATATCCATGTTCGTTCTACGCTAACGTAACCGATTTCACTTCAATATAGCAg                                                                                                                     >  3:165726/1‑149 (MQ=255)
                                 ggCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTg                                                                                                                 <  8:1435/149‑1 (MQ=255)
                                   cGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGa                                                                                                               <  7:204865/149‑1 (MQ=255)
                                   cGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGGGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGa                                                                                                               >  4:197777/1‑149 (MQ=255)
                                    gAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTTACCGCATCGGGAAAAGATCACTGTTCGTTCTCCGCTAACGTAACCGATTTAACTTCAATATAGCAGTCTGGAc                                                                                                              >  4:137387/1‑149 (MQ=255)
                                      aaaaGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTGCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGAcgc                                                                                                            <  3:144808/149‑1 (MQ=255)
                                            tAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATcc                                                                                                      <  7:226141/149‑1 (MQ=255)
                                                         gCCTTCTTCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGCCCTCGTTCAGTGACCGCATCGGG‑AAATATCCCTTTTCTTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATaaa                                                                                        <  4:165726/149‑1 (MQ=255)
                                                         gCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACAt                                                                                           <  5:2303/147‑1 (MQ=255)
                                                         gCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATaa                                                                                         <  6:178455/149‑1 (MQ=255)
                                                         gCCTGCAGCGGCTACGCTGATCTACTCCCGTAGGTGTTTCTGAGCTCGTTCAGTGACCGCATCGGG‑AATTATCCCTGTGCGTTGTCCGCTAACGTAACGGTTTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATaaa                                                                                        <  4:168818/149‑1 (MQ=255)
                                                          ccTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATaaa                                                                                        <  3:17228/149‑1 (MQ=255)
                                                            ttCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGt                                                                                      >  7:70990/1‑149 (MQ=255)
                                                                            atcaACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCACTAACGTAACCGATTTCACTTCAATATAg                                                                                                                        <  4:193913/99‑1 (MQ=255)
                                                                            atcaACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCACTAACGTAACCGATTTCACTTCAATATAg                                                                                                                        >  3:193913/1‑99 (MQ=255)
                                                                                    ccGAAGGTGTTTCTGCCCTCGTTCAGTGACCGCAGCGGGAAAATATCCCGGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGcc                                                                <  3:47581/147‑1 (MQ=255)
                                                                                    ccGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCg                                                               >  8:109336/1‑148 (MQ=255)
                                                                                    ccGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCg                                                               >  6:70956/1‑148 (MQ=255)
                                                                                     cGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTa                                                             <  1:166972/149‑1 (MQ=255)
                                                                                     cGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAACATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTa                                                             >  2:31966/1‑149 (MQ=255)
                                                                                     cGAAGGTGTTTCTGACCTCGTTCAGTGACAGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTa                                                             >  7:227048/1‑149 (MQ=255)
                                                                                        aGGTGTTTCTGACCTCGTTCAGTGACAGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTAttt                                                          <  8:227048/149‑1 (MQ=255)
                                                                                                cTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCg                                                   <  6:94580/148‑1 (MQ=255)
                                                                                                cTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGc                                                  >  5:156947/1‑149 (MQ=255)
                                                                                                 tGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCt                                                 >  7:140948/1‑149 (MQ=255)
                                                                                                 tGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCACTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCt                                                 >  3:227628/1‑149 (MQ=255)
                                                                                                   aCCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCttt                                               >  5:148175/1‑149 (MQ=255)
                                                                                                        gTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCa                                            <  5:32759/147‑1 (MQ=255)
                                                                                                        gTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGc                                          >  7:52094/1‑149 (MQ=255)
                                                                                                                aCCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACTTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCtata                                  <  5:201686/149‑1 (MQ=255)
                                                                                                                    cATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCt                              <  3:197777/149‑1 (MQ=255)
                                                                                                                     aTCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCt                              >  4:101339/1‑148 (MQ=255)
                                                                                                                        gggAAAATATCCCTGTTCTTTCACCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACa                          <  8:46193/149‑1 (MQ=255)
                                                                                                                        gggAAAATATCCCTGTTCGTTCTCCACTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGAc                           <  4:227628/148‑1 (MQ=255)
                                                                                                                          gAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGt                        >  3:51608/1‑149 (MQ=255)
                                                                                                                              atatCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCt                     >  1:83701/1‑148 (MQ=255)
                                                                                                                               tatCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTg                    >  2:162826/1‑148 (MQ=255)
                                                                                                                                 tCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGaaa                 <  7:109336/149‑1 (MQ=255)
                                                                                                                                          gttcTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACgg                                                                              >  7:216640/1‑79 (MQ=255)
                                                                                                                                          gttcTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACgg                                                                              <  8:216640/79‑1 (MQ=255)
                                                                                                                                               ccGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGAAATTGATTCATTAAg    >  3:112972/1‑148 (MQ=255)
                                                                                                                                                cGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGAAATTGATTCATTAAgg   <  4:51608/148‑1 (MQ=255)
                                                                                                                                                cGCTAACGTAACCGATTTCACTTCAAGATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGAAATTGATTCATTAAGGa  <  1:31966/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGGTGAAAACCGTGTAATGGCTGAATGCAGCACGGCGAAAAAGATAACCGCACGCTGGCCTTCAGCCGCTACGCTCATCAACTCCCGAAGGTGTTTCTGACCTCGTTCAGTGACCGCATCGGGAAAATATCCCTGTTCGTTCTCCGCTAACGTAACCGATTTCACTTCAATATAGCAGTCTGGACGCGAATCCGCCTGCAACATAAAGTCAATACGGCTGCGTTCTGCGCCGTATTTTACTTCGCTTTTCAGCGAGCTATAGCCTGACAGTTCTGAAATTGATTCATTAAGGA  >  NC_000913/160917‑161209

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