Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,830 Δ1 bp 11.1% coding (264/597 nt) PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,8300G.11.1% 97.1 / 8.7 36coding (264/597 nt)PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base G (18/13);  new base . (2/2);  total (20/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 = 6.97e-01

CAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489681‑5489972
                                                                                                                                                      |                                                                                                                                              
caaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATg                                                                                                                                                <  2:530891/150‑1 (MQ=255)
 aaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATgg                                                                                                                                               <  1:197296/150‑1 (MQ=255)
       ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGATGTACCTGTTCGTGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  1:485018/150‑1 (MQ=255)
         tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCCGTCCGGGTGATCGCTCTTCGTAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTTCTGGCTGCGGGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGc                                                                                                                                       >  2:4961/1‑150 (MQ=255)
              ggTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGccgaaa                                                                                                                                  <  1:511725/150‑6 (MQ=255)
                 cgcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggc                                                                                                                               <  2:533797/150‑1 (MQ=255)
                          aGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGc                                                                                                                      >  2:56154/1‑150 (MQ=255)
                           gCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCa                                                                                                                     <  1:751603/150‑1 (MQ=255)
                               tGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACa                                                                                                                 >  1:416608/1‑150 (MQ=255)
                                                 tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                               >  1:146634/1‑150 (MQ=255)
                                                    aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGACATgg                                                                                                                                               <  2:519373/99‑1 (MQ=255)
                                                    aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGACATgg                                                                                                                                               >  1:519373/1‑99 (MQ=255)
                                                     ccTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCAt                                                                                           <  1:489405/150‑1 (MQ=255)
                                                      cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACGGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc                                                                                          >  2:291161/1‑150 (MQ=255)
                                                       tGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCg                                                                                         >  2:659760/1‑150 (MQ=255)
                                                                        ttCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCCTCGGTGCTGGCATCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGACGCGACCATCGTACCCCTGACCATTCCg                                                                        >  2:172413/1‑150 (MQ=255)
                                                                           ggAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTg                                                                     >  2:674128/1‑150 (MQ=255)
                                                                            gAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGTAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGa                                                                    >  2:570142/1‑150 (MQ=255)
                                                                                  tCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTg                                                              <  1:321238/150‑1 (MQ=255)
                                                                                    ggCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGTCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGc                                                            >  2:499939/1‑150 (MQ=255)
                                                                                     gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                           >  2:836022/1‑150 (MQ=255)
                                                                                           ccgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCa                                                     >  1:31801/1‑150 (MQ=255)
                                                                                                      ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGcc                                                                                                                                      >  2:842393/1‑57 (MQ=25)
                                                                                                      ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGcc                                                                                                                                      <  1:842393/58‑2 (MQ=25)
                                                                                                      ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                          >  2:3165/1‑150 (MQ=255)
                                                                                                           cATCGGTGCTGTCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGTACCATCGGTGCCctgct                                     <  2:81364/150‑1 (MQ=255)
                                                                                                           cATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCGTGCGCATGGCCCAGGGCAAGCCGGCGGTACAGGTAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgct                                     <  1:4961/150‑1 (MQ=255)
                                                                                                           cATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgct                                     <  1:291161/150‑1 (MQ=255)
                                                                                                              cGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGt                                  <  1:674128/150‑1 (MQ=255)
                                                                                                                       cAGCGTGCTGTTCATTTCGCCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                         <  2:561625/150‑1 (MQ=255)
                                                                                                                         tcGTGCTCTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                       <  1:32824/149‑1 (MQ=255)
                                                                                                                                     tttCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                           <  1:811696/134‑1 (MQ=255)
                                                                                                                                     tttCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                           >  2:811696/1‑134 (MQ=255)
                                                                                                                                          gCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt      >  2:781990/1‑150 (MQ=255)
                                                                                                                                           cccTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                 >  1:760633/1‑106 (MQ=255)
                                                                                                                                           cccTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                 <  2:760633/106‑1 (MQ=255)
                                                                                                                                           cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGATAACGTTCAGCAGAATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATg     >  1:322719/1‑150 (MQ=255)
                                                                                                                                              tGGGGATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaa  >  2:250496/1‑150 (MQ=255)
                                                                                                                                                    aTGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                           >  2:202158/1‑88 (MQ=255)
                                                                                                                                                    aTGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                           <  1:202158/88‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489681‑5489972

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