Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,077,829 G→A Q87Q (CAG→CAA)  cbtA → cytoskeleton‑binding toxin CbtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,077,8290GA97.1% 122.9 / NA 35Q87Q (CAG→CAA) cbtAcytoskeleton‑binding toxin CbtA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/17);  minor base C (1/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACGGTTATTCACTCCG  >  NC_000913/2077635‑2078067
                                                                                                                                                                                                  |                                                                                                                                                                                                                                              
cAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGAcccc                                                                                                                                                                                        >  2:82737/1‑250 (MQ=255)
        gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGc                                                                                                                                                                                                   <  1:207582/232‑1 (MQ=255)
        gctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGc                                                                                                                                                                                                   >  2:207582/1‑232 (MQ=255)
                                         cacaCTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGAGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTaa                                                                                                                                                 <  2:69217/249‑1 (MQ=255)
                                            acTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCATCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATc                                                                                                                                            >  2:221770/1‑251 (MQ=255)
                                                 aTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGCTTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACCTTACCCTGGGTAAGTATCCgg                                                                                                                                         >  1:355412/1‑249 (MQ=255)
                                                    acacacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTTTGGTATGCGGTGAACTTTCTCGTGGAAAAATACGCCCTGGTGCGTACCGACCAGCCGGGATTCATCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGTCTCGCAGGGCGATCGGCCGGATGACCCGCGACAATTACAGGACGGTAAATAACATGTCCCTGTGTaattagtcggatgcg                                                                                                                                    <  2:207860/251‑14 (MQ=255)
                                                             ttGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGAGATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGGGGCGAAatgatg                                                                                                                            <  1:132662/250‑1 (MQ=255)
                                                                ccGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCCTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCAGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGaca                                                                                                                         >  1:42250/1‑250 (MQ=255)
                                                                           gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCACTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGCAGCCGAc                                                                                                             >  2:271894/1‑251 (MQ=255)
                                                                                   gcagcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTaa                                                                                                     <  2:153890/251‑1 (MQ=255)
                                                                                     agcaTATTGAGGCAGGCCTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCCGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGAGACAATTAAAGAAAGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGACGACAGCGTTAACg                                                                                                   >  2:302001/1‑251 (MQ=255)
                                                                                             gAGGCAGGCATTTAACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGAATCCGTAGGCGAAATGATGACACTGGAAGCCGACAGCGTTACCGTACAGgct                                                                                           >  2:115561/1‑250 (MQ=255)
                                                                                              aggcaggcATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGgcgc                                                                                          >  1:325047/1‑251 (MQ=255)
                                                                                                       tttCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGata                                                                                    <  2:247640/248‑1 (MQ=255)
                                                                                                       tttCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGata                                                                                    >  1:247640/1‑248 (MQ=255)
                                                                                                           aCTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGTCCAGCCGGGATTCAGCGCCTGTACCCGCCCTCAATTAATAAACAGCATCGATATCCTCCggg                                                                                                                                                                                                                    <  2:186983/116‑1 (MQ=255)
                                                                                                           aCTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCggg                                                                                                                                                                                                                    >  1:186983/1‑116 (MQ=255)
                                                                                                           aCTGTGTGAGGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGca                                                                             <  1:348118/251‑1 (MQ=255)
                                                                                                                gtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACAtt                                                                         <  2:193033/250‑1 (MQ=255)
                                                                                                                gtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTg                                                                        <  1:333255/251‑1 (MQ=255)
                                                                                                                              tttCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGAt                                                          >  1:179776/1‑251 (MQ=255)
                                                                                                                                   cGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGgtgt                                                     <  2:103614/251‑1 (MQ=255)
                                                                                                                                   cGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGgtgt                                                     <  2:408790/251‑1 (MQ=255)
                                                                                                                                    gTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGGAAATAACATGTCCCTGGGTAAGTATCCGGAGGCTAAATGATGACACTGGAAGCCTACAGCGTTAACGGACAGGCGTTGGATATTGGGCACATTGTCGTTGACATTGATGGTGtt                                                    <  1:259906/251‑1 (MQ=255)
                                                                                                                                         aaaaTACGCCCTGCGGCGTACCGCCCAGCCGGGCTTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCATCCGGGCTCGCAGGGCGACCGGCCTAATGACCCGCGACAATTACAGACCGGTAAATACCATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGACGACAGCGTTAACGTACAGGCGCTGGATATGGGGAACATTGTCGTTGACATTGATGGTGTTAatat                                               >  2:60687/1‑251 (MQ=255)
                                                                                                                                           aaTACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCa                                             <  1:242006/251‑1 (MQ=255)
                                                                                                                                                     ggTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCactgaactg                                     <  2:59322/249‑1 (MQ=255)
                                                                                                                                                          gTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACa                                                                                                                                                                                                                                      <  1:249152/51‑1 (MQ=255)
                                                                                                                                                          gTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAATTAATAAACa                                                                                                                                                                                                                                      >  2:249152/1‑51 (MQ=255)
                                                                                                                                                                              aGCGCCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACGGTTAt          >  1:315746/1‑251 (MQ=255)
                                                                                                                                                                                cgcCTGTACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGACCGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACGGTTATTc        >  2:142007/1‑251 (MQ=255)
                                                                                                                                                                                      tACCCGCTCTCAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACGGTTATTCACTCcg  >  1:380641/1‑251 (MQ=255)
                                                                                                                                                                                             tctcAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACgg              >  2:256882/1‑232 (MQ=255)
                                                                                                                                                                                             tctcAATTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACgg              <  1:256882/232‑1 (MQ=255)
                                                                                                                                                                                                  |                                                                                                                                                                                                                                              
CAGATACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGCGAAATGATGACACTGGAAGCCGACAGCGTTAACGTACAGGCGCTGGATATGGGGCACATTGTCGTTGACATTGATGGTGTTAATATCACTGAACTGATTAATAAGGCCGCTGAAAACGGTTATTCACTCCG  >  NC_000913/2077635‑2078067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑