Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,204,721 T→G 11.1% K274T (AAA→ACA)  U712_11055 ← Stage IV sporulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,204,7210TG11.1% 129.3 / 2.7 58K274T (AAA→ACA) U712_11055Stage IV sporulation protein A
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base G (5/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATATAGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTAGAGGGCCTCTC  >  CP006881/2204586‑2204853
                                                                                                                                        |                                                                                                                                    
tCATATAGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTaa                                                                                                                                   <  7:54709/139‑1 (MQ=255)
tCATATAGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTa                                                                                                                                    <  1:161385/138‑1 (MQ=255)
  atatAGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAat                                                                                                                                  <  3:14454/138‑1 (MQ=255)
   tataGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAata                                                                                                                                 >  3:55468/1‑138 (MQ=255)
   tataGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAata                                                                                                                                 >  2:241510/1‑138 (MQ=255)
    ataGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGATGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTAgtat                                                                                                                                >  5:14374/1‑138 (MQ=255)
            tCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCg                                                                                                                       >  1:85462/1‑139 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGCAAACCCTGTCTACGTCCCGGAGCCGTTTAATAAACTTAACCga                                                                                                                      >  7:364616/1‑138 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACACCCCTGTCTACGTCCCGGAGCCGTTGAATATGCTTagaact                                                                                                                      >  1:234258/1‑134 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTa                                                                                                                                    >  7:31547/1‑125 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTa                                                                                                                                    <  8:31547/125‑1 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAAACGt                                                                                                                      >  6:372972/1‑139 (MQ=255)
             cAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTAATATCCTTAACCGt                                                                                                                      >  1:86446/1‑139 (MQ=255)
              aGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGt                                                                                                                      <  1:249391/138‑1 (MQ=255)
                gcgcGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGt                                                                                                                      <  7:363072/136‑1 (MQ=255)
                gcgcGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATACTTAACCGt                                                                                                                      >  8:363072/1‑136 (MQ=255)
                         aaTCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGa                                                                                                          <  1:27162/139‑1 (MQ=255)
                         aaTCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGa                                                                                                          <  4:55468/139‑1 (MQ=255)
                                           ctaggCCCAGCTCAATTCCGGCGAATCCCGCACTTTCAATGAATTCACACTCGCTGACTTGGCCGACAACCCTGGCTACGTCCCGGAGCCGCTAAAAATCCTAAACCGACACCTTCACGGACTCCTGATAGCTATCACg                                                                                        <  6:239074/136‑1 (MQ=255)
                                                   aGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTAATATCCTTAACCGGTTGCTTCACGGACTCCTGATATCTTTAACGCAACCa                                                                                  >  6:378580/1‑137 (MQ=255)
                                                   aGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGACGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAAt                                                                                <  8:295734/139‑1 (MQ=255)
                                                    gCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATg                                                                               >  1:27675/1‑139 (MQ=255)
                                                     cTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATgg                                                                              >  8:124582/1‑139 (MQ=255)
                                                     cTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTAACGCAACCAATgg                                                                              >  7:288259/1‑139 (MQ=255)
                                                     cTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTACTTCACGGCCGCCTGGTAGCGTTCCCGCCACCAAAgg                                                                              >  7:214529/1‑139 (MQ=255)
                                                               gCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAat                                                                                                                                  >  6:125059/1‑77 (MQ=255)
                                                               gCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAat                                                                                                                                  <  5:125059/77‑1 (MQ=255)
                                                               gCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCATAACCCTTTCCTTCACGGACGCCTGATAGCTTTAACGCAACCAATtgttttgttt                                                                      >  7:206345/1‑137 (MQ=255)
                                                                cTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAg                                                                   <  6:362056/139‑1 (MQ=255)
                                                                   aTCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTTACCGTTTCCTTCACGGACGCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCAc                                                                >  5:139568/1‑139 (MQ=255)
                                                                    tacGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTTAATTCACAAAAAGTTTCAATGAAGGCCTCCAGAAAGAATTCACACAAACAATGGGTATTTTTCAacaca                                                               >  2:99559/3‑138 (MQ=255)
                                                                    tCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTTATATCCATAAACGGTTACTTAACGGACGCCTGATAGCTTTAACGCAACCAATTGTTTTGTTTCAGCAcc                                                               >  4:305379/1‑139 (MQ=255)
                                                                    tCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCAc                                                                >  8:37898/1‑138 (MQ=255)
                                                                        gCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTa                                                           <  6:33119/139‑1 (MQ=255)
                                                                          acattCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTAcc                                                         <  2:9347/136‑1 (MQ=255)
                                                                          aCTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTAcc                                                         <  5:378580/139‑1 (MQ=255)
                                                                          aCTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTAcc                                                         <  6:55191/139‑1 (MQ=255)
                                                                          aCTTTCAATGAAATCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTAcc                                                         <  8:364616/139‑1 (MQ=255)
                                                                                      ttCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTAATATCCTTAACCGGTTGCTTCACGGCCACCGGATAGCTTTAACCCAACCAATGGTTTTATTTAAGCACCATTACCCAGCTTgg                                                 >  5:94590/1‑135 (MQ=255)
                                                                                            cTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACa                                       <  2:86446/139‑1 (MQ=255)
                                                                                               gCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACAt                                      <  1:241510/137‑1 (MQ=255)
                                                                                               gCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACAt                                      <  7:37898/137‑1 (MQ=255)
                                                                                                cTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTTCCCAGCTTGGGAGATTGACATTCa                                   >  7:31732/1‑139 (MQ=255)
                                                                                                cTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATACATAACCGGTTACGTCCCGGTCTCCTGATAGCTTTCACGCACCCAAAGGTTTTGTTTAAGCCCCATTACCCAGCTTGGGAGATTTACAGTCa                                   >  4:338582/1‑139 (MQ=255)
                                                                                                cTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTGTTATATACTTAACCGTTTCCGTCACGGACGCCTGATAGCTTTAACGCAACCAATGGTTTTGTTTCAGCCCCAATTACCAGCGTGGGGGATTGTCAttaa                                   >  5:234330/1‑137 (MQ=255)
                                                                                                              aaCCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGa                     <  5:302760/139‑1 (MQ=255)
                                                                                                              aaCCCTGTCTACGTCCCGGAGCCGTGTAATATCCTTAACCGTTTCCATCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGa                     <  6:14374/139‑1 (MQ=255)
                                                                                                              aaCCCTGTCGACGTCCCGGAGCCGATTAATATCCTTAACCGTGTCCGTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGa                     <  2:287961/139‑1 (MQ=255)
                                                                                                                   tGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGGTTCCTTCACGGACTCCTGATAGCTTTAACGCAACCAAAGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTc                >  8:119713/1‑139 (MQ=255)
                                                                                                                      cTACGGCCCGGAGCCGTTTAATAACCTTACCCGTTACCTTCACGGACTCCTGATAGCTTTCCCGAAACCAATGGTTTTATTACAGCACAATTACCCTGCTTGGGAGATTGAAATTCAATTCTAGCACCGGAAACGCGTa             <  5:372972/139‑1 (MQ=255)
                                                                                                                       tACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTGCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTAAGCACCATTACCCAGCGTGGGAGCTTGACATTCACTTCTAGCACCGGAAACTCGTag            >  8:70825/1‑139 (MQ=255)
                                                                                                                         cGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTaga           <  7:124582/138‑1 (MQ=255)
                                                                                                                         cGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTa             <  3:305379/136‑1 (MQ=255)
                                                                                                                         cGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCATCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTaga           <  6:234330/138‑1 (MQ=255)
                                                                                                                          gTCCCGGATCCCTGTAATATCCGTAACCGTTTCCTTCACGGACTCCTGCTAGCTTGCACGCAACCAATTGTGTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCAAGCACCGGAAACTCGTAGAgg         <  2:372719/139‑1 (MQ=255)
                                                                                                                           tCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTAGAggg        <  5:351401/139‑1 (MQ=255)
                                                                                                                            cccGGAGGCGGTTAATATCCATAACCGTATCCTTCACGGTCTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTAGAGGGc       <  2:234258/139‑1 (MQ=255)
                                                                                                                                 aGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTAGAGGGCctctc  <  8:206345/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TCATATAGATGATCAGGCGCGTACAAATCAATTTCTGCCACCCCTTGGCCCAGCTCAATTCCGGCTAATCCGGCACTTTCAATGAATTCAAACTCGCTGAATTGGCCGACAACCCTGTCTACGTCCCGGAGCCGTTTAATATCCTTAACCGTTTCCTTCACGGACTCCTGATAGCTTTCACGCAACCAATGGTTTTCTTTCAGCACCATTACCCAGCTTGGGAGATTGACATTCACTTCTAGCACCGGAAACTCGTAGAGGGCCTCTC  >  CP006881/2204586‑2204853

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