Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,672,979 C→A 20.0% V26L (GTG→TTG)  PP_4136 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,672,9790CA20.0% 53.5 / 6.4 25V26L (GTG→TTG) PP_4136LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base A (3/2);  total (13/12)
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.80e-01

GCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGG  >  NC_002947/4672837‑4673121
                                                                                                                                              |                                                                                                                                              
gCCTTGATCACATCAGCCACCCGGGTAGGCACTTCACATTTGCGCAGGATGTTGGATAGGTGAAAATTCACCGTTGACTCCTTGCATTCAAGGATCTGGCCGATCTCCTACGACCTCTTGCCATAAGCACACCACAACAGCACTTGCTGt                                                                                                                                         >  2:71299/1‑150 (MQ=255)
gCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAACTCTTGCGAGAATCACGCCACAACAGCACTTCCTGt                                                                                                                                         >  2:649073/1‑150 (MQ=255)
 ccTTGATCACAGCAGCCACTCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGtt                                                                                                                                        >  1:227550/1‑150 (MQ=255)
      aTCACAGCAGCCACTCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtggt                                                                                                                                   <  2:227550/150‑6 (MQ=255)
           agcagcCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGt                                                                                                                              >  2:808680/1‑150 (MQ=255)
                  aCCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgt                                                                                                                                      >  1:300443/1‑133 (MQ=255)
                  aCCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgt                                                                                                                                      <  2:300443/135‑3 (MQ=255)
                       ggTAGGCACTTCAAATTTGCGCAGGATGTTGTATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCCATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCgg                                                                                                                  <  1:71299/150‑1 (MQ=255)
                                ttCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGc                                                                                                         >  1:837858/1‑150 (MQ=255)
                                         gcgcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtg                                                                                                <  2:837858/150‑1 (MQ=255)
                                             aGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGTATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtgctgc                                                                                            <  1:649073/150‑1 (MQ=255)
                                                                aaTTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCTGGGcc                                                                                                              >  1:789190/1‑112 (MQ=255)
                                                                aaTTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCTGGGcc                                                                                                              <  2:789190/113‑2 (MQ=255)
                                                                 aTTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgtcagt                                                                        >  1:187812/1‑150 (MQ=255)
                                                                    cACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCAt                                                                     >  2:170068/1‑150 (MQ=255)
                                                                     aacggagACTTCTTGCATTCAAGGATCTGGCCTGTCTCCCACGGCCGTTTTCGATAAACACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGGTGCGGCGAGTGTCGAGCGTTGTCATTCATg                                                                    >  2:335701/7‑150 (MQ=255)
                                                                        gTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGttt                                                                 <  1:808680/150‑1 (MQ=255)
                                                                           gACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                              <  2:505034/150‑1 (MQ=255)
                                                                                             aTCTGCCCGATCTCTCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGc                                            <  1:335701/150‑1 (MQ=255)
                                                                                                   ccGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCAc                                                        >  1:614441/1‑132 (MQ=255)
                                                                                                   ccGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCAc                                                        <  2:614441/132‑1 (MQ=255)
                                                                                                        ctcCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgct                                 <  2:221833/150‑1 (MQ=255)
                                                                                                               gAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTAcc                          >  2:1838/1‑150 (MQ=255)
                                                                                                                                      caacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAg   <  1:513920/150‑1 (MQ=255)
                                                                                                                                       aacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCAAGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAgg  >  1:479998/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGG  >  NC_002947/4672837‑4673121

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: