Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,231,088 A→G 15.0% A24A (GCA→GCG)  yohK → LrgB family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,231,0880AG15.0% 88.3 / 11.0 47A24A (GCA→GCG) yohKLrgB family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (4/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

GTTTTTCTGGTGGTGAGCTGGAGTTCGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATCTGCTGCAACCGGCGGTCGTCGC  >  NC_000913/2230940‑2231231
                                                                                                                                                    |                                                                                                                                               
gtTTTTCTGGTGGTGAGCTGGAGTTCGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCa                                                                                                                                                 >  1:544322/1‑149 (MQ=255)
   tttCTGGTGGTGAGCTGGAGTTCGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACgg                                                                                                                                              <  2:427224/149‑1 (MQ=255)
                         cGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTTAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAAcc                                                                                                                        >  2:437372/1‑149 (MQ=255)
                         cGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAACCTGCCGGCACGGTATAAGTTTCCCATGCTTAAcc                                                                                                                        >  1:233776/1‑149 (MQ=255)
                         cGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAAcc                                                                                                                        >  1:133427/1‑149 (MQ=255)
                           cAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCg                                                                                                                      >  1:25267/1‑149 (MQ=255)
                           cAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCg                                                                                                                      >  1:79053/1‑149 (MQ=255)
                           cAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCg                                                                                                                      >  1:430741/1‑149 (MQ=255)
                                    cacGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTa                                                                                                             <  1:638633/149‑1 (MQ=255)
                                     acGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAg                                                                                                            >  1:442645/1‑149 (MQ=255)
                                             cGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGGACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAAtggtg                                                                                                    <  2:88795/149‑1 (MQ=255)
                                                aaaGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGTTGTTTGTAAAATGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGgtggtg                                                                                                 >  2:72330/1‑149 (MQ=255)
                                                  aGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGAt                                                                                               <  2:590288/149‑1 (MQ=255)
                                                               aaaGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTa                                                                                  >  2:524183/1‑149 (MQ=255)
                                                                 aGGATCAGAAGAATCATGGCGAATATCTGGTGGTCAGTACCGTTGACATTGATTGTCTTTTTTGCCGCCCGCAACCTGGCGGCACGGTATAAGGTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCTTTCCTTTTTTAAt                                                                                <  1:552280/149‑1 (MQ=255)
                                                                         aagaaTGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTg                                                                        >  1:308266/1‑149 (MQ=255)
                                                                         aagaaTGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTg                                                                        >  2:188135/1‑149 (MQ=255)
                                                                          agaaTGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTTTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTCACGCTGACCgg                                                                       <  2:696342/149‑1 (MQ=255)
                                                                               gatgGCGAATATCTGGTGGTCATTACCGCTGACTTTGATTTTCTTTTTTGCCGCCCGCAAACTGGAGGCACGGTATAAGTTTCCCTTGCTTAAACCGTTGCTGGTAGCACTGGTGGTGATCATTCGTTTTTTAATGGTGCCTGGctagt                                                                  >  2:363083/1‑145 (MQ=255)
                                                                               gatgGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTTCCTTTCGCACACTGGCGGCACGGTATAAGTTTCCCTTGCTTTACCCGTTGCGGGTAGCGATGGTGGTGACCATTCCTTTTTTGATACTGACTTGCAtag                                                                  >  2:79252/1‑147 (MQ=255)
                                                                               gatgGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGAGGTGATCATTCCTTTTTTAATGATGACTGGCAtct                                                                  >  2:315975/1‑149 (MQ=255)
                                                                                atgGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTACTTCTGACTGGCAtctc                                                                 >  1:370473/1‑149 (MQ=255)
                                                                                 tgGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCGTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTCGCATctct                                                                >  1:266172/1‑149 (MQ=255)
                                                                                                ggtCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACttt                                                 <  2:25267/149‑1 (MQ=255)
                                                                                                ggtCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACttt                                                 <  2:79053/149‑1 (MQ=255)
                                                                                                ggtCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCGCGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACttt                                                 <  2:430741/149‑1 (MQ=255)
                                                                                                 gtCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTCCTTTa                                                >  1:668032/1‑149 (MQ=255)
                                                                                                     ttACCGCTGACATTGATTGTCTTTTTTTTTTTTTTTTAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTTTTACTTCTGCTACTTTAAAgg                                            >  2:124783/1‑149 (MQ=255)
                                                                                                         cGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTACTGCTGACTGGCATCTCTTACGAACGCTACTTTAAAGGTAgg                                        >  2:343045/1‑148 (MQ=255)
                                                                                                            tGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGaa                                     <  2:442645/149‑1 (MQ=255)
                                                                                                            tGACATTGATTGTCTTTTTTGCCGCCCGCAAACAGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGaa                                     >  1:655174/1‑149 (MQ=255)
                                                                                                                attgattgTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACGGCTACTTTAAAGGTAGCGAAGTGc                                 <  1:751137/149‑1 (MQ=255)
                                                                                                                attgattgTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGc                                 <  1:554677/149‑1 (MQ=255)
                                                                                                                attgattgTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGc                                 <  2:686500/149‑1 (MQ=255)
                                                                                                                attgattgTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGc                                 <  1:554681/149‑1 (MQ=255)
                                                                                                                       tgcgtttttGCCGCCCCGAAACTGTCGGCACGGTATTAGTTTCCCTTGCTTCACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCGGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGa                          <  2:419581/145‑1 (MQ=255)
                                                                                                                       gTCTTTTTGGCCGCCCGCAAACTTGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGa                          <  2:509205/149‑1 (MQ=255)
                                                                                                                        tCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGCAGTGCTAAACGAt                         <  2:544322/149‑1 (MQ=255)
                                                                                                                        tCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGAt                         <  2:232348/149‑1 (MQ=255)
                                                                                                                        tCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGAt                         <  1:124783/149‑1 (MQ=255)
                                                                                                                        gcTTTTTTGCCGCTCGCAAACTGGCGGCACGGTATAAGTTTTCCTTGCTGAACCCGTTGCTGGTAGCAATGGTGGAGATCAGTCCTTTTTGAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGAt                         <  1:87051/148‑1 (MQ=255)
                                                                                                                          ttttttGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATct                       >  1:624719/1‑149 (MQ=255)
                                                                                                                              ttGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATctgctg                   <  2:67229/149‑1 (MQ=255)
                                                                                                                                     ccgcAACCTGGCGGCACGCTATAACTATCCCTAGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTATAATGCTGACTGGCAACTCTGACGACAGCTACTTTAAAAGTAGCGAAGTGCTAAACGATCTGCTGCAACcgg            <  1:701275/149‑1 (MQ=255)
                                                                                                                                     ccgcAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATCTGCTGCAACcgg            <  1:715214/149‑1 (MQ=255)
                                                                                                                                         aaaCTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATCTGCTGCAACCGGCGgt        <  1:72330/149‑1 (MQ=255)
                                                                                                                                           aCTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATAGTGGTGAGTATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATCTGCTGCAACCGGCGgtcg      >  1:415666/1‑149 (MQ=255)
                                                                                                                                               gcggcACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAACAGATCTGCTGCAACCAGCGGTCGTcgc  >  2:462998/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTTTTTCTGGTGGTGAGCTGGAGTTCGCAACTGGTACACGGTGAACGTAAAGTCGTAGGTCAGAAAGGATCAGAAGAATGATGGCGAATATCTGGTGGTCATTACCGCTGACATTGATTGTCTTTTTTGCCGCCCGCAAACTGGCGGCACGGTATAAGTTTCCCTTGCTTAACCCGTTGCTGGTAGCAATGGTGGTGATCATTCCTTTTTTAATGCTGACTGGCATCTCTTACGACAGCTACTTTAAAGGTAGCGAAGTGCTAAACGATCTGCTGCAACCGGCGGTCGTCGC  >  NC_000913/2230940‑2231231

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