Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 901,418 G→A 9.2% intergenic (‑137/‑194) cysZ ← / → trxB sulfate transport protein CysZ/thioredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947901,4180GA 9.2% 147.8 / 4.5 54intergenic (‑137/‑194)cysZ/trxBsulfate transport protein CysZ/thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base A (2/3);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901273‑901566
                                                                                                                                                 |                                                                                                                                                    
ggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATg                                                                                                                                                   >  8:88926/1‑149 (MQ=255)
      cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTAt                                                                                                                                              <  3:70882/148‑1 (MQ=255)
        tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAg                                                                                                                                           >  7:151896/1‑149 (MQ=255)
        tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAg                                                                                                                                           >  5:182953/1‑149 (MQ=255)
         cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAgg                                                                                                                                          >  6:71910/1‑149 (MQ=255)
           agagACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTTTGAgggg                                                                                                                                        >  4:229460/1‑147 (MQ=255)
            gagaCTCTGCAGAGAAGGAAAACGCGCTTACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTTGGCGGGAAACGGGATGTAACAAAAGCGCCCGCGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgctg                                                                                                                                          <  1:252475/146‑4 (MQ=255)
            gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         <  5:355940/147‑1 (MQ=255)
            gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTgg                                                                                                                                       <  2:139883/149‑1 (MQ=255)
             agaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                      <  7:88926/149‑1 (MQ=255)
                 tCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGTCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCCCAGCATAGGCGATGCCTATGAGGTGGATTGt                                                                                                                                  >  8:82404/1‑149 (MQ=255)
                     cAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              <  5:71910/149‑1 (MQ=255)
                          aaGGAACACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCGTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAGTGGCTACGGata                                                                                                                         >  8:294910/1‑149 (MQ=255)
                            ggAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    >  2:146011/1‑136 (MQ=255)
                            ggAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    <  1:146011/136‑1 (MQ=255)
                               aaaCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTcc                                                                                                                    >  1:321865/1‑149 (MQ=255)
                                aaCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCt                                                                                                                   <  2:305473/149‑1 (MQ=255)
                                         aCCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGTGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           <  7:79197/148‑1 (MQ=255)
                                         aCCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTg                                                                                                          >  7:212378/1‑149 (MQ=255)
                                          ccTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc                                                                                                         >  8:217804/1‑149 (MQ=255)
                                                    tttGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTa                                                                                               <  7:74992/149‑1 (MQ=255)
                                                         cgcgCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTg                                                                                          >  7:36632/1‑149 (MQ=255)
                                                            gcCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTgcgc                                                                                       >  3:158666/1‑149 (MQ=255)
                                                             ccGACGAAAGGCACGCCGTGCGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCt                                                                                      <  6:240572/149‑1 (MQ=255)
                                                                  gaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt                                                                                 >  4:185445/1‑149 (MQ=255)
                                                                  gaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgct                                                                                                                                           >  7:69446/1‑89 (MQ=255)
                                                                  gaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgct                                                                                                                                           <  8:69446/91‑3 (MQ=255)
                                                                     aGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACtttt                                                                              <  4:158666/149‑1 (MQ=255)
                                                                     aGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACtttt                                                                              <  6:182953/149‑1 (MQ=255)
                                                                       gCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCg                                                                            >  6:202233/1‑149 (MQ=255)
                                                                        cACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGc                                                                           >  8:154107/1‑149 (MQ=255)
                                                                         aCGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCt                                                                          >  8:12782/1‑149 (MQ=255)
                                                                              gTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGc                                                                      <  7:82404/148‑1 (MQ=255)
                                                                                   cGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCt                                                                <  7:217804/149‑1 (MQ=255)
                                                                                     ggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgc                                                              <  8:36632/149‑1 (MQ=255)
                                                                                            ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCt                                                                                                                   >  1:220163/1‑89 (MQ=255)
                                                                                            ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCt                                                                                                                   <  2:220163/89‑1 (MQ=255)
                                                                                            ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGtttt                                                       <  8:212378/149‑1 (MQ=255)
                                                                                             gATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  7:12782/149‑1 (MQ=255)
                                                                                                   aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCtt                                                 <  7:154107/148‑1 (MQ=255)
                                                                                                   aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCtt                                                 <  5:15364/148‑1 (MQ=255)
                                                                                                      aaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcccc                                             <  5:202233/149‑1 (MQ=255)
                                                                                                        aGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGGCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           <  2:74566/149‑1 (MQ=255)
                                                                                                         gcgcCCGGGGCTTTTTTATCGATGAAATCGCAGGATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGGGGGCCTGCGCGTTTTAGCCTTCCCCAAg                                          <  5:273997/149‑1 (MQ=255)
                                                                                                         gcgcCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAg                                          <  3:185445/149‑1 (MQ=255)
                                                                                                              cGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGtgctg                                     >  4:211048/1‑149 (MQ=255)
                                                                                                                  gCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTg                                 >  5:186559/1‑149 (MQ=255)
                                                                                                                       tttATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTc                             <  3:211048/148‑1 (MQ=255)
                                                                                                                          aTCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      >  4:54138/1‑120 (MQ=255)
                                                                                                                          aTCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  3:54138/120‑1 (MQ=255)
                                                                                                                                  aaTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTc                 >  4:351743/1‑149 (MQ=255)
                                                                                                                                      acctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcc                                               >  3:67929/5‑115 (MQ=255)
                                                                                                                                      acctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcc                                               <  4:67929/111‑1 (MQ=255)
                                                                                                                                               gagATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTATGCGACATCGCTACAGTGCGCTCAACTTTTCGCTTTGAGGGCCTGCGCGTTTTAGACTACCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGTCGCaa    >  5:117551/3‑148 (MQ=255)
                                                                                                                                               gCGATGACTATGAGGTGGTTTGTCTAAGGTTATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGgcgc      <  6:117551/147‑1 (MQ=255)
                                                                                                                                                 gATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  >  4:77490/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901273‑901566

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