Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 389,046 T→G 17.1% R437R (CGT→CGG U712_01810 → Surfactin synthase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881389,0460TG17.1% 51.4 / 2.2 43R437R (CGT→CGGU712_01810Surfactin synthase subunit 2
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base G (6/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCACCG  >  CP006881/388934‑389182
                                                                                                                |                                                                                                                                         
tCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGGGGGGTTACGGGCTTGTGGGAGCCGaa                                                                                                                   >  4:351564/1‑137 (MQ=255)
tCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTGCGTGAGGTCGGGCGTTGGGGGGACGAagc                                                                                                                 >  7:34809/1‑138 (MQ=255)
 cAAAATAAAATATAACGGAGCCGCCTTTGATTCCGCTGTTATCGACCGGACAGCGGAGCACCATACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGGGAGTACGGGCTTGTGGGAGACGAagag                                                                                                                <  2:105142/139‑1 (MQ=255)
         aaTATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCa                                                                                                       <  3:222821/139‑1 (MQ=255)
           tataACGGAGCCGCCGTTGATTCCGCTTTTATCGAACGGACCGCGGAGCACCTAACCCGCATGATGGAAGCAGCTGGCGATCAGCCGGCCGCTTTTGTGCGGGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCaaa                                                                                                     <  8:364562/139‑1 (MQ=255)
                cGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTTTGGGGGCGTATGGGGGTGGGGGGGGCGAAGCGGCGCGGGCAAAtttg                                                                                                 >  1:232347/1‑137 (MQ=255)
                 ggAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                               <  1:18836/139‑1 (MQ=255)
                 ggAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                               <  8:300601/139‑1 (MQ=255)
                 ggAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                               <  8:146736/139‑1 (MQ=255)
                 ggAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                               <  4:430889/139‑1 (MQ=255)
                 ggAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGGGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                               <  8:301493/139‑1 (MQ=255)
                   aGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAg                                                                                             <  3:236302/139‑1 (MQ=255)
                     ccgccTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTGTTTTGCGTGAGTACCGGGGTGTGGGGGCCGAAGAGCCCCGGGAAAATGtgcaggg                                                                                           >  1:125262/1‑133 (MQ=255)
                        ccgttGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCGGTCGATCAGCCGGCCGCGTTTTTGGGGGGGTTAGGGGGGGGGGGGGGCGAAGAGGCGCGGCGAATTGTCGAGGTGtt                                                                                        >  4:354409/4‑139 (MQ=255)
                           ttGATTGCGCTTTTTTCGAACGGACAGCGGAGCACCGTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTTTGGGGGTGTTGGGGGGTGTTGGGGGCGCAGGGGCGCGGCCAAATTTGGAGGTGtttt                                                                                      >  6:376817/1‑137 (MQ=255)
                             gcttCCGCTTTTCTCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAAc                                                                                    <  3:295595/136‑1 (MQ=255)
                               ttCCGCTTTTATCGAACGGACAGCGGCGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTGCGGGTGTACGGGCGTGTGGGGGACGAAGAGCCGCGGCCAATTGGGGGGGGGTTTAACAGc                                                                                 >  8:326289/1‑139 (MQ=255)
                                  cGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTTCGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACg                                                                              >  6:73789/1‑139 (MQ=255)
                                       ttATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGaaa                                                                           <  5:73789/137‑1 (MQ=255)
                                         aTCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTGCGTGAGTAAGGGCTTGTGGGAGACGAAGAGCCGCGGCAAATTTTAGAGGTGTTTTACAGCACGAAAGCCg                                                                       >  2:102232/1‑139 (MQ=255)
                                           cGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGaa                                                                     <  3:80248/139‑1 (MQ=255)
                                                 gACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCt                                                               <  1:322342/139‑1 (MQ=255)
                                                      cGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTTTGGGTGAGTTAGGGCGTGGGGGGGTCGCAGAGCAGCGGCCAATTTTAGAGGGGTTTAAACTCACGAAAGCCGAACACCCAGAAgg                                                          >  4:80248/1‑139 (MQ=255)
                                                        gAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGaa                                                        <  1:281812/139‑1 (MQ=255)
                                                            aCCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGc                                                    <  3:215878/139‑1 (MQ=255)
                                                              cTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGGGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTg                                                  <  5:6860/139‑1 (MQ=255)
                                                                tACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTTGGTGAGTAAGGGCTTGTGGGTGTCGAAAAGCAGCGGCAAATTGTGGAGGTGGTTTAAAGCCCGGAACGCCAACAGCCTGAAGTAATGTCTGt                                                 >  4:433618/1‑138 (MQ=255)
                                                                   ccGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCAc                                             <  3:354409/139‑1 (MQ=255)
                                                                     gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTTTTCGTGAGTACGGGCTTTTGGGAGAAGAAGACCACCAGCAAAATTTAGAGGTGTTTAACAACACGAAATCCGAACTCCCTGAAGTAATTGCTTTTAAcc                                            >  6:6860/1‑138 (MQ=255)
                                                                     gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCa                                           <  1:102232/139‑1 (MQ=255)
                                                                     gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCa                                           <  8:34809/139‑1 (MQ=255)
                                                                     gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTAAGGGCTTGTGGGAGCCGAAGAGCAGCGGCAAAATGTCGAAGTGTTTTAAAGCACGCAAGCCGAACACCCAGAAGGAAAGGCGGttaccc                                            >  4:309439/1‑134 (MQ=255)
                                                                     gCATGATGGAACCAGCTGGCGAGCAGCCGGCCGCTTTTGTGCGTGAGTCCGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCa                                           <  6:197036/139‑1 (MQ=255)
                                                                        tgatgGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGGGGGAGTACGGGGGTGTGGGAGACGCAGAGCAGCGGCAAATTTTCGCGGTGTTTAACAGCACGAAAGCCGCACACCCGGAAGGGATGGCTGGTCACCAAgt                                        >  4:296177/1‑139 (MQ=255)
                                                                        tgatgGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTa                                                                                                                                      <  1:255327/46‑1 (MQ=255)
                                                                        tgatgGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTa                                                                                                                                      >  2:255327/1‑46 (MQ=255)
                                                                        tgatgGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAagag                                                                                                                <  7:348902/68‑1 (MQ=255)
                                                                        tgatgGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAagag                                                                                                                >  8:348902/1‑68 (MQ=255)
                                                                                gcagcTGTCGATCAGCCGGCCGCGTTTGTGCGTGAGGACGGGGGTGTGGGTGGCGAAGAGCCGCGGCAAATTGTCGCGGTGTTTTACAGCACGAAAGCCGAACTCCCGGCAGGAATGGCTGTTCACCAAGGGGTTGGag                                >  2:293654/1‑139 (MQ=255)
                                                                                  agcTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAATTTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAagag                              <  3:309439/139‑1 (MQ=255)
                                                                                              gccggccgCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCt                                                               <  2:234201/94‑1 (MQ=255)
                                                                                              gccggccgCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCt                                                               >  1:234201/1‑94 (MQ=255)
                                                                                                              cGTGAGTACGGGCTTGGGGGAGACGACGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccg  <  7:326289/139‑1 (MQ=255)
                                                                                                                |                                                                                                                                         
TCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCACCG  >  CP006881/388934‑389182

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