Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,437,224 C→T 13.8% R114C (CGC→TGC)  PP_3932 → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,437,2240CT13.8% 67.7 / 3.9 29R114C (CGC→TGC) PP_3932diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base T (2/2);  total (15/14)
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.38e-01

CCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTT  >  NC_002947/4437085‑4437372
                                                                                                                                           |                                                                                                                                                    
ccTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAg                                                                                                                                            >  1:266863/1‑150 (MQ=255)
          ctggcCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGc                                                                                                                                  <  1:83288/150‑1 (MQ=255)
            ggcCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACgg                                                                                                                                   >  1:196227/1‑147 (MQ=255)
            ggcCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACgg                                                                                                                                   <  2:196227/147‑1 (MQ=255)
            ggcCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGt                                                                                                                                <  1:579731/150‑1 (MQ=255)
                   gTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGACTTACTCAGc                                                                                                                         <  1:303004/150‑1 (MQ=255)
                     cGGCGCCGACCTTCTGTGTGGCAGCGCCTGAGTACCGGCTCTGTGGTTAAACCTTTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCga                                                                                                                       <  2:299797/150‑2 (MQ=255)
                                     tgtgGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCCGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGcgcc                                                                                                       <  1:232121/150‑1 (MQ=255)
                                          cAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGcgccgccgc                                                                                                  >  2:601354/1‑150 (MQ=255)
                                                        cGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTggcgg                                                                                    >  2:512364/1‑150 (MQ=255)
                                                          gCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCgt                                                                                                                       <  1:319392/113‑1 (MQ=255)
                                                          gCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCgt                                                                                                                       >  2:319392/1‑113 (MQ=255)
                                                            tctGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAg                                                                                <  1:512364/150‑1 (MQ=255)
                                                             ctatgGTTAATCCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGc                                                                               >  2:460354/4‑150 (MQ=255)
                                                                  gTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGc                                                                          >  1:245278/1‑150 (MQ=255)
                                                                      aaCCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCg                                                                      <  2:266863/150‑1 (MQ=255)
                                                                       aCCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCgg                                                                     >  2:252791/1‑150 (MQ=255)
                                                                            tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCt                                                                                                 >  1:149346/1‑117 (MQ=255)
                                                                            tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCt                                                                                                 <  2:149346/117‑1 (MQ=255)
                                                                                              gtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGaa                                              <  1:487588/150‑1 (MQ=255)
                                                                                                   ccTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTc                                         >  2:440578/1‑150 (MQ=255)
                                                                                                                             cACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGt               >  2:572790/1‑150 (MQ=255)
                                                                                                                              aCGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTa                      <  1:611148/142‑1 (MQ=255)
                                                                                                                              aCGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTa                      >  2:611148/1‑142 (MQ=255)
                                                                                                                              aCGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGtt              >  2:225126/1‑150 (MQ=255)
                                                                                                                               cGCATGCGCCGGTGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTgg                         >  2:238418/1‑138 (MQ=255)
                                                                                                                               cGCATGCGCCGGTGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTgg                         <  1:238418/138‑1 (MQ=255)
                                                                                                                                gCATGCGCCGGTGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggtgg                                                                                       >  2:489361/1‑75 (MQ=255)
                                                                                                                                gCATGCGCCGGTGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggtgg                                                                                       <  1:489361/75‑1 (MQ=255)
                                                                                                                                          gcgcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGtt  >  2:274779/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTT  >  NC_002947/4437085‑4437372

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