Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,231 G→T 13.8% intergenic (+99/+106) PP_5531 → / ← PP_2505 hypothetical protein/GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,852,2310GT13.8% 73.6 / 3.4 29intergenic (+99/+106)PP_5531/PP_2505hypothetical protein/GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (7/18);  new base T (2/2);  total (9/20)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGG  >  NC_002947/2852082‑2852367
                                                                                                                                                     |                                                                                                                                        
cTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTg                                                                                                                                          >  1:268575/1‑150 (MQ=255)
   ggAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCGAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGtt                                                                                                                                       <  2:723944/150‑1 (MQ=255)
              ttACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAgggg                                                                                                                            >  1:252874/1‑150 (MQ=255)
                 cGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGggcgg                                                                                                                         <  1:16428/150‑1 (MQ=255)
                  gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCCCTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGggcggc                                                                                                                        <  2:252874/150‑1 (MQ=255)
                    gcagcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGGCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGa                                                                                                                      <  1:910408/150‑1 (MQ=255)
                      agcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTTATTGCAACCAGGggcgggaatg                                                                                                                    <  2:416299/150‑6 (MQ=255)
                         acaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGt                                                                                                                 <  2:431101/150‑1 (MQ=255)
                         acaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGGGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGCGCGGCGATCAGt                                                                                                                 <  2:431118/150‑1 (MQ=255)
                          caTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGtt                                                                                                                <  2:268575/150‑1 (MQ=255)
                              ttCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGaaa                                                                                                            <  1:57956/150‑1 (MQ=255)
                                cGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGc                                                                                                          <  2:971313/150‑1 (MQ=255)
                                  ccTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGa                                                                                                        <  1:297868/150‑1 (MQ=255)
                                           gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg                                                                                                                  <  1:868451/131‑1 (MQ=255)
                                           gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg                                                                                                                  >  2:868451/1‑131 (MQ=255)
                                           gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAg                                                                                               >  1:1030201/1‑150 (MQ=255)
                                                tGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCaaa                                                                                          >  2:535293/1‑150 (MQ=255)
                                                                       tAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGAt                                                                   >  2:477884/1‑150 (MQ=255)
                                                                                           gCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCgg                                               >  2:501932/1‑150 (MQ=255)
                                                                                                     gCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAAc                                     <  1:34610/150‑1 (MQ=255)
                                                                                                             aGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAtt                                                                                        >  2:1037146/1‑91 (MQ=255)
                                                                                                             aGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAtt                                                                                        <  1:1037146/91‑1 (MQ=255)
                                                                                                                  gcttggCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCCATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATaa                        <  1:68759/148‑1 (MQ=255)
                                                                                                                      gtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCTGGATCCGGAGATCTTAACTTCCTCAGGATAACGCt                    <  2:1043556/150‑1 (MQ=255)
                                                                                                                        ccAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTgg                  <  1:501932/150‑1 (MQ=255)
                                                                                                                            cGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGcaca              <  1:535293/150‑1 (MQ=255)
                                                                                                                                   gACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACaa       <  1:722183/150‑1 (MQ=255)
                                                                                                                                        gatcgccGCCCCTTATTGCAACCAGGGGCGGCGATAAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTgg  >  2:419924/6‑150 (MQ=255)
                                                                                                                                            gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATa            <  1:589509/135‑1 (MQ=255)
                                                                                                                                            gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATa            >  2:589509/2‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
CTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGG  >  NC_002947/2852082‑2852367

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: