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,8590GA83.3% 99.1 / 0.8 42S77N (AGC→AAC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/19);  minor base C (0/7);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 3.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCG  >  NC_000913/2509711‑2510005
                                                                                                                                                    |                                                                                                                                                  
aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAAcc                                                                                                                                                  <  2:64014/151‑1 (MQ=255)
aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAAcc                                                                                                                                                  <  2:370643/151‑1 (MQ=255)
        gatgaAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCa                                                                                                                                           <  2:42538/150‑1 (MQ=255)
            aaaaTCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCg                                                                                                                                      >  2:82462/1‑151 (MQ=255)
                       ggTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGAcccc                                                                                                                           >  1:73163/1‑151 (MQ=255)
                        gTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCg                                                                                                                          >  2:208340/1‑151 (MQ=255)
                             aCAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCtg                                                                                                                      <  1:507857/150‑1 (MQ=255)
                             aCAAAATTTGCTCGGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATCCCCCCCTCTGGATCATCGGTGTATAAACGCTAACGGGTATTGCGGTGGGGTTGTTTATCCGTTTCAACCAGGGTCAGGCCGGACCAGCCCCCGCCtgt                                                                                                                     <  1:182579/151‑1 (MQ=255)
                               aGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCtgtg                                                                                                                    >  2:578148/1‑150 (MQ=255)
                                gAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGa                                                                                                                   >  2:579715/1‑150 (MQ=255)
                                   tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCg                                                                                                               >  2:207398/1‑151 (MQ=255)
                                       ctctGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTg                                                                                                            >  1:157678/1‑150 (MQ=255)
                                             cAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGgcg                                                                                                     >  1:40529/1‑151 (MQ=255)
                                             cAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATAAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCTTTTCACCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGgc                                                                                                      <  2:486730/150‑1 (MQ=255)
                                                      ccGGGAACTGTGGGGATAGCCCAGGATCCCCCCCTCTGGATCTTGGGTGTATAAACGCTAACGGGTATGGCGGTGGGGTTGTTTCTCCGTTTCAACCAGGGTCGGGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGtt                                                                                             <  1:161984/150‑1 (MQ=255)
                                                        gggACCTCTGGGGAGACCCCAGGATCCCCCCCTCTGAGTCATCGGTGAATTACCGCAAAGGGGTATTGCGGGGGGTTTGTTTACCCGTTCCACCCAGGGTCATGCGGAACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCg                                                                                          <  1:378398/151‑1 (MQ=255)
                                                        gggAACTCGGGGGAAACCCTGGGTTTCCCCCCTCGGGATCGTCGGTGTATAAAGGCTAGCGGGTATTGCGGGGGGTTTGGTTACCAATTTCACCCAGGTTCATGCCGTGCAAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTcc                                                                                           <  2:86564/150‑1 (MQ=255)
                                                              tctGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctct                                                                                    >  1:440011/1‑151 (MQ=255)
                                                                 ggggATAGCCCAGGATTCACCCCTCGGAATCATCGGTGTATTAACGCTAAGGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGTTCTTCCCGGACAAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctct                                                                                    <  2:55337/148‑1 (MQ=255)
                                                                 gggaATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAAGGGGTATTGCGGTGGGGTTGTTTACCCGTTTAAACAGGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctct                                                                                    <  2:479898/148‑1 (MQ=255)
                                                                  gggATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTgcgc                                                                                >  1:225273/1‑151 (MQ=255)
                                                                     aTAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTAc                                                                             <  2:169544/151‑1 (MQ=255)
                                                                         cccAGAATTCCCCCCTCGGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGTTTATCCGTTTCAACCAGGGTCAGGCCGCCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTgg                                                                         <  1:724829/151‑1 (MQ=255)
                                                                          ccAGGATCCCCCCCTCTGGAGCATCGGTGTATTAACGCTAACGGTTATTGCGGTGGGGTTGTTTATCCGTTTCAACCAGGGCCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCt                                                                           <  2:306566/148‑1 (MQ=255)
                                                                                ttCCCCCCTCTGGAACATCGGTGTATAAACCCTAACGGGTATGGCGGTGGGGTTGGTTATCCGTTTCACCCAGGGCCATGCCGAACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTa                                                                    <  1:509612/149‑1 (MQ=255)
                                                                                ttCACCCCTCTGGATCATCGGTTTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCGGGGCCAGGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTAt                                                                   <  1:654988/150‑1 (MQ=255)
                                                                                  cACCCCTCTGGATCATCGGTGTATTAACGTTACCGGGTATGGCGTGGGGGTTGGTTATCCGTTTCACCCAGGTTCATGCCGGCCCACACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATc                                                                  <  2:523938/149‑1 (MQ=255)
                                                                                               catcGGTGTATTAACGCTAACGGGTATTGGGGTGGGGTTGGTTATCCGTTTCAACCGGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATtct                                                    <  2:106908/150‑1 (MQ=255)
                                                                                                           aaCGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTgg                                        <  2:7139/150‑1 (MQ=255)
                                                                                                                 aaCGGGTATTGCGGTGGGGTTGTTTATCCGTTTCAACCGGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGgcggcg                                   <  2:775484/149‑1 (MQ=255)
                                                                                                                     ggTATTGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGcc                             >  2:314491/1‑151 (MQ=255)
                                                                                                                         ttGCGGTGGGGTTGGTTATCCTTTTCACCCAGGGTCAGGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTg                           <  1:496361/149‑1 (MQ=255)
                                                                                                                         ttGCGGTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg                          <  2:300727/150‑1 (MQ=255)
                                                                                                                             ggTGGGTTTGTTTACCCTTTCCACCAAGGTCCAGGCCGGACCAACCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTggggg                       <  1:566097/149‑2 (MQ=255)
                                                                                                                             ggTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCg                     >  1:265346/1‑151 (MQ=255)
                                                                                                                             ggTGGGGTTGGTTATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCg                     >  2:771574/1‑151 (MQ=255)
                                                                                                                             ggTGGGGTTGGTTACCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGcc                      <  2:608523/150‑1 (MQ=255)
                                                                                                                                ggggTTGTTTACCCTTTCCAACCAGGGTCATCCCGACCCACCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCgg                    <  1:496954/149‑1 (MQ=255)
                                                                                                                                      ggttACCCGTTTCAACCAGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATc              <  1:738870/149‑1 (MQ=255)
                                                                                                                                         tATCCGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATc         >  2:391013/1‑151 (MQ=255)
                                                                                                                                            ccGTTTCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATg      >  1:340070/1‑151 (MQ=255)
                                                                                                                                                ttCAACCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCg  >  1:318282/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCG  >  NC_000913/2509711‑2510005

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