Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,672,979 C→A 18.4% 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,9790CA18.4% 55.1 / 7.6 28V26L (GTG→TTG) PP_4136LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base A (2/3);  total (12/16)
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.38e-01

GCGAATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGG  >  NC_002947/4672830‑4673121
                                                                                                                                                     |                                                                                                                                              
gCGAATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAc                                                                                                                                                >  2:245491/1‑150 (MQ=255)
  gAATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACtt                                                                                                                                              >  1:193871/1‑150 (MQ=255)
            gATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtgg                                                                                                                                    >  2:259665/1‑146 (MQ=255)
             aTCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGt                                                                                                                                   <  2:293757/150‑1 (MQ=255)
             aTCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtggt                                                                                                                                   <  1:259665/150‑6 (MQ=255)
               cacaGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTgg                                                                                                                                 >  2:534001/1‑150 (MQ=255)
               cacaGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtggtgt                                                                                                                                 <  1:263909/150‑8 (MQ=255)
                 cagcagCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCg                                                                                                                               >  1:614750/1‑150 (MQ=255)
                    cagcCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACCTGGCGTCa                                                                                                                            <  1:119650/150‑1 (MQ=255)
                               gTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                 <  1:430350/150‑1 (MQ=255)
                                     aCTTCAAATTTGCGCGGGATGTTGGATATGTTAAAATGCACCGTTGTCTCCTTGCATTCAAGGATTTGGCCGATCTCGCACGAGCTATTCCCATAAGCACACAACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTcc                                                                                                           <  2:25386/150‑1 (MQ=255)
                                        tCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtg                                                                                                                                     <  2:36788/121‑4 (MQ=255)
                                        tCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtg                                                                                                                                     >  1:36788/1‑118 (MQ=255)
                                            aTTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                    <  1:245491/150‑1 (MQ=255)
                                                      gATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGa                                                                                          <  1:546898/150‑1 (MQ=255)
                                                                                tCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTg                                                                >  1:495281/1‑150 (MQ=255)
                                                                                  gACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                              >  2:119201/1‑150 (MQ=255)
                                                                                              ccAAGGATCTGCCTGTTCTCCTACGCGCTTTTGGCATAAGCACACCACAACAGCATTTGCTGTTCACGTGGCGTCAGGCGAACCGGGTCGTCCGCCAGTGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTCAACTCCACGGGTGa                                                  <  2:193871/149‑1 (MQ=255)
                                                                                                       tGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                         <  1:119201/150‑1 (MQ=255)
                                                                                                       tGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                         >  1:402140/1‑150 (MQ=255)
                                                                                                                    aCGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGTCGAACCggg                                                                                                                 >  1:352050/1‑65 (MQ=255)
                                                                                                                    aCGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                 <  2:352050/65‑1 (MQ=255)
                                                                                                                          tctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                      <  2:402140/150‑1 (MQ=255)
                                                                                                                                       acacCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTCCGaa         <  2:614750/150‑1 (MQ=255)
                                                                                                                                          ccacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca      >  2:524687/1‑150 (MQ=255)
                                                                                                                                           cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGt                                                                   <  2:525087/88‑1 (MQ=255)
                                                                                                                                           cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGt                                                                   >  1:525087/1‑88 (MQ=255)
                                                                                                                                           cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGcac     >  1:550484/1‑150 (MQ=255)
                                                                                                                                              aacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAgg  <  2:302627/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCGAATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGG  >  NC_002947/4672830‑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: