Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,866,103 C→T 12.5% A259T (GCT→ACT)  PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,866,1030CT12.5% 76.9 / 3.2 32A259T (GCT→ACT) PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base T (2/2);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

CCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCC  >  NC_002947/3865956‑3866251
                                                                                                                                                   |                                                                                                                                                    
cccgccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGc                                                                                                                                                      >  6:33518/1‑148 (MQ=255)
                     gtgtGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACagaagaag                                                                                                                                <  1:245460/149‑1 (MQ=255)
                           gAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCg                                                                                                                          <  5:46017/149‑1 (MQ=255)
                                      cgCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGAGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCgg                                                                                                               >  5:188643/2‑149 (MQ=255)
                                       gCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGgc                                                                                                              <  8:22276/149‑1 (MQ=255)
                                          gTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGcc                                                                                                           >  1:246377/1‑149 (MQ=255)
                                               aCCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      >  5:157611/1‑149 (MQ=255)
                                                ccATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      >  7:159319/1‑148 (MQ=255)
                                                  atatGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCAt                                                                                                   <  5:186052/149‑1 (MQ=255)
                                                       cgcCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAgc                                                                                               <  1:243784/148‑1 (MQ=255)
                                                         ccGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAgctgg                                                                                            >  1:14193/1‑149 (MQ=255)
                                                         ccGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAgctgg                                                                                            <  5:33518/149‑1 (MQ=255)
                                                                 gggAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACAGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAg                                                                                    >  4:248924/1‑149 (MQ=255)
                                                                           tGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGcc                                                                           >  2:170183/1‑148 (MQ=255)
                                                                               aCCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAggg                                                                      <  2:246377/149‑1 (MQ=255)
                                                                                ccGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAggg                                                                      >  7:87655/1‑148 (MQ=255)
                                                                                 cGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGTCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      <  4:101738/115‑1 (MQ=255)
                                                                                 cGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGTCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      <  4:97772/115‑1 (MQ=255)
                                                                                 cGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGTCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      >  3:97772/1‑115 (MQ=255)
                                                                                 cGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGTCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCctgctg                                                                                                      >  3:101738/1‑115 (MQ=255)
                                                                                    ccATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATc                                                                 <  8:159319/149‑1 (MQ=255)
                                                                                     cATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACa                                                                                     <  5:224769/128‑1 (MQ=255)
                                                                                     cATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACa                                                                                     >  6:224769/1‑128 (MQ=255)
                                                                                            aCCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACag                                                                                                                                      <  4:149293/72‑1 (MQ=255)
                                                                                            aCCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACag                                                                                                                                      >  3:149293/1‑72 (MQ=255)
                                                                                            aCCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATgg                                                         <  8:87655/149‑1 (MQ=255)
                                                                                                   gTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCggg                                                  <  1:170183/149‑1 (MQ=255)
                                                                                                     tgCACCACCACCCCTTGACGGAAGTCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGGTACAGAAGAAGATACCGTCGACATCCGGCGCCTCCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGc                                                <  6:188643/148‑1 (MQ=255)
                                                                                                      gCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGc                                                <  4:188432/148‑1 (MQ=255)
                                                                                                                            ggCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCt                         >  6:63360/1‑149 (MQ=255)
                                                                                                                                    gcaCCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCa                 <  3:248924/149‑1 (MQ=255)
                                                                                                                                        cgcgCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGAt                                                           >  2:94968/1‑103 (MQ=255)
                                                                                                                                        cgcgCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGAt                                                           <  1:94968/103‑1 (MQ=255)
                                                                                                                                                   ccAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGcc  >  6:59317/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCC  >  NC_002947/3865956‑3866251

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