Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,280 C→T 12.9% P18P (CCG→CCA)  PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,2800CT12.9% 75.3 / 4.1 31P18P (CCG→CCA) PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (15/12);  new base T (2/2);  total (17/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 = 9.62e-01

CACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGG  >  NC_002947/4240133‑4240426
                                                                                                                                                   |                                                                                                                                                  
cACGAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCgg                                                                                                                                                  >  2:426449/1‑150 (MQ=255)
  ccAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCa                                                                                                                                                <  2:343380/150‑1 (MQ=255)
  ccAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCACTGCAACGCCACCTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCa                                                                                                                                                <  1:700080/150‑1 (MQ=255)
         gCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGc                                                                                                                                         <  2:759352/150‑1 (MQ=255)
            gATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTg                                                                                                                                      >  2:800301/1‑150 (MQ=255)
                                ccAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGcacac                                                                                                                  >  2:777254/1‑149 (MQ=255)
                                 cAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCATGCCCGCGCACAGc                                                                                                                 >  2:697563/1‑150 (MQ=255)
                                        tCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTgc                                                                                                          >  2:806402/1‑150 (MQ=255)
                                              gTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGACGCCTGGGTCATGCCCGACAGTCTAGCCTGCCAGGCACGCGCGCGCACAGCCCTTTGCGCCAgg                                                                                                    >  1:184852/1‑150 (MQ=255)
                                                     cATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCg                                                                                             >  1:635139/1‑150 (MQ=255)
                                                           ccAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTg                                                                                       <  1:806402/150‑1 (MQ=255)
                                                                          aCGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGcc                                                                                                                                    >  1:692107/1‑90 (MQ=255)
                                                                          aCGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGcc                                                                                                                                    >  1:692250/1‑90 (MQ=255)
                                                                          aCGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGcc                                                                                                                                    <  2:692250/90‑1 (MQ=255)
                                                                          aCGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGcc                                                                                                                                    <  2:692107/90‑1 (MQ=255)
                                                                             ccACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGt                                                                     <  1:426449/150‑1 (MQ=255)
                                                                               aCGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCATGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGa                                                                   <  1:697563/150‑1 (MQ=255)
                                                                                    ccACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGa                                                                   <  1:206110/145‑1 (MQ=255)
                                                                                    ccACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGa                                                                   >  2:206110/1‑145 (MQ=255)
                                                                                    ccACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCa                                                              >  1:547196/1‑150 (MQ=255)
                                                                                     cACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTTATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAc                                                             <  1:618138/150‑1 (MQ=255)
                                                                                      aCTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcg                                                            >  1:778088/1‑150 (MQ=255)
                                                                                          ctgcAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTa                                                                                                                                           <  2:575426/67‑1 (MQ=255)
                                                                                          ctgcAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTa                                                                                                                                           >  1:575426/1‑67 (MQ=255)
                                                                                                     ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                             <  2:635139/150‑1 (MQ=255)
                                                                                                           tgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                       >  1:306706/1‑150 (MQ=255)
                                                                                                                  aaGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAgg                                <  1:800301/150‑1 (MQ=255)
                                                                                                                    gCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGt                              >  1:267169/1‑150 (MQ=255)
                                                                                                                         gTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTcga                         >  2:674450/1‑150 (MQ=255)
                                                                                                                                            tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTgcgc      >  2:830706/1‑150 (MQ=255)
                                                                                                                                                gCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATgg  <  2:778088/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGG  >  NC_002947/4240133‑4240426

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