Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,009,200 T→C 19.5% intergenic (+41/+41) norM → / ← PP_5263 multidrug resistance protein NorM/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,009,2000TC19.5% 38.0 / 2.4 20intergenic (+41/+41)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (9/7);  new base C (2/2);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGG  >  NC_002947/6009070‑6009340
                                                                                                                                   |                                                                                                                                            
gcggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGTCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGc                                                                                                                           <  1:717871/150‑1 (MQ=255)
        gcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGt                                                                                                                     >  2:658584/1‑148 (MQ=25)
        gcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGt                                                                                                                     <  1:658584/148‑1 (MQ=25)
                ttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                 >  2:93633/1‑128 (MQ=255)
                ttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                 <  1:93633/128‑1 (MQ=255)
                             gCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                     <  1:97329/127‑1 (MQ=255)
                             gCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                     >  2:97329/1‑127 (MQ=255)
                                     aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGt                                                                                                                     >  1:555642/1‑119 (MQ=18)
                                     aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGt                                                                                                                     <  2:555642/119‑1 (MQ=18)
                                                                      aaGCTGCTGTTTCGGCCTGACTTCCGTGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcg                                                                                                                                           <  1:609966/64‑1 (MQ=255)
                                                                      aaGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcg                                                                                                                                           >  2:609966/1‑64 (MQ=255)
                                                                                    gCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCggg                                                                                                                        <  1:619179/69‑1 (MQ=255)
                                                                                    gCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCggg                                                                                                                        >  2:619179/1‑69 (MQ=255)
                                                                                                             gCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAg                                                                                                <  2:217616/68‑1 (MQ=255)
                                                                                                             gCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAg                                                                                                >  1:217616/1‑68 (MQ=255)
                                                                                                                  cccGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTg         >  2:529443/1‑150 (MQ=255)
                                                                                                                       tCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCt    >  1:138212/1‑150 (MQ=255)
                                                                                                                         ccACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTgg  >  1:440277/1‑150 (MQ=255)
                                                                                                                           acaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGcc            <  2:232428/138‑1 (MQ=255)
                                                                                                                           acaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGcc            >  1:232428/1‑138 (MQ=255)
                                                                                                                                   |                                                                                                                                            
GCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGG  >  NC_002947/6009070‑6009340

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: