Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,861 C→A 7.1% R310R (CGG→CGT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,8610CA 7.1% 196.0 / 2.9 56R310R (CGG→CGT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base A (2/2);  total (30/26)
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.57e-01

ACCATGGTTGAAAAAAGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGAT  >  NC_000913/1579724‑1579994
                                                                                                                                         |                                                                                                                                      
aCCATGGTTGAAAAAAGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCg                                                                                                                                      >  1:1247464/1‑140 (MQ=255)
aCCATGGTTGAAAAAAGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCg                                                                                                                                      >  1:1247491/1‑140 (MQ=255)
           aaaaaGGATTTTATAACCTTCGCTAAATTTAGATACAGTCTCATTGTTTAATTTAGTAATTCGTTGATTAGGACAACGGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTgcg                                                                                                                          <  1:1673078/140‑1 (MQ=255)
            aaaaGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTgcgc                                                                                                                         >  2:559561/1‑140 (MQ=255)
               aGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:1093942/140‑1 (MQ=255)
               aGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:685856/140‑1 (MQ=255)
               aGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:1265581/140‑1 (MQ=255)
               aGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:1532207/140‑1 (MQ=255)
               aGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:150428/140‑1 (MQ=255)
                ggATTTTATAACCTTCGCAAAATTTTTTAAGAGTCTCATTGTTTACTTTAGTACTACGATGCTTAGGTGAACCGCCATTGCAGTTAGGTTTTTTTGCACATTGCTGACATTTCGCTGGATTCCGTTTTTTTTTGCGCTGt                                                                                                                      <  1:163478/140‑1 (MQ=255)
                ggATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTTGCGCTGt                                                                                                                      <  1:163432/140‑1 (MQ=255)
                ggATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTc                                                                                                                     >  2:978903/1‑140 (MQ=255)
                ggATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTc                                                                                                                     >  2:803828/1‑140 (MQ=255)
                 gATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCa                                                                                                                    >  2:109780/1‑140 (MQ=255)
                 gATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCa                                                                                                                    >  2:260637/1‑140 (MQ=255)
                         aaCCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTaca                                                                                                            <  1:1446628/140‑1 (MQ=255)
                         aaCCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTaca                                                                                                            <  1:27875/140‑1 (MQ=255)
                           ccTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTATTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACt                                                                                                          >  2:528117/1‑140 (MQ=255)
                           ccTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACt                                                                                                          <  2:1114546/140‑1 (MQ=255)
                           ccTTCGCAAAAATAAGAAACAGTCTCATTGTTTAATTTAGTAATACGATGCTTAGGAAAACAGCCATTGCAGATAAGTTAATTTGCACATTTCTTACATTTCTCTTGTATCCG‑TTTTATTTGCGAATACAGTTGTAAACt                                                                                                          >  2:838624/1‑140 (MQ=255)
                               cGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTc                                                                                                      >  1:1680842/1‑140 (MQ=255)
                                      aTAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGttt                                                                                               >  2:815225/1‑140 (MQ=255)
                                               aGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAg                                                                                      <  1:1317437/140‑1 (MQ=255)
                                               aGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAg                                                                                      <  1:109780/140‑1 (MQ=255)
                                                 tctcATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGAt                                                                                    >  1:1267036/1‑140 (MQ=255)
                                                              tttAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTc                                                                       <  2:1115220/140‑1 (MQ=255)
                                                              tttAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTc                                                                       <  1:1013101/140‑1 (MQ=255)
                                                                 aGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTTTACACAGTTCATCGTTTTGATTTAAGATTTATTAAAGTTTACaa                                                                    >  1:404636/1‑140 (MQ=255)
                                                                 aGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCaa                                                                    >  1:404154/1‑140 (MQ=255)
                                                                       aCGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAAtgttgt                                                              >  2:966539/1‑140 (MQ=255)
                                                                       aCGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGt                                                              >  2:872676/1‑140 (MQ=255)
                                                                         gATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATTCAAATTGCTGACATTTCGCTGGAATCCT‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTAc                                                            <  2:880686/140‑1 (MQ=255)
                                                                         gATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTAc                                                            <  1:247484/140‑1 (MQ=255)
                                                                           tGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACtg                                                          >  1:263118/1‑140 (MQ=255)
                                                                                aGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGAt                                                     >  2:1395282/1‑140 (MQ=255)
                                                                                aGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTACATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGAt                                                     >  2:1395312/1‑140 (MQ=255)
                                                                                       ccGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATTTTTCCAATTTTGTACTGTGGATAGACaaa                                              >  2:856332/1‑140 (MQ=255)
                                                                                       ccGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACaaa                                              >  1:667743/1‑140 (MQ=255)
                                                                                               gCAGTTAGGTTTATATGCACTTTGATGACATTTGGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCTTCGTTTTGTGTTCAGATTTATTTATGTTTCCTATTTTGTACTGTGGTTAGACAAAATGGTCGc                                      <  1:1308807/140‑1 (MQ=255)
                                                                                                 aGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATACG‑TTTTTTTTGCGCTGt                                                                                                                      >  2:1100683/1‑58 (MQ=255)
                                                                                                 aGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATACG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:1100665/58‑1 (MQ=255)
                                                                                                 aGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATACG‑TTTTTTTTGCGCTGt                                                                                                                      >  2:1100665/1‑58 (MQ=255)
                                                                                                 aGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATACG‑TTTTTTTTGCGCTGt                                                                                                                      <  1:1100683/58‑1 (MQ=255)
                                                                                                   aTAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTc                                  <  2:1247491/140‑1 (MQ=255)
                                                                                                   aTAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTc                                  <  2:1247464/140‑1 (MQ=255)
                                                                                                          tataTGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATg                           >  1:958889/1‑140 (MQ=255)
                                                                                                          tataTGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATg                           >  1:958033/1‑140 (MQ=255)
                                                                                                                caccTTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGTGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCa                     <  1:1364464/140‑1 (MQ=255)
                                                                                                                 acaTTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCAt                    >  2:914965/1‑140 (MQ=255)
                                                                                                                 acaTTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTCATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCAt                    >  2:303134/1‑140 (MQ=255)
                                                                                                                     tGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACt                >  1:746320/1‑140 (MQ=255)
                                                                                                                             tttCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTcaacaa        >  2:621773/1‑140 (MQ=255)
                                                                                                                               tCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAAcc      >  1:963860/1‑140 (MQ=255)
                                                                                                                                   tGGTATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGAt  <  2:1404861/140‑1 (MQ=255)
                                                                                                                                   tGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGAt  <  1:1438874/140‑1 (MQ=255)
                                                                                                                                   tGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGAt  <  1:1484545/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
ACCATGGTTGAAAAAAGGATTTTATAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTACTTTAGTAATACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCG‑TTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGAT  >  NC_000913/1579724‑1579994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑