Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,825,023 T→A 11.0% intergenic (‑86/+117) ves ← / ← spy cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0230TA11.0% 93.6 / 3.0 37intergenic (‑86/+117)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base T (12/20);  new base A (2/2);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACC  >  NC_000913/1824878‑1825166
                                                                                                                                                  |                                                                                                                                               
ttCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGAt                                                                                                                                              <  2:849340/149‑1 (MQ=255)
   aCCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCTTTCAGGTAACCGTGGATGGTTATCTGCAACTCGGGTGATTTACGTTAGTGGTGATCAg                                                                                                                                           <  1:559803/149‑1 (MQ=255)
          gcaTTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATTCTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCttttt                                                                                                                                    <  1:851996/149‑1 (MQ=255)
              ttCGCCACAGATTCACCGACATTTTACGCTTATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGAGGGTTATCGGCAACTCGGGTGGTGTACATTAGGGGTGATCAGGCTTTTTTTAc                                                                                                                                >  1:521326/1‑149 (MQ=255)
                 gCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGATCCCCAGCGAGTT‑TTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAAc                                                                                                                             <  1:623897/149‑1 (MQ=255)
                 gCCACAGATCCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCA‑CCAAGTACCCCGTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAAc                                                                                                                             <  1:199949/149‑1 (MQ=255)
                       gATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTg                                                                                                                       <  1:269395/149‑1 (MQ=255)
                         ttCACCGACATTTTACGCAGACCAAAGTATTCCATACTCCCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGGGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgag                                                                                                                     <  2:119575/149‑1 (MQ=255)
                          tCACCGACATTTTACGCATATCACAGTATTCCATCCTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgaga                                                                                                                    <  1:125644/149‑1 (MQ=255)
                          tCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGAGATCAGGCTTTTTTTACAACTGTTTgaga                                                                                                                    >  1:349169/1‑149 (MQ=255)
                          tCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGAGATCAGGCTTTTTTTACAACTGTTTgaga                                                                                                                    >  1:349224/1‑149 (MQ=255)
                                aCATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGCCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTCTACAACTCTTTgagaagagaa                                                                                                              <  1:606811/149‑1 (MQ=255)
                                   ttttACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGAGGGGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                           >  1:755999/1‑149 (MQ=255)
                                     gtACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                         <  1:250667/148‑1 (MQ=255)
                                             atCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAGCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGACAAAAGAAAAccgccg                                                                                                 >  2:578905/1‑149 (MQ=255)
                                              tCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                >  2:682782/1‑149 (MQ=255)
                                               cAAGGTATTCGAGGCTGGCT‑CCCGGCCCGTTCTTTTGATGTTTATCGGTCATGGAACCGTGGATGGGTATCGGCAACTCGTGTGATTTACGTTAGTGTTGATCAGGCTTTTTTTACAACTGTATGATAAGAGAAAAGAAAACCGCCGAt                                                                                               <  2:245399/149‑1 (MQ=255)
                                                          aTACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGACCCGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg                                                                                    <  1:455492/149‑1 (MQ=255)
                                                               cgcg‑CCGGGCCCGTTCTTTTCATTGTTACCTCTCAGGTAACCGTGGATTGTTATCGGCAACTTGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                               <  1:413410/149‑1 (MQ=255)
                                                               cGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                               <  2:565549/149‑1 (MQ=255)
                                                               cGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                               <  1:436507/149‑1 (MQ=255)
                                                               cGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGGGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                               <  2:23042/149‑1 (MQ=255)
                                                                      cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTACCCGTGGGTGGTTATCAGCAACTCGGGTGACTTACGTTAGTGGTGATCCGGCCTTTTTTACAACTGATTGAGATGATAAACCAAAACCGCCGCTTCCATCCACCGCACTACTgaac                                                                         >  2:622786/1‑146 (MQ=255)
                                                                           cGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCACGGAAg                                                                    >  1:384427/1‑149 (MQ=255)
                                                                               cTTTTCATTGTTATCGATCAGGTAACCGTGGCTGGTTATCGGCAACTCGGGTGATTTACGATAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAAAAGAGAAAAGAAAAGAGCCGATCCTGTCCACCGCATTACTGGAAGGTAGTGGa                                                                >  2:422851/1‑149 (MQ=255)
                                                                                 tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                              >  1:427460/1‑149 (MQ=255)
                                                                                 tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGAGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                              <  2:349169/149‑1 (MQ=255)
                                                                                 tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGAGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                              <  2:349224/149‑1 (MQ=255)
                                                                                              gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGCTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                 >  1:67448/1‑149 (MQ=255)
                                                                                              gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGCTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                 >  1:350392/1‑149 (MQ=255)
                                                                                               ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                >  1:837971/1‑149 (MQ=255)
                                                                                                    gTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGCGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGGTCCTGTCCACCGCATTACTGGAAGGTAGTGGACCAGACCGCCGGCCtcacgtg                                           >  2:760822/1‑143 (MQ=255)
                                                                                                                   tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGCGAGGAGGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                            <  1:847247/149‑1 (MQ=255)
                                                                                                                                ggggttttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg               <  2:755999/144‑1 (MQ=255)
                                                                                                                                       tACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt        <  1:136558/149‑1 (MQ=255)
                                                                                                                                           ttAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCGGGCAttttg    >  1:368488/1‑148 (MQ=255)
                                                                                                                                             aGTGGTGATCAGGCTTTTTTTTCACCTGTTTGAGAAGAGAAATGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGGCAAGACCGGCGGTCTTATGTTTTTTGGCTGAAAGATTATTCAGCAGTTCCAGGCATTTTAcc  <  2:524903/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCT‑CCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACC  >  NC_000913/1824878‑1825166

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