Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,833,384 C→A 14.8% P2P (CCC→CCA vgrG‑II → type VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,833,3840CA14.8% 60.5 / 4.7 27P2P (CCC→CCAvgrG‑IItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base A (2/2);  total (13/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 = 3.52e-01

GGCCCTGGTGGACCCCGCTGGTCATCTCGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCATCAGCTCCCCGTTCCGACTACAACTGGAGCTGGTCAGTTTCGAAGACGATGTCGACTTCGCCCA  >  NC_002947/3833238‑3833533
                                                                                                                                                  |                                                                                                                                                     
ggCCCTGGTGGACCCCGCTGGTCATCTCGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGc                                                                                                                                                    <  2:405838/150‑1 (MQ=255)
                     tCATCTCGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCAGAAGTCGTCGAATGGCTCAACGCCATTCTGCAGCCGTAAAGGGGAACGCCACACGCACAGCCACTCCGCTCTGCGCtgcggg                                                                                                                               >  2:191400/1‑145 (MQ=255)
                      cATCTCGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACtt                                                                                                                              >  2:548906/1‑150 (MQ=255)
                        tctcGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTc                                                                                                                            <  1:30817/150‑1 (MQ=255)
                          tcGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCa                                                                                                                          >  2:249990/1‑150 (MQ=255)
                                 tCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTg                                                                                                                   >  1:490668/1‑150 (MQ=255)
                                     cATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTg                                                                                                               <  2:490668/150‑1 (MQ=255)
                                                      gctactacagctacCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGCTTggg                                                                                              >  1:101140/1‑148 (MQ=255)
                                                      gctaccacagctacCCCCTCCACCTGCTAACCGCCGAAGTCGTCAAATCTCTCAACCCCATTCTGCAGCAGTAACGGGGGATGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAg                                                                                              <  1:229357/150‑1 (MQ=255)
                                                           tacagctacCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAggtggt                                                                                         <  1:459554/150‑1 (MQ=255)
                                                           tacagctacCGCCTCGACCTGCTAACCGCCCAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGATTGGCGGTAGCGAGTTTGtggtggg                                                                                         >  1:470004/1‑149 (MQ=255)
                                                                              tGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGaa                                                                      <  1:249990/150‑1 (MQ=255)
                                                                                     cgccgAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCATCa                                                               >  1:34724/1‑150 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAggtggt                                                                                         <  2:30531/114‑1 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAggtggt                                                                                         >  1:30531/1‑114 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCAAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCa                                                                  >  2:203575/1‑137 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCAAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCa                                                                  >  2:203549/1‑137 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCAAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCa                                                                  <  1:203549/137‑1 (MQ=255)
                                                                                               gtcgAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCAAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCa                                                                  <  1:203575/137‑1 (MQ=255)
                                                                                                   aaTCGCGCAACGCAATTCTGAAGAAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGGGGTGGTGTCGTTCGAACTGGACGAAGCCATCAGCTCCCCGTTCCGa                                                 <  2:101140/150‑1 (MQ=255)
                                                                                                           aaCGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATcc                                                                                                                                              >  1:457867/1‑49 (MQ=255)
                                                                                                           aaCGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATcc                                                                                                                                              <  2:457867/49‑1 (MQ=255)
                                                                                                           aaCGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGa                                                                                                     <  1:407031/90‑1 (MQ=255)
                                                                                                           aaCGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGa                                                                                                     >  2:407031/1‑90 (MQ=255)
                                                                                                                            gTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTg                                                                            <  2:572465/98‑1 (MQ=255)
                                                                                                                            gTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTg                                                                            >  1:572465/1‑98 (MQ=255)
                                                                                                                             tAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCATCAGCTCCCCGTTCCGACTACAACTGGAGCTGGTCAGTTTCGa                       <  1:528452/150‑1 (MQ=255)
                                                                                                                                                  cAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCCGTTGGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCATCAGCTCCCCGTTCCGACTACAACTGGAGCTGGTCAGTTTCGAAGACGATGTCGACTTCGCCCa  >  1:204763/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCCCTGGTGGACCCCGCTGGTCATCTCGGTCATCGCCATGGTCTGCGCCTACAGCTACTACAGCTACCGCCTCGACCTGCTAACCGCCGAAGTCGTCGAATCGCTCAACGCCATTCTGCAGCAGTAACGGGGGACGCCACATGCCCAGCCAATCCGATCTGCGCTTCACTTTCCAGCCGCTGGTTGGCCGTAGCGAGTTTGAGGTGGTGTCGTTCGAACTGGACGAAGCCATCAGCTCCCCGTTCCGACTACAACTGGAGCTGGTCAGTTTCGAAGACGATGTCGACTTCGCCCA  >  NC_002947/3833238‑3833533

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: