Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,189 G→A 8.0% A713A (GCG→GCA)  ydbH → putative membrane‑anchored protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,445,1890GA 8.0% 147.6 / 3.0 50A713A (GCG→GCA) ydbHputative membrane‑anchored protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (2/2);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

ACAGCGAAAAATATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGAAAGATGGCTGGCTGGC  >  NC_000913/1445046‑1445332
                                                                                                                                               |                                                                                                                                               
aCAGCGAAAAATATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGtt                                                                                                                                            >  5:137021/1‑149 (MQ=255)
     gAAAAATATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGc                                                                                                                                       <  8:167913/149‑1 (MQ=255)
      aaaaaTATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCt                                                                                                                                      >  8:178657/1‑149 (MQ=255)
        aaaTATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTg                                                                                                                                     >  5:35970/1‑148 (MQ=255)
         aaTATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGaa                                                                                                                                   >  1:20443/1‑149 (MQ=255)
                cGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTACCGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAAtct                                                                                                                             <  7:120505/148‑1 (MQ=255)
                  gCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAAtctctc                                                                                                                          <  3:204070/149‑1 (MQ=255)
                         tGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAg                                                                                                                    <  3:93442/148‑1 (MQ=255)
                           cAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCg                                                                                                                  <  2:85163/148‑1 (MQ=255)
                             aGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAAc                                                                                                               <  2:176291/149‑1 (MQ=255)
                                   ttATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCGGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTa                                                                                                         <  6:137021/149‑1 (MQ=255)
                                         gTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCg                                                                                                   >  3:172464/1‑149 (MQ=255)
                                               tGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATc                                                                                             >  2:82042/1‑149 (MQ=255)
                                               tGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACAGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCACTCATCCAGCGAACTGGTTAGCGCCGTCAATc                                                                                             >  2:111417/1‑149 (MQ=255)
                                                 aagaagaaCCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATcc                                                                                            >  2:213738/1‑148 (MQ=255)
                                                 aagaagaaCCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCg                                                                                           <  8:73560/149‑1 (MQ=255)
                                                         cccTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCACTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAAtt                                                                                   >  5:181744/1‑149 (MQ=255)
                                                             tGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGACCTGGTTAGCGCCGTCAATCCGCAACAATTCGc                                                                                >  4:15770/1‑148 (MQ=255)
                                                              gCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGc                                                                                <  2:20443/147‑1 (MQ=255)
                                                                tGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGCCCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCTATCCGAAACACTTCGCCATg                                                                            >  3:141011/1‑149 (MQ=255)
                                                                       ggATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGCAACAATTACGTATGCCGACACATGCGCCGGCGCTGTTGCGGCGGAAGAATCGCTCATCCAGCGCACTGGGGCGCGCTGTCACTCCGCAACACTTCGCCAGGGCCgggg                                                                     >  2:201231/2‑149 (MQ=255)
                                                                         aTGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCgggg                                                                     <  4:172464/147‑1 (MQ=255)
                                                                            ttAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCAtt                                                                 <  7:178657/148‑1 (MQ=255)
                                                                                 gTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTCAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAgtgg                                                           >  1:211958/1‑149 (MQ=255)
                                                                                   cGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAgtggtg                                                         <  6:143618/149‑1 (MQ=255)
                                                                                                ggTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATg                                            <  6:35970/149‑1 (MQ=255)
                                                                                                ggTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATg                                            >  4:62089/1‑149 (MQ=255)
                                                                                                ggTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACACGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTAGCCATGTCCGGGGCAGTTAGTGGTGCACTGCCGTTATg                                            >  4:58945/1‑149 (MQ=255)
                                                                                                     cGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCgg                                                                       >  5:213540/1‑117 (MQ=255)
                                                                                                     cGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCgg                                                                       <  6:213540/117‑1 (MQ=255)
                                                                                                     cGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGa                                       <  4:114296/149‑1 (MQ=255)
                                                                                                     cGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTAATCCAGCGAACTGGTTAGCGCCGTCTATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTg                                        >  5:76014/1‑148 (MQ=255)
                                                                                                      gTACTGATGAAACAATTTCGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAAGTGGTTAGCGCCGTGAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGGCGTTATGGCTGaa                                      <  1:111417/149‑1 (MQ=255)
                                                                                                            atgaAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTAt                                             <  1:76497/136‑1 (MQ=255)
                                                                                                            atgaAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTAt                                             >  2:76497/1‑136 (MQ=255)
                                                                                                            atgaAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATg                                 >  2:173580/1‑148 (MQ=255)
                                                                                                            atgaAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGa                                >  5:46478/1‑149 (MQ=255)
                                                                                                            atgaAACAATTACGTATGCCGCAACATGACCCGGCGCTGATGCGGCTGAATAATCTCTCATACAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATg                                 >  4:45120/1‑148 (MQ=255)
                                                                                                               aaaCAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGaaa                              <  1:82042/148‑1 (MQ=255)
                                                                                                               aaaCAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGaaa                              >  7:144385/1‑148 (MQ=255)
                                                                                                                 aCAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAAt                           <  7:240498/149‑1 (MQ=255)
                                                                                                                   aaTTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATg                          >  8:16951/1‑148 (MQ=255)
                                                                                                                       aCGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTCAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGAtt                      <  2:211958/148‑1 (MQ=255)
                                                                                                                           aTGCCGCAACATGACCCGGCACTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGa                        >  3:214594/1‑142 (MQ=255)
                                                                                                                           aTGCCGCAACATGACCCGGCACTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGa                        <  4:214594/142‑1 (MQ=255)
                                                                                                                             gCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGa                  >  8:200995/1‑146 (MQ=255)
                                                                                                                              ccGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGa                  >  8:80580/1‑145 (MQ=255)
                                                                                                                              ccGCAACATGACCCGGCACTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGAAAGa              <  6:181744/149‑1 (MQ=255)
                                                                                                                                       gACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGAAAGAtggctggct     >  5:181144/1‑149 (MQ=255)
                                                                                                                                          ccGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTTCCGTTATGGCTGAACAATGAAAAATGGATTGTGAAAGAtggctggctggc  >  5:134137/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ACAGCGAAAAATATTACGGCTGATTTGCAAGGGCGTTATCCGTGGACTGAAGAAGAACCCTTGCTGTTGACTGATGTTAGCGTCGATGTGTTAGGCGGTAACGTACTGATGAAACAATTACGTATGCCGCAACATGACCCGGCGCTGTTGCGGCTGAATAATCTCTCATCCAGCGAACTGGTTAGCGCCGTCAATCCGAAACAATTCGCCATGTCCGGGGCATTTAGTGGTGCACTGCCGTTATGGCTGAACAATGAAAAATGGATTGTGAAAGATGGCTGGCTGGC  >  NC_000913/1445046‑1445332

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: ‑