Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,818 G→C G214A (GGC→GCC)  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,8180GC93.3% 79.7 / ‑2.3 30G214A (GGC→GCC) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (13/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

TCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAA  >  NC_000913/4638569‑4639051
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
acAACGCGCCAAAGATCACCTTAGGCAAGCCACTTAGTGTGATTAGCGTCGGGGATGCGCGTGGGAGTGGTGTGGTGAAAGCGCCGGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTggc                                                                                                                                                                                                                                           <  2:210066/250‑1 (MQ=255)
              atCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTACCTGAAGGGGGCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGGGTAATAGCgg                                                                                                                                                                                                                             >  1:142481/1‑250 (MQ=255)
                catcTTAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAGGGGGACTGGCGGGAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTAAGCGGAACCGGCGATCATTCTTGGGTGCCTGCCGGGCGGAATAGctac                                                                                                                                                                                                                            >  1:214374/4‑247 (MQ=255)
                  tcTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGGAATg                                                                                                                                                                                                                         >  1:333462/2‑251 (MQ=255)
                     tAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAAt                                                                                                                                                                                                                          <  1:285870/247‑1 (MQ=255)
                     tAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGa                                                                                                                                                                                                                        <  1:189292/249‑1 (MQ=255)
                      aGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAACGGGACCGGCGATATTTCTGTGGTGCCTTCCGGGCGTAATAGCGAAATTACCt                                                                                                                                                                                                                     >  1:283933/1‑251 (MQ=255)
                            gccctttATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGGAGGGGCCTGGCGGAAGCGGGACCTGGAGATGAATTAGGCCCCGGAGTTATACGGGACCGGCGGTCTTTTTGTGGTGCCTGCCGTGGATAATAGCGagagacacttaacc                                                                                                                                                                                                                >  2:59510/5‑237 (MQ=255)
                                      gtgATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGACTATGGCCCTGAATTTAAGCGGTACCGGCGATATTTTTGTGGTGCCTGCCGGGCGTAATAGCGGAATGACCTTAACCAGCAACTGGcc                                                                                                                                                                                                     >  1:241807/1‑251 (MQ=255)
                                              gggcGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTTAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCa                                                                                                                                                                                             <  2:333462/248‑1 (MQ=255)
                                              cGTCGGGGTTGAGCGTGGTATTGGTGTTGTGAAAGCGCCGTAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAATCAGCACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTTCCGGCGATCTTTCTGTGGTGCCTGCCGGTCGTAATAGCGAAATGACCTTAACCAGCAACAGGCCGCATCCCa                                                                                                                                                                                             <  2:291218/251‑1 (MQ=255)
                                               gTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGGAATGACCTTAACCAGCAACTGGCCGCCTccc                                                                                                                                                                                              >  1:118128/1‑249 (MQ=255)
                                                gcGGGGATGCGCGGGGTATTGGTGTGGTGAAAGCGCCTGAATGTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAATCATAACCTGAAGCTGAATATGGCCCTGAATTGAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGt                                                                                                                                                                                           <  2:230212/250‑1 (MQ=255)
                                                      aTGCGAGTGGTATGGGTGTTGTGAAAGCGCATGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCTTTCCCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTTAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCAGCCCAGttttgta                                                                                                                                                                                     <  2:193414/251‑4 (MQ=255)
                                                                 aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAtttttt                                                                                                                                                                          >  1:235578/1‑249 (MQ=255)
                                                                 aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                                                                                           <  1:22136/250‑2 (MQ=255)
                                                                                 cgcCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGaa                                                                                                                                                         <  1:351465/251‑1 (MQ=255)
                                                                                      gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAAACAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACggggaggt                                                                                                                                                      >  1:243596/1‑244 (MQ=255)
                                                                                           tAACGGAACGGCGCTGCCCATTGACCCCTTCCCCGGGTTAGAGCACGGCGGGCAGAGCGTGCATATCCCTTTACCCGAAGGGGACTGGCTGAAGCAGAACCTGAATGTGAATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGTCGTAATAGCGCAATGACCTTAACCAGCAACTGGCCGAATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAACCGGTAAGTTAGCgcg                                                                                                                                                <  2:284779/250‑3 (MQ=255)
                                                                                               ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGGGGTACCGGCGATCTTTCTGTGGATCCTGCCGGGCGTAATAGCGAAATTGCCTTAACCAGCAACTGGCCGCATCCCAGGTTTTTAGGTGATTTTTTGCCCGCCAAACGGGAAATTaacggggaggg                                                                                                                                           >  2:309966/1‑241 (MQ=255)
                                                                                                     gcgcTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCAGAACCGGAATCTGAATATTGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTTATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCa                                                                                                                                     <  2:243596/251‑1 (MQ=255)
                                                                                                          gACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGgcgc                                                                                                                                <  1:93614/251‑1 (MQ=255)
                                                                                                           cccATTGAACCCGGCACCGGGTTAGAGCAAGACGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAATCTGAATATGGCCCTTAAGTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCa                                                                                                                               <  2:328199/250‑1 (MQ=255)
                                                                                                                  aaCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCa                                                                                                                         >  1:109061/1‑250 (MQ=255)
                                                                                                                            cGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTGTCCCGAAGGTGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCATGCGCAATGGCAAAGCAGCTGGt                                                                                                              <  2:258279/251‑1 (MQ=255)
                                                                                                                                  agagCAAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCta                                                                                                        <  1:299413/251‑1 (MQ=255)
                                                                                                                                      cAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaataa                                                                                                    <  1:59510/251‑1 (MQ=255)
                                                                                                                                                                  tACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGAACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTTGGTGATTTTTTGCACACCAAACGGGAAGGTAGCGCGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTTATAATAATCTCGGTGAGCGTTTTGCTTCAgg                                                                            >  2:126504/1‑247 (MQ=255)
                                                                                                                                                                                            aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTccccgc                                              >  1:62440/1‑248 (MQ=255)
                                                                                                                                                                                                                                    gCGATATTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCa       >  1:128722/1‑250 (MQ=255)
                                                                                                                                                                                                                                           ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa  <  1:126504/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
TCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAA  >  NC_000913/4638569‑4639051

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: