Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,265,530 C→A 12.9% intergenic (‑105/+22) PP_1106 ← / ← quiP hypothetical protein/acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,265,5300CA12.9% 80.4 / 3.3 31intergenic (‑105/+22)PP_1106/quiPhypothetical protein/acyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base A (2/2);  total (15/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 = 5.42e-01

AGTACAGGCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCAACCCAT  >  NC_002947/1265389‑1265661
                                                                                                                                             |                                                                                                                                   
aGTACAGGCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGccc                                                                                                                             <  1:67881/150‑1 (MQ=255)
  tACAGGCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTAATGCGCCCCt                                                                                                                           >  1:551642/1‑150 (MQ=255)
  tACAGGCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTAATGCGCCCCt                                                                                                                           >  1:551593/1‑150 (MQ=255)
       gCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTgaaga                                                                                                                      >  2:610812/1‑150 (MQ=255)
         cTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAAt                                                                                                                    <  1:310695/150‑1 (MQ=255)
                               gTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCaa                                                                                              <  2:369745/150‑1 (MQ=255)
                                  cATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGt                                                                                           >  1:9210/1‑150 (MQ=255)
                                   aTGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTa                                                                                          >  2:418139/1‑150 (MQ=255)
                                                      tGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTAATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACg                                                                       <  2:551642/150‑1 (MQ=255)
                                                      tGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTAATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACg                                                                       <  2:551593/150‑1 (MQ=255)
                                                               ctgTTAAGAGTCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGa                                                              <  2:358174/148‑1 (MQ=255)
                                                               ccGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGa                                                              >  2:36003/1‑150 (MQ=255)
                                                               ccGTTAAGAGGCATGGTGCTGGACGATATATCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGa                                                              >  1:137076/1‑150 (MQ=255)
                                                                  ttAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGt                                                           >  2:101227/1‑150 (MQ=255)
                                                                    aaGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGTGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGTGCCCCTGAAGAATGCCGTTACCGcccc                                                                                                      <  1:137257/105‑1 (MQ=255)
                                                                    aaGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGcccc                                                                                                      >  2:137257/1‑105 (MQ=255)
                                                                    aaGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTc                                                         <  2:137076/150‑1 (MQ=255)
                                                                           aTGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCg                                                                                                                <  2:244495/88‑1 (MQ=255)
                                                                           aTGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCg                                                                                                                >  1:244495/1‑88 (MQ=255)
                                                                                gCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGt                                             <  2:330395/150‑1 (MQ=255)
                                                                                     aCGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAg                                        <  1:371766/150‑1 (MQ=255)
                                                                                      cGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAgg                                       <  1:101227/150‑1 (MQ=255)
                                                                                            ttCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATg                                 >  2:266370/1‑150 (MQ=255)
                                                                                            ttCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATg                                 <  1:610812/150‑1 (MQ=255)
                                                                                               ggatggatTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAg                              >  2:113678/1‑150 (MQ=255)
                                                                                                   ggatTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGcc                          <  2:9210/150‑1 (MQ=255)
                                                                                                                       gAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCaa       >  2:448410/1‑149 (MQ=255)
                                                                                                                       gAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCaa       <  1:448410/149‑1 (MQ=255)
                                                                                                                        aaaaaCTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTATTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCAAcc     >  2:442764/1‑150 (MQ=255)
                                                                                                                           aaCTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCAACCCAt  <  1:418139/150‑1 (MQ=255)
                                                                                                                           aaCTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCAACCCAt  >  1:17974/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                   
AGTACAGGCCTTCGGGGCGGGCGATGACGAGGTCCATGGCACGGGCACTGGTTCTGGGGTTACCCGTTAAGAGGCATGGTGCTGGACGTTAGTTCGGATGGATTTGGGAGGGAGCTGCTGAAAAACTTCATCAAGCCCCCTCATGCGCCCCTGAAGAATGCCGTTACCGCCCCGGCACCAATGTAAGCCGCTGGTTGCCATACGCCCTGCCGAAGTTCTGCTGTTGCAGTGGCAGGCTCATGAAGCGCCCTTTGAACCAGGCATCCAACCCAT  >  NC_002947/1265389‑1265661

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: