Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 901,420 T→C 8.9% intergenic (‑139/‑192) 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,4200TC 8.9% 153.9 / 4.9 56intergenic (‑139/‑192)cysZ/trxBsulfate transport protein CysZ/thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base C (2/3);  total (26/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 = 9.22e-01

GGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCCGG  >  NC_002947/901273‑901568
                                                                                                                                                   |                                                                                                                                                    
ggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATg                                                                                                                                                     >  8:88926/1‑149 (MQ=255)
      cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTAt                                                                                                                                                <  3:70882/148‑1 (MQ=255)
        tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAg                                                                                                                                             >  5:182953/1‑149 (MQ=255)
        tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAg                                                                                                                                             >  7:151896/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                                                                                                                                      <  1:146011/136‑1 (MQ=255)
                            ggAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                      >  2:146011/1‑136 (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                                                   <  5:15364/148‑1 (MQ=255)
                                                                                                   aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCtt                                                   <  7:154107/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                                                 <  4:67929/111‑1 (MQ=255)
                                                                                                                                      acctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcc                                                 >  3:67929/5‑115 (MQ=255)
                                                                                                                                               gagATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTATGCGACATCGCTACAGTGCGCTCAACTTTTCGCTTTGAGGGCCTGCGCGTTTTAGACTACCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGTCGCaa      >  5:117551/3‑148 (MQ=255)
                                                                                                                                               gCGATGACTATGAGGTGGTTTGTCTAAGGTTATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGgcgc        <  6:117551/147‑1 (MQ=255)
                                                                                                                                                 gATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc    >  4:77490/1‑149 (MQ=255)
                                                                                                                                                   tGCCTATGAGTTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTAGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCCgg  <  3:229460/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCCGG  >  NC_002947/901273‑901568

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