Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,567 A→G E130E (GAA→GAG creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,5670AG90.3% 67.4 / 0.6 33E130E (GAA→GAGcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (17/11);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

AGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGT  >  NC_000913/4638325‑4638809
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
cggtgtccGCAAAAACTCGGTTGGCCGTTTTGCGCTATTCCTTTGTCTGAGCGTTATACGGTTCAGGAAGAGGATAAAACTGTGGTGCTGAAACGAAGTTTTATCCATTTTTGGTTACCTTAGTCAGTGATTGTTGATGGCATTCAGAACGTGGAAGAACGCAAGATAGCGACTTATGCCGGTCAGGTCTTTCTTAGTGCTTTAACGCTAAAAGTTGATTTAGATTTTTCGCGTCTTTGCGATCTCAAcgt                                                                                                                                                                                                                                              <  2:127580/245‑2 (MQ=255)
              aCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCatct                                                                                                                                                                                                                                >  1:193414/1‑248 (MQ=255)
                tCGTTGGGCAGCTTATCGCTATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACGTGAGTCATTGATGGTTGAGGGATATCAGAACGTGGATGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCTATTACGATGTTTCGCGTCTTAGCGGGCTCAACGCGCAAAATATCATCtt                                                                                                                                                                                                                               >  2:93614/1‑250 (MQ=255)
                    tGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGc                                                                                                                                                                                                                           >  1:328199/1‑250 (MQ=255)
                                                       ataCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGcgcg                                                                                                                                                                                        <  1:308408/250‑1 (MQ=255)
                                                         aCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGAGGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTg                                                                                                                                                                                      <  1:8672/250‑1 (MQ=255)
                                                                    agagGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATtggtgtgg                                                                                                                                                                          >  1:251194/1‑251 (MQ=255)
                                                                           aaaaCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTGTTGTGATTAGCGTCGGGGATGCGCGTGGTAGTGGTGTGGTGAAAg                                                                                                                                                                    <  1:135621/250‑1 (MQ=255)
                                                                                        ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaac                                                                                                                                                      >  1:58784/1‑251 (MQ=255)
                                                                                                        ccagttgtGGTTGCCTGAGTCATTGATTGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGGCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAcc                                                                                                                                       <  2:378130/243‑1 (MQ=255)
                                                                                                                 gTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCg                                                                                                                             <  1:100345/251‑1 (MQ=255)
                                                                                                                      cTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCAcc                                                                                                                        >  1:58432/1‑251 (MQ=255)
                                                                                                                       tGAGTCATTGATGGTTAATGGCCATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTAGCACAGGGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGTAATTAGCGTCGGGGATACGCGAGGTATTGGTGTGGTGAAAGCGCATGAAGTTAACAGAACGGGGCTGaaattttaagccggcgccg                                                                                                                       >  2:250878/1‑233 (MQ=255)
                                                                                                                                       gatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCgg                                                                                                       <  1:369332/251‑1 (MQ=255)
                                                                                                                                        atggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTAAGGTCTGGCACAGTGATTTACCGCTAAACGCCGATTTAGATGTTTCGAGTCTTAGCGAGCTCAACGCGCCACATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGGGGAGGAAGCGCCTGAAATTAACGGAACGGCGCCGACGATTTAACCCGGCACCGGGGTAGAGCAAAGCggg                                                                                                      >  2:384838/1‑251 (MQ=255)
                                                                                                                                             aaTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAggg                                                                                                 >  1:230212/1‑251 (MQ=255)
                                                                                                                                                    aCGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGAGCCAAATATCATCTTAGGCAAGACCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGAGATGACCATTGCACCCGACACCGGTTTAGAGCCAGGCGGGCAggggggggat                                                                                          >  2:259013/1‑244 (MQ=255)
                                                                                                                                                                 cAAGTTAGTGATTTATGCCGGTCAGGTCTGGTACAGTGATTTAACGCTAAAAACTGAGTTCGATGTTTCGCGTCTTAGCGAGCTCAACTCGCCAAATATCATCTTAGGCAAGCCCCTTATTGTGATTATCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCGGAAGTTAACGGAACTGCGCTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCttttcccg                                                                             <  2:102508/251‑6 (MQ=255)
                                                                                                                                                                       aGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAgggg                                                                      <  1:195437/251‑1 (MQ=255)
                                                                                                                                                                                            cTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTg                                                 >  1:258279/1‑250 (MQ=255)
                                                                                                                                                                                                       aTTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCCGAACCTGAAGCTGAATATg                                     >  1:210066/1‑251 (MQ=255)
                                                                                                                                                                                                           aaCGCTAAAAGCCGATTTCTATGTGTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCATGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAatattgcac                                 <  2:58784/251‑6 (MQ=255)
                                                                                                                                                                                                             cactAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGTAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGACAGGGGGTGCATATCCCTTTACCCGAAGGAGAACGGCGGAAACAGaaacagacaattgaaatgctccc                                 >  1:219066/4‑229 (MQ=255)
                                                                                                                                                                                                               ctAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGCAGTTAACGGAACGGCGCTGACCCTTGACCCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGAGAGCATATCCCTTTACCCGAAGGGGACCGGCGGGAGCAGGATCTGGAGCTggaaatgggcatgac                             >  2:22136/2‑237 (MQ=255)
                                                                                                                                                                                                               ctAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGGACCCAAAGCTGAATATGGCCCTGa                              >  1:74445/2‑250 (MQ=255)
                                                                                                                                                                                                                         aTTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAg                                              <  1:166537/224‑1 (MQ=255)
                                                                                                                                                                                                                           ttctaTGTTTCGCGTCTTAGCTATCCCATCGCGTCAAATATCAGCTTAGGCTTTTCCTTTTGTGTGATTAGCGTCGGGGATGCGCGTGGTATTGTTGTGGTGAAAGCGCCCGAAGTTAACGGAACGGCGCTTACCATTGAACTCGGCGCCGGGTTAGAGCAAGGCTGGCTGGGCGTTCATATCTCTTTATCCTAAGGGGACTGGCGTAAGCAGATCTTGAAGCGGAATATGGCTCTGAATTTaggtgata                  <  2:346976/246‑8 (MQ=255)
                                                                                                                                                                                                                           ttCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGGAGCTGGATAAGGCCCTGAATTTAAgggggac                 >  1:292915/1‑245 (MQ=255)
                                                                                                                                                                                                                            tCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCgggaac                >  1:318478/1‑247 (MQ=255)
                                                                                                                                                                                                                                 gTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCg            <  1:259013/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cttCTTAGCGAGCACAACGCGCCAAATATCATCTTAGGCAAGCTCTTTATTGTGATTAGCTTCGGGGATGCTCGAGGTATTGGCGTGTTGAAAGCGCCGGAAGTTAACGGATCGGCGCTGTCCATTGAACCCGGCTCCTGTTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCTTTCCCCGGAGTGGACTGGCGGAAGCAGTACCTTAAGCTGAATATGGCCCTGAATTTAAGCGGTACCTGTGATCtt       <  2:228434/247‑1 (MQ=255)
                                                                                                                                                                                                                                        gTCTCAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCTGAAGGGGACTGGCGGGAGCAGGACCATGAGCTGGATATGGCCCTGGATTTTAGCGGGACCAGCGGTCttttt    >  2:189292/1‑249 (MQ=255)
                                                                                                                                                                                                                                          cTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCCGAATTTTAGCGGTACCGGCGATCTTTCtgt  >  1:179235/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGT  >  NC_000913/4638325‑4638809

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