Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,826:1 +C 18.2% coding (260/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,8261.C18.2% 90.1 / 21.5 33coding (260/597 nt)PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base . (10/17);  new base C (4/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATG  >  NC_002947/5489683‑5489969
                                                                                                                                                |                                                                                                                                               
aCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTTCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGc                                                                                                                                           <  1:210270/150‑1 (MQ=255)
  tGGCCTGTAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGccatgccc                                                                                                                                         <  1:331986/150‑8 (MQ=255)
     ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                      <  1:139832/150‑1 (MQ=255)
      cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                     <  1:43512/150‑1 (MQ=255)
         gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCaa                                                                                                                                  <  2:143511/150‑1 (MQ=255)
                gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                           <  1:374953/150‑1 (MQ=255)
                gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                           <  1:487892/150‑1 (MQ=255)
                 cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcgg                                                                                                                          <  2:382853/150‑1 (MQ=255)
                  gTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATAGATCCGACATCGGTGCTGGCGGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGt                                                                                                                         >  1:261531/1‑150 (MQ=255)
                     cGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACa                                                                                                                      <  2:225194/150‑1 (MQ=255)
                          ccTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTACTTCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGTCAGGGCCCAGGGCAAGCCGGCGGTACAGGCAg                                                                                                                 >  1:348826/1‑150 (MQ=255)
                                  gCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCCTCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCCGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                         <  2:467476/150‑1 (MQ=255)
                                               tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGACGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCCGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                            <  1:539494/150‑1 (MQ=255)
                                                  aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                <  2:530842/143‑1 (MQ=255)
                                                  aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                >  1:530842/1‑143 (MQ=255)
                                                    cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACg                                                                                                           >  2:580531/1‑130 (MQ=255)
                                                    cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACg                                                                                                           <  1:580531/130‑1 (MQ=255)
                                                                  gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAt                                                                         >  2:278951/1‑150 (MQ=255)
                                                                                   gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGcc                                                                                                                                   <  1:438556/75‑2 (MQ=38)
                                                                                   gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGcc                                                                                                                                   >  2:438556/1‑74 (MQ=38)
                                                                                   gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc                                                                                       <  1:395150/119‑1 (MQ=255)
                                                                                   gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc                                                                                       >  2:395150/1‑119 (MQ=255)
                                                                                                  ccTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTg                                         <  2:261531/150‑1 (MQ=255)
                                                                                                       cGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg                                    >  1:312318/1‑150 (MQ=255)
                                                                                                            gggTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGCCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGt                               <  2:348826/149‑1 (MQ=255)
                                                                                                              gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                        >  2:284023/1‑75 (MQ=255)
                                                                                                              gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                        <  1:284023/75‑1 (MQ=255)
                                                                                                                          tGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                   <  2:347804/100‑1 (MQ=255)
                                                                                                                          tGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                   >  1:347804/1‑100 (MQ=255)
                                                                                                                               ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt                                                                                              >  1:578794/1‑68 (MQ=255)
                                                                                                                               ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt                                                                                              <  2:578794/68‑1 (MQ=255)
                                                                                                                                      ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg     >  2:573526/3‑150 (MQ=255)
                                                                                                                                        gCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt   >  1:548357/1‑150 (MQ=255)
                                                                                                                                         cccTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATg  >  2:369824/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATG  >  NC_002947/5489683‑5489969

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