Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,162,657:1 +T 13.3% intergenic (+132/+69) PP_4544 → / ← PP_4545 hypothetical protein/beta‑ketodecanoyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,162,6571.T13.3% 137.4 / 20.9 45intergenic (+132/+69)PP_4544/PP_4545hypothetical protein/beta‑ketodecanoyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base . (16/23);  new base T (2/4);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAG  >  NC_002947/5162510‑5162806
                                                                                                                                                    |                                                                                                                                                     
ggAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                                     <  1:905445/150‑1 (MQ=255)
  gcagggcctgttccgcagtgggcatgacgtttttgttcaccaggatttccgaacgcaggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggccttcgcgggcacgcccgctcccacaggaaCACCACAacctttga                                                                                                                                                  >  2:1136499‑M2/136‑147 (MQ=255)
  aaGCCAAGTTCTAAGCAACCCCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                                   <  2:1249742/150‑1 (MQ=255)
     ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                                <  2:128613/150‑1 (MQ=255)
     ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                                >  2:1187552/1‑150 (MQ=255)
         gTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtgt                                                                                                                                            <  2:1155074/150‑1 (MQ=255)
               cgcagtgggcatgacgtttttgttcaccaggatttccgaacgcaggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggccttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTa                                                                                                                                      <  1:1136499‑M2/28‑1 (MQ=255)
                    cacCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                                   <  2:829999/132‑1 (MQ=255)
                    cacCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                                   >  1:829999/1‑132 (MQ=255)
                     gcatgacgtttttgttcaccaggatttccgaacgcaggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggccctttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTACCtgtg                                                                                                                                <  1:517440‑M2/34‑1 (MQ=255)
                            tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgg                                                                                                                         >  1:324372‑M2/1‑138 (MQ=255)
                             ctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcggg                                                                                                                        >  2:415037‑M2/1‑137 (MQ=255)
                                ctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaa                                                                                                                     >  2:805113‑M2/1‑134 (MQ=255)
                                 ttCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCAAATTAGCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaag                                                                                                                    <  1:1243229‑M2/150‑18 (MQ=255)
                                       tGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgg                                                                                                                               <  2:1023692‑M2/133‑7 (MQ=255)
                                       tGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgg                                                                                                                               >  1:1023692‑M2/1‑127 (MQ=255)
                                           gCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAaga                                                                                                          <  1:716400/150‑1 (MQ=255)
                                            cGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAagag                                                                                                         <  2:195728/150‑1 (MQ=255)
                                                 aaCGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGt                                                                                                    <  2:550054/150‑1 (MQ=255)
                                                    gCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgc                                                                                                               >  1:1203165‑M2/1‑114 (MQ=255)
                                                    gCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgc                                                                                                               <  2:1203165‑M2/136‑23 (MQ=255)
                                                    aacgcaggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggcctttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGccc                                                                                                                 <  2:1267835‑M2/49‑1 (MQ=255)
                                                    aacgcaggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggcctttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGccc                                                                                                                 >  1:1267835‑M2/86‑134 (MQ=255)
                                                        aCCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGgc                                                                                             <  1:1187552/150‑1 (MQ=255)
                                                          aggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggccttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGcat                                                                                                                     <  2:768273‑M2/45‑2 (MQ=255)
                                                          aggaccatggccgttcctctgtgtttattcagtggggtaggggtttgggccttcgcgggcacgcccgctcccacaggaaCACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGcat                                                                                                                     >  1:768273‑M2/80‑123 (MQ=255)
                                                              tCCGGCGTTTCCCGTTTCTTTTTGTCTGTGTGCACGGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtga                                                                                       <  1:6461/150‑1 (MQ=255)
                                                              tCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtga                                                                                       >  1:752501/1‑150 (MQ=255)
                                                               ccGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGTCCCCATCGCGGGCTTGTCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                                      <  2:875295/150‑1 (MQ=255)
                                                                     tttCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCa                                                                                <  2:752501/150‑1 (MQ=255)
                                                                     tttCCCGCTTCTTTTTGTCCGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGCCCTTGAAAACTGGTCTGTACCTGTGGGAGCGGCCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCa                                                                                <  1:1242118/150‑1 (MQ=38)
                                                                      ttCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCATAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAg                                                                               >  1:80573/1‑150 (MQ=255)
                                                                                                aTGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAAt                                                     <  2:282399/150‑1 (MQ=255)
                                                                                                          cccTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGTCCGTAACAGgcgc                                                                                           >  2:127780/1‑102 (MQ=255)
                                                                                                          cccTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGTCCGTAACAGgcgc                                                                                           <  1:127780/102‑1 (MQ=255)
                                                                                                                gcgGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCa                                     >  2:752802/1‑150 (MQ=255)
                                                                                                                   ggCTTGCCCGCTCCCAACGTATCACACCAGAGCTTCAGAAATGTGATGTACGGGTGGGAGCTGGCATGCCCGCGAAGAGGCCGTAACAGGCGAATGATGCATCAGCGCTCGAGAAGAAACATACTACCAATCGAGCACCCCGCACCAAAc                                  <  2:1243229/150‑1 (MQ=255)
                                                                                                                          ccGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgcgaaggggccggcgcatttcacacagacaaaaagaagcgggaaacgccggatggggtcggcgtttcccgcttcaggtgaagcaga                           <  2:677527‑M2/150‑107 (MQ=255)
                                                                                                                          ccGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGTGAACGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                           >  1:4036/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAAt                                                     >  1:678013/1‑124 (MQ=255)
                                                                                                                          ccGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAAt                                                     <  2:678013/124‑1 (MQ=255)
                                                                                                                          ccGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                           >  1:68173/1‑150 (MQ=255)
                                                                                                                            gCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgcgaaggggccggcgcatttcacacagacaaaaagaagcgggaaacgccggatggggtcggcgttt                                               <  2:1145070‑M2/128‑87 (MQ=255)
                                                                                                                            gCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgcgaaggggccggcgcatttcacacagacaaaaagaagcgggaaacgccggatggggtcggcgttt                                               >  1:1145070‑M2/1‑42 (MQ=255)
                                                                                                                              tCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAacac                       <  1:752802/150‑1 (MQ=255)
                                                                                                                                          accacAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGcc           >  1:233934/1‑150 (MQ=255)
                                                                                                                                           ccacAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCt          <  2:80573/150‑1 (MQ=255)
                                                                                                                                                   cTTGAGAACTGTGCTGTACctgtgggagcgggcaagcccgcgaaggggccggcgcatttcacacagacaaaaagaagcgggaaacgccggatggggtcggcgtttcccgcttcaggtgaagcagagcagggtgttgcttagaacttggct  >  1:677527‑M2/1‑19 (MQ=255)
                                                                                                                                                   cTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAg  >  2:698486/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAG  >  NC_002947/5162510‑5162806

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