Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 389,998 A→C 14.3% T755P (ACA→CCA)  U712_01810 → Surfactin synthase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881389,9980AC14.3% 104.6 / 2.0 54T755P (ACA→CCA) U712_01810Surfactin synthase subunit 2
Reads supporting (aligned to +/- strand):  ref base A (15/29);  new base C (5/3);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGCAAAGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTG  >  CP006881/389865‑390136
                                                                                                                                        |                                                                                                                                          
tgctgcAAAGGAGGCGGGTTTGTTCGGCACACTGCGCCATTTACTCCTCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGTATACGGCCCGACata                                                                                                                                        <  2:96870/138‑3 (MQ=255)
       aaGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGCGTTTAATTG‑GGAAAGGCTAAGGACCCCCAAAAAACaaga                                                                                                                                >  4:50481/1‑139 (MQ=255)
             gCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTCCTTGTGGAACGGATACGGCCCGCCAGAAAACacgacggtgt                                                                                                                          >  8:268083/1‑135 (MQ=255)
             gCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTCATTGGGGAACGGATAAGGCCCGACAGAAAACACGACGttt                                                                                                                           >  2:12223/1‑138 (MQ=255)
                ggTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAAAGGATAAGGCCCGACACAAAACACCCCGTTTTaaa                                                                                                                       >  8:39473/1‑139 (MQ=255)
                    tGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTAATTGTGGAACGGATACGGCCCGACAGAAAACACCACGTTTTCAACCAg                                                                                                                   >  7:196152/1‑139 (MQ=255)
                    tGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTGCTTGTGGAACGGCAGAGGCCCCCCCGAAAACACCGCCGTTTGAACCAg                                                                                                                   >  2:349458/1‑139 (MQ=255)
                    tGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTCATTGTGGAACGGTTACGGCCCGGCAGAAAACACGCCGTTTTCAACCAg                                                                                                                   >  4:276507/1‑139 (MQ=255)
                           cacaCTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTGATTGTGGGACGGTTACGGCCCGTCAGCAAACCCGCCCTTTTGCACCAGgtttgt                                                                                                             >  4:157982/1‑136 (MQ=255)
                           cacaCTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GGTAATTTTGGGAAGGGTACGGCCCCCCCGAAAACACCACCGTGTGAACCAGTTTTTTg                                                                                                            >  6:202666/1‑139 (MQ=255)
                            acacTCCGCCATTTAATCAGCGCGGGAGATGCCC‑TTGCCCCGCATATGGACAGCAAAGAAAAC‑CAGGCAGCGCCGGCGC‑TTTCATTGTGGAACGGATACGGCCCCCCAGAAAACACGACGTTTTCAACCAGTTTTCTGa                                                                                                           <  4:33583/139‑1 (MQ=255)
                             cacTGCGCAATTTATTCATCGGCGTAGATGCCC‑TTGTCCCGAATATTGTCAGAAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCTTTGTGGAACGGATACGGCCCGACATAAAACACGACGTTTTCAACCAGTTTTCTGAt                                                                                                          <  4:306724/139‑1 (MQ=255)
                             aaatGCGCCATTTATTAATCGGCGGTGAGGCGC‑TTGTCCCGCATATTGTCAGCAAAGTAAAC‑CCGGCATCGGCGGCGG‑TTTCATTGGGCCCCGGTTACGGCCCGAGAGAACACACGACGTTTTCAACCAGTTTTCTGAt                                                                                                          <  3:308063/136‑1 (MQ=255)
                                tGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGCGTTTATTTGTGGAACGGTTACGGCCCCCCTCAACACACGCCGTTTTAACCCAGTTTtgtgtaag                                                                                                        >  2:251288/1‑135 (MQ=255)
                                  cgcCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑GTTAAATGTTGAACGGATACCGCCCGACAGAAAACACGCCCTTTTAAAACAGGTTTATGATCGAcc                                                                                                     >  4:335429/1‑139 (MQ=255)
                                  cgcCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGAATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGAc                                                                                                      <  4:331198/138‑1 (MQ=255)
                                       tttAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAg                                                                                                <  6:160349/139‑1 (MQ=255)
                                       tttAATCATCGGCGG‑TTTGCCCATTGACCCGAATTTGGACAG‑AAATTAAAACCCGCCACCGCGGTCGC‑TTTATTTGTGGAAGGTATCCGGCCCGCCATAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAg                                                                                                <  8:124231/139‑1 (MQ=255)
                                            tcatcGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  3:155917/139‑1 (MQ=255)
                                            gcatcGGCGCAGAGGCCC‑TTGACCCTCATATTGGCAGCAAAGAAAAA‑CTGCCCTCGCCGTCCC‑ATTCATTGTGGACCGCATCCGGCCCGACAGAAACCACGACGGTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  8:61792/138‑1 (MQ=255)
                                            ccatgGGCGGAGAGGCCG‑TTGCCCCGCATATTGACAGCAAAGTAAAC‑CCGGCCTCGCCGTCGC‑TTGAATTGTGGACCGTATCCGCCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  4:217712/138‑1 (MQ=255)
                                             cttcGGCGGAGCGGCCC‑TTGTCCCGCATATTGGCAGAAAAGTAAAA‑CAGGCATCGCGGGCGC‑TTTCATTGTGGACCGGATCCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  8:330530/136‑1 (MQ=255)
                                                cggcggAGAGGCCC‑TTGTCCCGAATTTTGTCCGCAAAGTAAAA‑CCGGCACCCCCGCCGC‑TTTCATTGTGCAACGCATACGGCCCGACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  2:311487/135‑1 (MQ=255)
                                                 ggcggAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGCCGGTTTCAACCAGTTTTCTGATCGACCGCGGGTATggcggcg                                                                                      >  4:174628/1‑138 (MQ=255)
                                                                       cATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGT‑GAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGct                                                                                      <  6:160850/117‑1 (MQ=255)
                                                                       cATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGT‑GAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGct                                                                                      >  5:160850/1‑117 (MQ=255)
                                                                          atTGTAAGAAAAGTAAAA‑CAGGCTGCGCCGGCGC‑TTTCATTGGGGAACGGATCCGGCCCGACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATc                                                               <  8:235814/138‑1 (MQ=255)
                                                                                aGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGGACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAAt                                                         >  1:203110/1‑138 (MQ=255)
                                                                                   aacgTAAAA‑CAGGCTGCGCCGGCGA‑TTTCTTTGTGGAACGGATCCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  6:312499/136‑1 (MQ=255)
                                                                                   aaaGTAAAC‑CAGGCATCGCCGGCGC‑TTTCATTGTGGAACGGATCCGCCCCGACAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  7:268083/139‑1 (MQ=255)
                                                                                   aaaGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  1:12223/139‑1 (MQ=255)
                                                                                    atgTAAAA‑CTGGCATCGGCGTCGA‑TTTCATTGTGTAACGGATACGGCCCGACATAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCAc                                                    <  7:39473/137‑1 (MQ=255)
                                                                                        aaaa‑CTGGCATCGCCGTCGC‑TTTCATTGTGCAACGTATACGGCCCGCCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGc                                                 <  8:8197/138‑1 (MQ=255)
                                                                                                  aTCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACCACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATg                                       >  5:118234/1‑139 (MQ=255)
                                                                                                     gcgggcGC‑TTTCTTTGTGCAACGCATACGCCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATgg                                      <  5:246857/132‑1 (MQ=255)
                                                                                                      cggTCGG‑TTTCATTGTGGAAGGGATCCGGCCCGACCGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATg                                   <  4:163680/137‑1 (MQ=255)
                                                                                                        gTCGC‑TTTCATTGGGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGag                                 <  6:256453/139‑1 (MQ=255)
                                                                                                         ttgtgtttCTTTGTGGAACGGATACGGCCCTAAAAAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTATATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGag                                 <  5:43285/134‑1 (MQ=255)
                                                                                                             atttCATTGTGCACCTGATCCGGCCCCAAACAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGagcagc                             <  4:75568/138‑1 (MQ=255)
                                                                                                              tttCTTTGGGCAACGGAAACGGCCCGACATAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  8:196152/139‑1 (MQ=255)
                                                                                                              gtccATTGTGCAACGGACCCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCGGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  5:83533/136‑1 (MQ=255)
                                                                                                               ttCTTTGTGCAAGGGATCCGGCCCCACAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  3:131047/138‑1 (MQ=255)
                                                                                                                tCAATGTGGAAGGTAAACGGCCCGACACAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  7:231770/137‑1 (MQ=255)
                                                                                                                      tgGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGcc                                                <  7:9281/111‑1 (MQ=255)
                                                                                                                      tgGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGcc                                                >  8:9281/1‑111 (MQ=255)
                                                                                                                        gAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAg                  >  1:131077/1‑139 (MQ=255)
                                                                                                                           cGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCa               <  3:174628/139‑1 (MQ=255)
                                                                                                                                  ggCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGcc                       <  7:171845/124‑1 (MQ=255)
                                                                                                                                  ggCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGcc                       >  8:171845/1‑124 (MQ=255)
                                                                                                                                   gCCCGACAGAAAACACGACGTTTTCACCCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGgcgc       >  1:72578/1‑139 (MQ=255)
                                                                                                                                   gCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCgg          >  4:264221/1‑136 (MQ=255)
                                                                                                                                   gCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCgg          <  3:264221/136‑1 (MQ=255)
                                                                                                                                       gACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCCGCAATGCCTGCAGCCAATCGGCGCGCCt   >  6:293982/1‑139 (MQ=255)
                                                                                                                                       gACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCt   >  1:96943/1‑139 (MQ=255)
                                                                                                                                        aCAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  6:85186/113‑1 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           >  6:246496/1‑50 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  5:246496/50‑1 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            >  5:85186/1‑113 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTg  >  5:74484/1‑139 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTg  >  1:316814/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGCTGCAAAGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCC‑TTGTCCCGCATATTGTCAGCAAAGTAAAA‑CAGGCATCGCCGTCGC‑TTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTG  >  CP006881/389865‑390136

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: ‑