Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 6,084,764 A→C 29.1% N178T (AAC→ACC)  paiB → Protease synthase and sporulation protein PAI 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1226,084,7640AC29.1% 77.4 / 2.4 52N178T (AAC→ACC) paiBProtease synthase and sporulation protein PAI 2
Reads supporting (aligned to +/- strand):  ref base A (26/9);  new base C (4/13);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 8.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCT  >  ACB122/6084631‑6084896
                                                                                                                                     |                                                                                                                                    
cACGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGCAc                                                                                                                                     >  2:30989/1‑135 (MQ=255)
 actaaCTGGGTCGCCGGGAAACATGGAATGCCACGGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAAcc                                                                                                                                    <  4:101350/132‑1 (MQ=255)
 aCGAACAGGGTCGCAGCGAACATTGTAAAGCCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAAcc                                                                                                                                    <  8:280249/135‑1 (MQ=255)
             gCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAg                                                                                                                        >  7:197503/1‑135 (MQ=255)
                 cGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACAt                                                                                                                    >  1:135179/1‑135 (MQ=255)
                      cATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTg                                                                                                               >  5:140785/1‑135 (MQ=255)
                       aTGGAACGTCACCGACGCCCCTGCCGACTACATCGAGGGCTTGCTCCGCGCCATCGGCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTgg                                                                                                              <  3:55276/135‑1 (MQ=255)
                            aaGTCACCCACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTgc                                                                                                         <  2:60737/135‑1 (MQ=255)
                              gTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTgcgc                                                                                                       >  5:159006/1‑135 (MQ=255)
                                  ccGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAg                                                                                                   >  4:186290/1‑135 (MQ=255)
                                         cccTGCCGCCTACATCGTGGGCATGCTCCTCGCCATGGTTGGTTGGCTTCTGGCCATCGCGCGCCTGCAAGGTGGGCGCAAGCTCATCCAGAACCGCTCCGGGCAAGACATCGCGGGTGGGCGCGAAGGCCCGac                                                                                            <  8:29211/135‑2 (MQ=255)
                                            tGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAgcg                                                                                          >  4:263781/1‑134 (MQ=255)
                                                      aTCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCCAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGAt                                                                               >  3:276203/1‑135 (MQ=255)
                                                      aTCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGAt                                                                               >  3:12053/1‑135 (MQ=255)
                                                           gggCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACt                                                                          >  2:147113/1‑135 (MQ=255)
                                                                tGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACa                                                                     >  8:247879/1‑135 (MQ=255)
                                                                   tCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGCGCGCCAGCCCCGATCAACTGGACAAt                                                                   >  1:49288/1‑134 (MQ=255)
                                                                   tCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGACCCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCCGCCCCGATCAACTGGACAAt                                                                   >  6:118419/1‑134 (MQ=255)
                                                                   tCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAAt                                                                   >  1:80014/1‑134 (MQ=255)
                                                                   tCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATc                                                                  >  1:174970/1‑135 (MQ=255)
                                                                    ccGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCa                                                                 >  4:187939/1‑135 (MQ=255)
                                                                               tcgGTATTCGTCAGCCCAGCGCGCGCCTCCACGGGGCGCGCAAGCTCCGCCAGAACCGCGCCGGGCAAGACATCGCTGGGGTCCGCCACGCCCTGAGCCCCAGCCCCGATCAACTGGACAATCAACTCGCGgcgc                                                      <  8:226588/135‑1 (MQ=255)
                                                                                cgGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGgcgct                                                     >  5:183211/1‑134 (MQ=255)
                                                                                cgGTATCCGTCTGCCCATCGAGCGCCTTCAAGGTGCTCGCAAGCTCAGCCATCACCGCTCAGGGCAAGACCTCGCTTGTGTGCGCGAAGGCCTGAGAGCCAGCCCCGATCAACTCGACAATCAACACGCGgcgc                                                      >  4:217355/1‑134 (MQ=255)
                                                                                            gCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAAct                                         >  5:289454/1‑135 (MQ=255)
                                                                                             ccgcgcgcgcgcgGGCACGGGGCGCGACAGCGCACCCAGACCCGCCCCGGGCAAGACATCGCGGGGGGGCGCGCCGGCCGGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAActc                                        <  1:144582/130‑1 (MQ=255)
                                                                                             ccccgcgcgccCGTGCCTGGTGCGCGCACGCCCCCCCAGACCCCCTCCCGGCAAGCCTGGGCGGGTGCGCGCGAAGGCCAGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAActc                                        <  8:91262/130‑1 (MQ=255)
                                                                                                aTCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTCCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCACCTCGCGGCGCACATGCGCCACCTCTGa                                     >  3:184934/1‑135 (MQ=255)
                                                                                                  ggcgcgcCGGCACGGGGCGCGCAAGCGCACCCAGCCCCGCGCCGGGCAAGACCGCGGTGGTGGGCGCGACGGCCGCCCCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAg                                   <  6:140785/134‑1 (MQ=255)
                                                                                                   gcgcgcCGGCCAGGTGTGCGCAAGCTCCGCCAGACCCGCCGGGGGCAACCCATCGCGTGTGGGCGCGATGGCCTGCCCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAgg                                  <  8:135206/135‑1 (MQ=255)
                                                                                                    cgcgcCGGCAAGGTGCGCGCAAGCTCAGCCAGACCCGCGCCGGGCACGACATCGCTGGTGTGCGCGCAGGCCGGCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAg                                   <  8:3167/133‑1 (MQ=255)
                                                                                                     gcgcCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCGGAAgca                                 >  6:197197/1‑132 (MQ=255)
                                                                                                       gcCTCCATGGTGCGCGCAAGCTCAGCCAGAACCGCGCCTCGCAAGACAAGGCTGGGGTTCTCGAAGGCCTTAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATa                              <  8:81777/135‑1 (MQ=255)
                                                                                                       gcCGGCCAGGTGCGCGCTAGCGCAGCCAGCCCCGCTCCGGGCAAGACAGCGGGGGGGTGCGCGCAGGCCCGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATa                              <  2:39159/135‑1 (MQ=255)
                                                                                                       gcCGGCACGGGGCGCGCAAGCTCACCCACACCCGCTCCGGGCAAGACCGCGCGGGTGTGCGCGACGGCCGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATa                              <  8:218196/135‑1 (MQ=255)
                                                                                                       gcCGGCAAGGGCCGCGCAACCGCCCCCCGCCCCGCGCCGCGCAAGCCAGCGCGGGTGCGCGCGCAGCCCGCACCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATa                              <  6:125478/135‑1 (MQ=255)
                                                                                                             tcggcgcgcccACGCTCAGCCAGCACCGCGCCGGGCAAGACCTGGCGGGTGTGCGCGAAGGCCGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  1:253906/130‑1 (MQ=255)
                                                                                                             caGGTGCGCACATGCTGCGCAAGACCCGCTCGTGCCCAGACTTCGCTGGTGTTCGCTAAGCCGGGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  8:64431/134‑1 (MQ=255)
                                                                                                             acggggcgctcAAGCCCAGCCACACCCGCCCCGGGCAAGACATCGCGGGTGGGCGCGAAGGCCTGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  4:166936/130‑1 (MQ=255)
                                                                                                             acggggcgcgcACGCTCCGCCAGCCCCGCTCCGGGCACGACAGCGCTGGGGGGCGCGACGGCCGGCCCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  3:88218/130‑1 (MQ=255)
                                                                                                             acggggcgcgcACGCGCAGCCAGACCCGCGCCGGGCAAGACAGCGCGGGGGGGCGCGACGGCCGCCCCCCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  4:12053/130‑1 (MQ=255)
                                                                                                             aaGGGGCGCGCAAGCGCCTCCCGACCCGCTCGGGGCAAGAGTTCGCTGGGGTGCGCGAAGGCCAGCGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCAt                        <  3:64039/135‑1 (MQ=255)
                                                                                                                   cgcgcAAGCTCCGCCAGAACCGCGCCGGGCACGCCCGCGGGGGGGTGCGCGACGCCCGGACCCCCAGCCCCGATCAACTAGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCaa                  <  5:161196/135‑1 (MQ=255)
                                                                                                                    gcgcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACCTGCGCCACCTCTGAAGGAATAGCCCATGCCCAAc                 >  5:186183/1‑135 (MQ=255)
                                                                                                                    gcgcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAc                 >  8:2361/1‑135 (MQ=255)
                                                                                                                    gcgcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGACCAAc                 >  7:303201/1‑135 (MQ=255)
                                                                                                                      gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAcc                >  2:166203/1‑134 (MQ=255)
                                                                                                                      gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAcc                >  6:155799/1‑134 (MQ=255)
                                                                                                                      gcAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAAcc                >  1:185567/1‑134 (MQ=255)
                                                                                                                              aGCCAGCACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGACGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGAACATGCGCCAACTCTGAAGGAATAGCCCATGCCCCACCTTCCACTgc       >  2:203492/1‑135 (MQ=255)
                                                                                                                              aGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTgc       >  4:125894/1‑135 (MQ=255)
                                                                                                                                   gAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCACGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCt  >  1:303596/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CACGAACAGGGTCGCAGCGAACCATGGAAAGTCACCGACGCCCCTGCCGACTACATCGAGGGCATGCTCCGCGCCATCGTCGGTATCCGTCTGCCCATCGCGCGCCTGCAAGGTGCGCGCAAGCTCAGCCAGAACCGCTCCGGGCAAGACATCGCTGGTGTGCGCGAAGGCCTGAGCGCCAGCCCCGATCAACTGGACAATCAACTCGCGGCGCACATGCGCCAACTCTGAAGGAATAGCCCATGCCCAACCTTACCCTGCGCCCT  >  ACB122/6084631‑6084896

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