Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,509,859 G→A S77N (AGC→AAC)  yfeO → putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,8590GA82.1% 150.8 / ‑3.0 56S77N (AGC→AAC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/20);  minor base C (0/10);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 9.93e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCG  >  NC_000913/2509716‑2510005
                                                                                                                                               |                                                                                                                                                  
ggtgATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGt                                                                                                                                             <  1:660780/151‑1 (MQ=255)
   gatgaAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCa                                                                                                                                           <  1:278064/150‑1 (MQ=255)
    atgaAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATg                                                                                                                                         <  1:143795/151‑1 (MQ=255)
       aaaaTCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCg                                                                                                                                      <  1:12692/151‑1 (MQ=255)
             gCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGAcc                                                                                                                                  >  1:604450/1‑149 (MQ=255)
                 cGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGAccc                                                                                                                            <  1:131911/151‑1 (MQ=255)
                  ggTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGAcccc                                                                                                                           >  1:460354/1‑151 (MQ=255)
                  ggTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCATCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTCATCCGTTTCAACCAGGGTCATGGCGGTCCAGAccc                                                                                                                            >  2:623119/1‑150 (MQ=255)
                   gTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCg                                                                                                                          >  2:653878/1‑151 (MQ=255)
                       tACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCtg                                                                                                                      >  1:507280/1‑151 (MQ=255)
                                   tctGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGAt                                                                                                          >  1:637223/1‑151 (MQ=255)
                                      ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCgg                                                                                                       >  2:456107/1‑151 (MQ=255)
                                      ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCgg                                                                                                       >  1:787758/1‑151 (MQ=255)
                                         aaCGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTGGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGgc                                                                                                      <  1:497467/149‑1 (MQ=255)
                                          aCGACTGCCGGGAACTCTGGGGATACCCCAGGATTCACCCCTCGGGATCATCGGTGTATTAACGCTACCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCa                                                                                                   <  1:805877/151‑1 (MQ=255)
                                           cGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCa                                                                                                   <  1:446480/150‑1 (MQ=255)
                                             aCTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCg                                                                                                >  1:642617/1‑151 (MQ=255)
                                                gCCGGGAACTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGTTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCgg                                                                                               <  1:567664/149‑1 (MQ=255)
                                                 ccGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTc                                                                                            >  1:553773/1‑151 (MQ=255)
                                                 ccGGAACCTCTGGGGATAGCCCAGGATTCCCCCCTCGGGATCATCGGTTTATTAACGCTAACGGGTATTGGGGGGGGGTTGGTTCTCCGTTTCAACCAGGGCCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCgg                                                                                               <  2:460307/148‑1 (MQ=255)
                                                 ccGGAAACTCTGGGGATAGCCCAGGTTTCACCCCTCTGAATCAGCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCACCCGGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTc                                                                                            <  2:445944/151‑1 (MQ=255)
                                                      aaCTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGTTTGGTTATCCGTTTCACCCAGGGTCAGGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGcc                                                                                        <  2:764712/150‑1 (MQ=255)
                                                      aaCTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGccc                                                                                       >  2:179895/1‑151 (MQ=255)
                                                        ctctGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctc                                                                                     >  1:25025/1‑151 (MQ=255)
                                                          ctGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTg                                                                                   >  2:757399/1‑151 (MQ=255)
                                                            ggggATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctct                                                                                    >  2:644480/1‑148 (MQ=255)
                                                                     ccAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTgg                                                                         >  1:515200/1‑150 (MQ=255)
                                                                       aGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACt                                                                      >  1:376579/1‑151 (MQ=255)
                                                                        ggtttCCCCCCTCTGGATCATCGGTGTATTAACGCTACCGGTTATTGGGGTGGGGTTGTTTATCCGTTTCACCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACtt                                                                     <  2:369379/148‑1 (MQ=255)
                                                                        ggATTCCCCCCTCTGGATCATCGGTTTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACt                                                                      <  2:25003/150‑1 (MQ=255)
                                                                                cccTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGc                                                             >  1:691275/1‑151 (MQ=255)
                                                                                      ggATCATCGGTGTTTTAACGCTACCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCACCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAt                                                           <  1:904636/147‑1 (MQ=255)
                                                                                        atcatcGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATtc                                                     >  1:658199/1‑151 (MQ=255)
                                                                                              ggTTTATTAACGCTAACGGGTATTGGGGGGGGGTTGGTTACCCGTTTCACCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTc                                               <  2:505584/151‑1 (MQ=255)
                                                                                              ggTGTATTAACGCTAACGGGTATTGCGGTGGGTTTGGTTATCCTTTTCAACCAGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTc                                               <  1:266528/151‑1 (MQ=255)
                                                                                                        cGCTAGCGGGTATTGCGGTGGGCTTGTTGACCCTTTTCACCCAGGGTCATGCCGCACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCGCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcg                                      <  2:835084/150‑1 (MQ=255)
                                                                                                            aaggggTATGGCGGGGGGGTTGTTTATCCGTTTCAACCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcggc                                    <  2:226946/145‑1 (MQ=255)
                                                                                                            aaCGGGTATGGGGGTGGGGTTGGTTAGCCGTTACACCAAGGGTAAGGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCt                                          <  1:204429/142‑1 (MQ=255)
                                                                                                             aCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTc                                 >  2:865713/1‑150 (MQ=255)
                                                                                                               gggTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGc                              <  2:108030/151‑1 (MQ=255)
                                                                                                               gggTATGGCGGTGGGGTTGTTTACCCGTTTCAACCAGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCa                                <  1:757449/149‑1 (MQ=255)
                                                                                                                        ggTGGGGTTGGTTATCCGTTTCACCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGc                       <  2:903217/149‑1 (MQ=255)
                                                                                                                        ggTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGcc                      >  1:199368/1‑150 (MQ=255)
                                                                                                                           ggggTTGGTTTTCCTTTTCACCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGa                   <  1:586904/150‑1 (MQ=255)
                                                                                                                            gggTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCACCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCg                     <  2:456915/147‑1 (MQ=255)
                                                                                                                                tggttATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc             >  1:701784/1‑151 (MQ=255)
                                                                                                                                tggttATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc             >  2:674340/1‑151 (MQ=255)
                                                                                                                                tggttATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc             >  1:427511/1‑151 (MQ=255)
                                                                                                                                 gtttATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc             <  1:867071/148‑1 (MQ=255)
                                                                                                                                   ttATCCGTTCCACCCAGGGTCGGGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGa           <  2:515150/150‑1 (MQ=255)
                                                                                                                                   ttACCCGTTTCAACCAGGGTCATCCCGGTCCAGCCCCCGCCTGTGAACCGCTGATCGTCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGTCCGGAACATc              <  2:449133/147‑1 (MQ=255)
                                                                                                                                     aTCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATc         >  1:645908/1‑150 (MQ=255)
                                                                                                                                     aTCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATc         >  1:226243/1‑150 (MQ=255)
                                                                                                                                       ccGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCAt       <  1:182207/150‑1 (MQ=255)
                                                                                                                                         gTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGAc    >  2:700167/1‑151 (MQ=255)
                                                                                                                                            tCAACCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCg  <  2:137669/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCG  >  NC_000913/2509716‑2510005

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