Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,370 C→T 28.6% V183I (GTA→ATA)  pctB ← methyl‑accepting chemotaxis protein PctB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,3700CT28.6% 17.8 / 5.7 14V183I (GTA→ATA) pctBmethyl‑accepting chemotaxis protein PctB
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base T (2/2);  total (8/6)
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

CATCACCTGGTCTTTATTGGGGCTGACGATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGCCCGGCGTTGGCGGCGGCGCCATACCAAGGGCGCTG  >  NC_002947/2568236‑2568505
                                                                                                                                       |                                                                                                                                        
catcaCCTGGTCTTTATTGGGGCTGACGATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTaca                                                                                                                                         >  1:128649/1‑136 (MQ=255)
                  ggggCTGACGATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTg                                                                                                                      >  2:557696/1‑136 (MQ=255)
                     gCTGACGATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTgcgc                                                                                                                   <  1:154395/136‑1 (MQ=255)
                           gATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGctc                                                                                                             >  2:505575/1‑136 (MQ=255)
                           gATCACCTGGCCATTGGCATCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGctc                                                                                                             >  2:106954/1‑136 (MQ=255)
                                   ggCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCATTATTATCCGATCAGTTCGCCATTGCGCTTGCTCTTCACCg                                                                                                      >  2:89405/1‑136 (MQ=255)
                                   ggCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCATTATTATCCGATCAGTTCGCCATTGCGCTTGCTCTTCACCg                                                                                                      >  2:317091/1‑136 (MQ=255)
                                                                    tGCCGACGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGc                                                                    <  1:444331/136‑1 (MQ=255)
                                                                                     gAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTaca                                                                                                                                         <  1:119729/51‑1 (MQ=255)
                                                                                     gAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTaca                                                                                                                                         >  2:119729/1‑51 (MQ=255)
                                                                                            tgatCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCAc                                            <  2:60679/136‑1 (MQ=255)
                                                                                                   gacgaGGGTGTCGAGGGTCACGTCACCGCCATTATTATCCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTg                                      <  1:89405/136‑1 (MQ=255)
                                                                                                   gacgaGGGTGTCGAGGGTCACGTCACCGCCATTATTATCCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTg                                      <  1:317091/136‑1 (MQ=255)
                                                                                                    acgaGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGcc                                    >  2:564503/1‑136 (MQ=255)
                                                                                                                          cacCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGCCCGGCGTTGGCGGCGGCGCCATa              >  2:96339/1‑136 (MQ=255)
                                                                                                                               ccgACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGc                                                                    >  2:480523/1‑77 (MQ=255)
                                                                                                                               ccgACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGc                                                                    <  1:480523/77‑1 (MQ=255)
                                                                                                                                       cacCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGCCCGGCGTTGGCGGCGGCGCCATACCAAGGGCGCTg  <  1:505575/136‑1 (MQ=255)
                                                                                                                                       cacCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGCCCGGCGTTGGCGGCGGCGCCATACCAAGGGCGCTg  <  1:106954/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CATCACCTGGTCTTTATTGGGGCTGACGATCACCTGGCCATTGGCGTCGGCCAGGAAGGCGTGACCGATGCCGCCGAAGTCGACCGAGTTGATGATCTCGACGAGGGTGTCGAGGGTCACGTCACCGCCGACTACACCGATCAGTTCGCCATTGCGCTTGCTCTTCACCGGCGTGGCGATGGTCACCATCAGGCCACCGACTGCACCCTGGTACGGGGCGGTCAGCACGGTCTGCCCGGCGTTGGCGGCGGCGCCATACCAAGGGCGCTG  >  NC_002947/2568236‑2568505

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: