Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,207 C→G 9.1% intergenic (+30/+38) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,2070CG 9.1% 128.2 / 3.3 44intergenic (+30/+38)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base G (2/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAG  >  NC_002947/5790061‑5790355
                                                                                                                                                  |                                                                                                                                                    
ggggTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtg                                                                                                                                                   >  2:247446/1‑150 (MQ=255)
 gggTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtgt                                                                                                                                                  <  2:575071/150‑1 (MQ=255)
   gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGGGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                                <  1:412855/150‑1 (MQ=255)
       tCTGGCACCGGTGGGGCTGGAACGGTTTCGTGCTGTTCATTGGCTTCTTGCTGGCGGTGGCCTTGATGGTGGCGGTGCACCTTAGCAAGTTGCCACCCAAAACTGTGTGATGCATTCGTGGGCACGCCCGCATCCACAACCTGTGGGAGc                                                                                                                                            <  2:177663/150‑1 (MQ=255)
         tGGCACCAGTGGGGCGGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                          <  2:219657/150‑1 (MQ=255)
                                gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatca                                                                                                                   >  1:702613/1‑150 (MQ=255)
                                   cGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagat                                                                                                                <  2:290821/150‑1 (MQ=255)
                                    gggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagctg                                                                                                               <  1:247446/150‑4 (MQ=255)
                                       cTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGt                                                                                                            <  1:38383/150‑1 (MQ=255)
                                          ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                                         >  2:315171/1‑150 (MQ=255)
                                          ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                                         <  2:3719/150‑1 (MQ=255)
                                           tCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgat                                                                                                        <  2:702613/150‑1 (MQ=255)
                                             aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgatga                                                                                                      >  2:471392/1‑150 (MQ=255)
                                                ggCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCTGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                   <  2:283295/150‑1 (MQ=255)
                                                 gCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtctc                                                                                                  <  2:364214/150‑1 (MQ=255)
                                                   aGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCa                                                                                                                      <  2:522471/128‑1 (MQ=255)
                                                   aGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCa                                                                                                                      >  1:522471/1‑128 (MQ=255)
                                                            gcgcTGGCCTTGATGGTGGTGCGCCACCTGAGCAATTTGCCACCCAAAACGGTGTGTTGCATTCGTGGTCACGCCCTCAATCTCAACCTGTGGGAGCAGATGTCACCACGAAGCGGCCATCAGATCAGTTGATGATCTCGACAATatcca                                                                                       <  2:327634/150‑1 (MQ=255)
                                                            gCGGTGGCCTTGATGGTGGCGTTGCCCCTCAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGCGCACGCCCGCACCCACAACCTGTGAGAGCAGGTCCGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcca                                                                                       <  2:327635/150‑1 (MQ=255)
                                                                    cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCGCGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                               >  1:197753/1‑150 (MQ=255)
                                                                        aTGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGTCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTccc                                                                           >  2:824318/1‑150 (MQ=255)
                                                                             ggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtc                                                                      <  1:471392/150‑1 (MQ=255)
                                                                                     aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTg                                                              >  1:739366/1‑150 (MQ=255)
                                                                                               gTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAcc                                                    >  1:680479/1‑150 (MQ=255)
                                                                                                    cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                                    >  2:832912/1‑81 (MQ=255)
                                                                                                    cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                                    <  1:832912/81‑1 (MQ=255)
                                                                                                    cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCGCGGGTCATCAGGTGCTTTTCAACCTCAcc                                               >  1:483670/1‑150 (MQ=255)
                                                                                                     aCCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                   <  2:477527/129‑1 (MQ=255)
                                                                                                     aCCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                   >  1:477527/1‑129 (MQ=255)
                                                                                                     aCCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATACCAACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCg                                              >  1:264398/1‑150 (MQ=255)
                                                                                                       ccAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGt                                            >  2:507259/1‑150 (MQ=255)
                                                                                                                      gCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGt                             >  2:274072/1‑150 (MQ=255)
                                                                                                                            gTGGGCACGCCCGCACCCACAACCTGGGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTc                       <  1:507259/150‑1 (MQ=255)
                                                                                                                               ggCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGtt                    >  1:99305/1‑150 (MQ=255)
                                                                                                                                  tctcccGCTCCCACAGGTTGTGGGTGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATAACCACATCCACTTCCCGGCTCATCANGNGNTTTTCAACCTCACCNNTCAGCTTCACCTTGGTCTTGTCGTTGaa                 <  1:486149/147‑1 (MQ=255)
                                                                                                                                  acacccGCTCCCACAGGTTGTGGGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTg                   >  1:266588/4‑148 (MQ=38)
                                                                                                                                  acacccGCTCCCACAGGTTGTGGGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTg                   <  2:266588/145‑1 (MQ=38)
                                                                                                                                   cGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAc                <  1:824318/150‑1 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGTTGTGGGTGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTTGTTGAAcaga             >  2:486149/1‑147 (MQ=255)
                                                                                                                                        gCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCg                                              <  1:709421/115‑1 (MQ=255)
                                                                                                                                        gCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCg                                              >  2:709421/1‑115 (MQ=255)
                                                                                                                                            ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                   <  1:809310/90‑1 (MQ=255)
                                                                                                                                            ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                   >  2:809310/1‑90 (MQ=255)
                                                                                                                                                 aCCTGTGGGGTCGTGTGTGACCGCCAAGAGGCCATTAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTTCTTTTCAACCCCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAg  <  2:264398/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCAG  >  NC_002947/5790061‑5790355

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