Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,193 G→A 14.1% intergenic (+16/+52) 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,1930GA14.1% 65.6 / 2.7 28intergenic (+16/+52)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (2/2);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAC  >  NC_002947/5790060‑5790341
                                                                                                                                     |                                                                                                                                                    
cgggGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGAGCACGCCCGCACCcacaggtt                                                                                                                                      <  1:76149/150‑5 (MQ=255)
     tgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACACCCGCACCCACAACCTGTggg                                                                                                                                 >  1:135602/1‑150 (MQ=255)
       ttCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAg                                                                                                                               >  2:141812/1‑150 (MQ=255)
           ggCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCggg                                                                                                                           >  2:426839/1‑150 (MQ=255)
              aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCGGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                        >  2:284625/1‑150 (MQ=255)
              aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                        >  1:45958/1‑150 (MQ=255)
              aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                        >  1:160443/1‑150 (MQ=255)
                        cTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACACCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGaat                                                                                                              <  2:135602/150‑2 (MQ=255)
                            aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaatgcat                                                                                                          <  1:363155/150‑6 (MQ=25)
                                                     gCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGAc                                                                                 <  2:428897/150‑1 (MQ=255)
                                                          cTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAatat                                                                            <  2:134725/150‑1 (MQ=255)
                                                           tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGc                                                                                                                              >  1:167821/1‑99 (MQ=255)
                                                           tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGc                                                                                                                              <  2:167821/99‑1 (MQ=255)
                                                            ggCGGTGGCCTTTATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                          <  2:160443/150‑1 (MQ=255)
                                                            ggCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                          <  2:45958/150‑1 (MQ=255)
                                                              cGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccac                                                                        <  2:162631/150‑1 (MQ=255)
                                                                   gCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagat                                                                                                  >  1:445777/1‑119 (MQ=255)
                                                                   gCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagat                                                                                                  <  2:445777/119‑1 (MQ=255)
                                                                          gggtggCGTTGCCCCTGAGCAAGTTGCCCCCCACAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCg                                                            <  2:348508/149‑1 (MQ=255)
                                                                                 gTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                     >  1:22181/1‑150 (MQ=255)
                                                                                  ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatca                                                    >  2:492181/1‑150 (MQ=255)
                                                                                          gAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                            <  1:284625/150‑1 (MQ=255)
                                                                                              aaGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCaa                                        <  1:141812/150‑1 (MQ=255)
                                                                                                     cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAcc                                 >  1:340730/1‑150 (MQ=255)
                                                                                                                 tgtgATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAc                     <  2:340730/150‑1 (MQ=255)
                                                                                                                        tcttCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTAATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTc              >  1:166340/3‑150 (MQ=21)
                                                                                                                                   aCGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGaa   >  2:321374/1‑150 (MQ=255)
                                                                                                                                    cGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAc  <  1:321374/150‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
CGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAC  >  NC_002947/5790060‑5790341

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