Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 594,198 A→G 11.7% intergenic (‑56/+54) oprG ← / ← PP_0505 outer membrane protein OprG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947594,1980AG11.7% 86.5 / 3.4 34intergenic (‑56/+54)oprG/PP_0505outer membrane protein OprG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/18);  new base G (2/2);  total (14/20)
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.05e-01

CGCAGGATCATGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCTTTGACCTGCATCTTTCCCTC  >  NC_002947/594051‑594347
                                                                                                                                                   |                                                                                                                                                     
cgcAGGATCATGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAAc                                                                                                                                                     <  2:596791/150‑1 (MQ=255)
     gATCATGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGaaa                                                                                                                                                >  2:58088/1‑150 (MQ=255)
        cATGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGtt                                                                                                                                             <  1:99801/150‑1 (MQ=255)
         aTGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGttt                                                                                                                                            <  2:40398/150‑1 (MQ=255)
            tCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGaaa                                                                                                                                                >  1:508567/1‑143 (MQ=255)
            tCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGaaa                                                                                                                                                <  2:508567/143‑1 (MQ=255)
               cccGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCa                                                                                                                                      <  2:687265/150‑1 (MQ=255)
                ccGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAATGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAg                                                                                                                                     >  1:96621/1‑150 (MQ=255)
                           gCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGc                                                                                                                          >  2:465245/1‑150 (MQ=255)
                             gTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGt                                                                                                                        <  1:465245/150‑1 (MQ=255)
                              tGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAAGGAAGCGTc                                                                                                                       >  1:516702/1‑150 (MQ=255)
                                                  cgcAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGtt                                                                                                   <  2:296356/150‑1 (MQ=255)
                                                     aaGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTtctc                                                                                                <  1:58088/150‑1 (MQ=255)
                                                     aaGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTtctc                                                                                                >  1:642183/1‑150 (MQ=255)
                                                                               aGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTtc                                                                      <  2:642183/150‑1 (MQ=255)
                                                                                   cTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCt                                                                                                                      >  2:434434/1‑98 (MQ=255)
                                                                                   cTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCt                                                                                                                      <  1:434434/98‑1 (MQ=255)
                                                                                         cATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAGACGGAAAGtt                                                                                                                                             <  2:570915/69‑1 (MQ=255)
                                                                                         cATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAGACGGAAAGtt                                                                                                                                             <  2:571278/69‑1 (MQ=255)
                                                                                         cATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAGACGGAAAGtt                                                                                                                                             >  1:571278/1‑69 (MQ=255)
                                                                                         cATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAGACGGAAAGtt                                                                                                                                             >  1:570915/1‑69 (MQ=255)
                                                                                            ggAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGa                                                                                                                                          >  2:101243/1‑69 (MQ=255)
                                                                                            ggAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGa                                                                                                                                          <  1:101243/69‑1 (MQ=255)
                                                                                                  tCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGAc                                                                                                                                         >  2:443258/1‑64 (MQ=255)
                                                                                                  tCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGAc                                                                                                                                         <  1:443258/64‑1 (MQ=255)
                                                                                                    ccAAATGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCg                                                 <  2:96621/150‑1 (MQ=255)
                                                                                                     cAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCgt                                                <  2:516702/150‑2 (MQ=255)
                                                                                                            cGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGc                                         <  2:342509/150‑1 (MQ=255)
                                                                                                                   aaTCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCg                                  >  2:376989/1‑150 (MQ=255)
                                                                                                                        agtagtCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCgctgct                             >  2:525691/1‑150 (MQ=255)
                                                                                                                             tCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAg                        >  2:6248/1‑150 (MQ=255)
                                                                                                                                 aTGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCttt                    <  1:6248/150‑1 (MQ=255)
                                                                                                                                    gTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCTTTGAc                 <  1:376989/150‑1 (MQ=255)
                                                                                                                                              ggTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCTTTGACCTGCATCttt       <  1:680057/150‑1 (MQ=255)
                                                                                                                                                   aaCGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCTTTGACCTGCATCTTTCCctc  >  1:308697/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCAGGATCATGTCCCCCGCCTGATAGGCGTGGGCGGCAGGGGCGGCGAGCGCAAGGGCCACAAGCGAGGCACCCAGCAAGGACTTGTTCATGGAAGCTCCCAAAGGACGTTATTAATCGAGTAGTCCAATGGTACGGAGCTGGTTAAACGGAAAGTTTGACCCAGCTCAACGAAGCGTCTTCACGAATTCGAAACAGTTCTCATCTCAGCTCATAAGCGTAGATTTTCTCGGCTTCCATCTGGTAGCCGGCGTCGGCCAGCTCGCTGCTGGAAGCTTTGACCTGCATCTTTCCCTC  >  NC_002947/594051‑594347

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: