Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 6,089,343 A→C 18.5% S192R (AGC→CGC)  paiB → Protease synthase and sporulation protein PAI 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1436,089,3430AC18.5% 93.3 / 7.4 59S192R (AGC→CGC) paiBProtease synthase and sporulation protein PAI 2
Reads supporting (aligned to +/- strand):  ref base A (35/12);  new base C (2/10);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 4.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTGGCACGCCCTTTGGC  >  ACB143/6089212‑6089474
                                                                                                                                   |                                                                                                                                   
cggCCGACTACCGCGACCGAATGCTCCGCGCCCTCGTCGGCAGCCGCCTGCCCATCGCGCGCCTGCACGGTGCGCGCAAGCTCCGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGCGCGCTAAGGCGTCAgccc                                                                                                                                 <  6:50823/134‑3 (MQ=255)
 tGCCGACTACATCGAGGGCATGCTCCGCCCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAgcgc                                                                                                                                 >  3:262024/1‑135 (MQ=255)
         aCATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGa                                                                                                                        >  7:252564/1‑136 (MQ=255)
         aCATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCAGACCATCGCCCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGACTGAGCGCCAGCCACg                                                                                                                         >  8:182031/1‑135 (MQ=255)
             cGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCaa                                                                                                                    >  5:288030/1‑136 (MQ=255)
                   cATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCCGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGa                                                                                                               >  1:20864/1‑135 (MQ=255)
                   cATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCACCTGGa                                                                                                               >  2:120992/1‑135 (MQ=255)
                     tGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACaa                                                                                                            >  7:162710/1‑136 (MQ=255)
                               cATCGTCGGTATCCGTCTGCCCATCTCGCGCCTGCAAGGTGCGCGCAAGCTCAGACAGAACCGATCCGGGCAAGCCATCGCTGGTGTGAGCGAAGGCCTGCGCGCCAGCCCCGAGCACCTGGACACTCAACTcgcg                                                                                                  >  1:110806/1‑136 (MQ=255)
                               cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCCGCCCCGATCAACTGGACAATCAACTcgcg                                                                                                  >  5:62207/1‑136 (MQ=255)
                               cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTcgcg                                                                                                  >  8:210405/1‑136 (MQ=255)
                               cATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTcgcg                                                                                                  >  3:163789/1‑136 (MQ=255)
                                    tcgGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGca                                                                                             >  2:213398/1‑136 (MQ=255)
                                     cgGTCTCCGCCTCCCCATCGCGCGCCTGCAAGGTGCGCTCACGCTCAGCCAGAACCGCGCCTGGCAATACATCGCTGGTGTGCGCGATGGCCTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGgcgc                                                                                              <  8:96282/134‑1 (MQ=255)
                                             gTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGCGCGCTAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGcc                                                                                     >  3:48294/1‑135 (MQ=255)
                                                tGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAAct                                                                                 >  8:99628/1‑136 (MQ=255)
                                                tGCCCATAGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCACCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTCGACAATCAACTCGCGGCGCACATGCGCCAAca                                                                                 >  1:293281/1‑135 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACCTGCGCCAACTCTg                                                                              >  4:27656/1‑134 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTg                                                                              >  7:136233/1‑134 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTg                                                                              >  7:63064/1‑134 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTg                                                                              >  7:93980/1‑134 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTg                                                                              >  1:199173/1‑134 (MQ=255)
                                                     aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTg                                                                              >  4:119637/1‑134 (MQ=255)
                                                        gcgggcCTCCACGGGGCGCGCAAGCTCCGCCAGCACCGCTCCGGGCACGACTTCGCGGGGGGGCGCGACGCCCGCCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAg                                                                           <  4:121210/134‑1 (MQ=255)
                                                          gcgcCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGaa                                                                        >  3:111313/1‑135 (MQ=255)
                                                          gcgcCTGCAAGGTGAGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCAGGTGTGCGCGAAGGCCTGAGCGACAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCACCTCTGAAGGaa                                                                        >  8:245900/1‑135 (MQ=255)
                                                           agcCTGCCAGGTGAGGGCAAGCTCACCCACACCGGCTCGGGGCTATCCTTGGCGTGGGTGCGGGAAGGCCTGACCGCCAGCCCAGATCAACTGGACAATCAACTCGCGGCGCAAATGCGCCAACTCTGAAGGAATa                                                                      <  8:84969/135‑1 (MQ=255)
                                                            tccTCCATGGTGGGGGCAAGCTCCGCCATACCCGCTCCGGCCAAGACAGGGGGGGGGGGGGCGACGGCCTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAg                                                                     <  4:189718/135‑1 (MQ=255)
                                                            gcCTGGAAGGGGCGCGCAAGCTCACCCCGAACCGCTCCGGCCAAGACAGCGCTGGTGTGCGCGACGGCCTGAGCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAg                                                                     <  2:53502/136‑1 (MQ=255)
                                                               ggcttggTGCGCGCAAGCTCACCCAGACCCGCGCCGTGCAAGACTTCGCTGTGGTGCGCGACGGCCTGCGCGCCAGCCCCGATCAACTTGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGccc                                                                  <  3:44347/131‑1 (MQ=255)
                                                                 ctgggTGCGGGACAGCTCAGCCAGCACCGCTCCGGGCAAGACTGCGCGGGGGTGCGCGCCGCCCGCACCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                                                                <  5:3368/133‑1 (MQ=255)
                                                                 cacggggcgcgcACGCCCCGCCAGACCCGCTCCGGGCACGCCCTCGCGGGTGCGCGCGCCCCCCCCCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                                                                <  7:210405/130‑1 (MQ=255)
                                                                 cacggggcgcgcAACCTCAGCCAGAACCGCTCCGTGCAAGCCAGCGCGGGTGTCCCCCCCGGCGCCCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                                                                <  6:111503/130‑1 (MQ=255)
                                                                 caagcggcgctcAAGCTCAGCCACAACCGCTCCGGTCAAGACTTCGCTGGTGTGCGCTAAGGCCGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                                                                <  8:103691/130‑1 (MQ=255)
                                                                        cgcgcCAGCTCACCCCGACCCGCTCCGGGCAAGACATCGGGGGGGGGCGCGACGGCCGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAc                                                         <  3:191445/136‑1 (MQ=255)
                                                                         gcgcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAcc                                                        >  8:27179/1‑136 (MQ=255)
                                                                           gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGTACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCt                                                       >  7:24308/1‑135 (MQ=255)
                                                                                       aGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCt                                          >  8:163561/1‑136 (MQ=255)
                                                                                        gacccGCTCCGGGCATGACAGCGCTGGTGTGCGCTAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTg                                         <  6:154068/133‑1 (MQ=255)
                                                                                        gAACCGCTCCGGGCACGACTTCGGTGGGGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTg                                         <  1:220090/136‑1 (MQ=255)
                                                                                            cGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCAc                                     >  8:268035/1‑136 (MQ=255)
                                                                                            cGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCAc                                     >  6:26991/1‑136 (MQ=255)
                                                                                             gCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACg                                    >  1:76649/1‑136 (MQ=255)
                                                                                                ccGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGcc                                 >  8:154986/1‑136 (MQ=255)
                                                                                                 cGGGCAAGACATCGGTGGTGTGCGCGAAGGCCGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGccc                                <  2:285124/136‑1 (MQ=255)
                                                                                                 cGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGccc                                >  2:297300/1‑136 (MQ=255)
                                                                                                 cGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGcc                                 >  6:10343/1‑135 (MQ=255)
                                                                                                 cGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCACCCTGTCACGGcc                                 >  4:92668/1‑135 (MQ=255)
                                                                                                 cGGGCAAGAAATCGCTGGTGTGCTCGAAGGCCTGAGCGCCAGCCCCGATAAACTGGACAATCAACTCTCGGCGCACATGAGCCAACTCTGAAGGAATAGCCCATGCCCACCCTTACCCTGCGCCCTGTCACGGcc                                 >  2:98294/1‑135 (MQ=255)
                                                                                                  gggCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGccc                                >  6:272290/1‑135 (MQ=255)
                                                                                                  gggCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCa                               >  2:102214/1‑136 (MQ=255)
                                                                                                       aGACCTCGGTGGGGGGCGCGAAGGCGGGACCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATc                          <  4:172898/136‑1 (MQ=255)
                                                                                                       aGACATCGCTGGTGGCCGCGACCGCCGCCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATc                          <  7:8378/136‑1 (MQ=255)
                                                                                                          cATCGCTGGGGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACg                       <  2:249307/136‑1 (MQ=255)
                                                                                                           aTCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACg                       >  7:122991/1‑135 (MQ=255)
                                                                                                                 ggggggcgcGAAGGCCTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTg                <  8:252564/131‑1 (MQ=255)
                                                                                                                 ggTGGGCTCGAAGGCCTTAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTg                <  2:87174/136‑1 (MQ=255)
                                                                                                                      gcgcgcAGGCCCGACCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTGGCACg           <  6:62207/136‑1 (MQ=255)
                                                                                                                               ccTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTGGCACGCCCTTTGGc  <  3:66883/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
CTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCTGTCACGGCCCAAGATCACGCAGCCTGGCACGCCCTTTGGC  >  ACB143/6089212‑6089474

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