Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 883,807 C→A 13.7% intergenic (+45/‑283) PP_0766 → / → PP_0767 hypothetical protein/LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947883,8070CA13.7% 73.8 / 3.5 29intergenic (+45/‑283)PP_0766/PP_0767hypothetical protein/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (10/15);  new base A (2/2);  total (12/17)
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.89e-01

CCCTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGG  >  NC_002947/883659‑883873
                                                                                                                                                    |                                                                  
cccTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTc                                                                    >  7:46333/1‑149 (MQ=255)
 ccTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCg                                                                   >  4:316743/1‑149 (MQ=255)
       gATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCgggg                                                             >  4:367295/1‑149 (MQ=255)
          gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGt                                                          >  8:154420/1‑149 (MQ=255)
          gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCGGGGc                                                            >  4:159253/1‑147 (MQ=255)
          gaagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCACCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCGGGGc                                                            <  3:159253/147‑1 (MQ=255)
           aagaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGt                                                          <  7:154420/148‑1 (MQ=255)
             gaaCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTc                                                       <  7:154052/149‑1 (MQ=255)
               aCGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTa                                                     >  6:12432/1‑149 (MQ=255)
                 gAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTAtt                                                   <  3:316743/149‑1 (MQ=255)
                   agagAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTTCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 <  2:203924/149‑1 (MQ=255)
                   agagAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 >  3:357271/1‑149 (MQ=255)
                    gagaAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTcc                                                >  4:107479/1‑149 (MQ=255)
                               cTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCct                                     >  7:124736/1‑149 (MQ=255)
                                tACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCgggg                                                             <  7:141848/124‑1 (MQ=255)
                                tACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCGAGCCgggg                                                             >  8:141848/1‑124 (MQ=255)
                                  cGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACACCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCtt                                  <  8:113537/149‑1 (MQ=255)
                                   gACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTc                                 <  6:242384/149‑1 (MQ=255)
                                   gACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTGCTATTTCCTGACCCGGCCTCTTc                                 >  2:283210/1‑149 (MQ=255)
                                     cTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTcgc                               <  4:357271/149‑1 (MQ=255)
                                      ttCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTcgcg                              >  3:280690/1‑149 (MQ=255)
                                        cGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCgg                             <  5:142123/148‑1 (MQ=255)
                                        cGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCgg                             <  5:12432/148‑1 (MQ=255)
                                               aGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAAccc                     <  8:124736/149‑1 (MQ=255)
                                                   ccAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCt                  <  8:64392/148‑1 (MQ=255)
                                                   ccAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTc                 <  3:367295/149‑1 (MQ=255)
                                                      gCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTa              <  4:280690/149‑1 (MQ=255)
                                                            ccGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 >  6:342369/1‑108 (MQ=255)
                                                            ccGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTc                                                 <  5:342369/108‑1 (MQ=255)
                                                                   ccAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCgg  <  3:107479/148‑1 (MQ=255)
                                                                                                                                                    |                                                                  
CCCTTGGGATGAAGAACGAAGAGAAACGTGCCTACGACTTCGGCTTCAGTGCCAGCAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGG  >  NC_002947/883659‑883873

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: