Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,140,564 A→C 7.6% Q373P (CAG→CCG)  sfnB → SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,140,5640AC 7.6% 161.1 / 4.6 53Q373P (CAG→CCG) sfnBSfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (26/23);  new base C (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

GCGCCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACGCCA  >  NC_002947/3140420‑3140701
                                                                                                                                                |                                                                                                                                         
gcgcCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAggg                                                                                                                                        <  5:59094/148‑1 (MQ=255)
 cgcCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGc                                                                                                                                      <  1:37363/149‑1 (MQ=255)
  gcCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCt                                                                                                                                     >  2:320512/1‑149 (MQ=255)
  gcCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGGGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCt                                                                                                                                     <  3:296150/149‑1 (MQ=255)
     gTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGa                                                                                                                                  <  2:199102/149‑1 (MQ=255)
      tCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGAt                                                                                                                                 >  1:12584/1‑149 (MQ=255)
      tCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGAt                                                                                                                                 >  1:18380/1‑149 (MQ=255)
          cGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCt                                                                                                                             <  4:134391/149‑1 (MQ=255)
            aGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTc                                                                                                                           <  2:79448/149‑1 (MQ=255)
               cTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCCctggctggc                                                                                                                                  <  6:303861/139‑2 (MQ=255)
               cTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCCctggctggc                                                                                                                                  >  5:303861/1‑138 (MQ=255)
               cTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTAGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCCGGATCGCTTCTgg                                                                                                                        <  2:167000/149‑1 (MQ=255)
                               tGCTGACGCCGCAGTGGTGGGCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCAccc                                                                                                        <  8:68749/149‑1 (MQ=255)
                                gCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCca                                                                                                       >  5:62272/1‑149 (MQ=255)
                                        cTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCCctggctggc                                                                                                                                  >  2:145954/1‑113 (MQ=255)
                                        cTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCCctggctggc                                                                                                                                  <  1:145954/114‑2 (MQ=255)
                                         tCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCa                                                                                              >  4:13013/1‑149 (MQ=255)
                                         tCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGAGTCCAGTCTCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGACAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCAAACGCTGCa                                                                                              >  4:18046/1‑149 (MQ=255)
                                             tggtggTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGAc                                                                                          <  4:61490/149‑1 (MQ=255)
                                               gtggtCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGAccc                                                                                        >  4:113093/1‑149 (MQ=255)
                                               gtggtCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGATGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGATGGATCGCTTCTGGCGCAATTCGCGCACCCACACCCTGCATGaa                                                                                          >  6:247474/1‑146 (MQ=255)
                                                tggtCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCg                                                                                       <  6:62272/149‑1 (MQ=255)
                                                         aggccaAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCGCACGCTGCATGACCCGGTGCGCTg                                                                               >  8:3511/1‑148 (MQ=255)
                                                           gccaAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGaa                                                                            >  3:25021/1‑149 (MQ=255)
                                                                      gTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTgg                                                                 >  2:270618/1‑149 (MQ=255)
                                                                       tCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTgg                                                                 <  1:144826/148‑1 (MQ=255)
                                                                         tACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCg                                                              >  7:350253/1‑149 (MQ=255)
                                                                           ccGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTg                                                                               <  4:135396/130‑1 (MQ=255)
                                                                           ccGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTg                                                                               >  3:135396/1‑130 (MQ=255)
                                                                           ccGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCgg                                                             <  6:137673/148‑1 (MQ=255)
                                                                           ccGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCgg                                                             >  5:265501/1‑148 (MQ=255)
                                                                           ccGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCgg                                                             <  3:113093/148‑1 (MQ=255)
                                                                                   tCCCTAGAGGCGTCCAGTCGCGTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTac                                                                        <  1:97841/129‑1 (MQ=255)
                                                                                   tCCCTAGAGGCGTCCAGTCGCGTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTac                                                                        >  2:97841/1‑129 (MQ=255)
                                                                                       tAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCg                                                              >  6:48308/1‑135 (MQ=255)
                                                                                       tAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCg                                                              <  5:48308/135‑1 (MQ=255)
                                                                                        agagGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTc                                                >  8:46923/1‑148 (MQ=255)
                                                                                                  aGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCaa                                     <  8:350253/149‑1 (MQ=255)
                                                                                                   gTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCCCCGGGGGGGCAAGCACGGGGTCAACgggagcagg                                    >  7:235595/1‑142 (MQ=255)
                                                                                                           ctgGAGCTGGGCGTCACGCGCAGCGTTTCGGCCAGCAAGGGTCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACg                            <  3:18046/149‑1 (MQ=255)
                                                                                                           ctgGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACg                            <  3:13013/149‑1 (MQ=255)
                                                                                                            tgGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACg                            >  4:185762/1‑148 (MQ=255)
                                                                                                            tgGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGt                           >  1:185178/1‑149 (MQ=255)
                                                                                                                  tGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTcc                     >  1:298741/1‑149 (MQ=255)
                                                                                                                   gggCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGAACCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCt                    >  2:300151/1‑149 (MQ=255)
                                                                                                                    ggcggcACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTg                   <  1:320512/149‑1 (MQ=255)
                                                                                                                         cACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCGCACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACt              <  7:3511/149‑1 (MQ=255)
                                                                                                                           cgcgCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGa            <  3:185762/149‑1 (MQ=255)
                                                                                                                              gcAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGt         >  8:90264/1‑149 (MQ=255)
                                                                                                                                  cGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACg     <  7:46923/149‑1 (MQ=255)
                                                                                                                                    tttCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACGc    >  6:326207/1‑148 (MQ=255)
                                                                                                                                    tttCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACGc    >  6:247239/1‑148 (MQ=255)
                                                                                                                                    tttCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCT‑GAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACGCCa  >  6:338268/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                         
GCGCCGTCGCCGAGCCTAGCGAAGAGCACGTTGCTGAAGCCTCAGTGGTGGTCGGCCAGGCCAAGGTGCTGTCTACCGAGATTTCCCTAGAGGCGTCCAGTCGCCTGCTGGAGCTGGGCGGCACGCGCAGCGTTTCGGCCAGCCAGGGGCTGGATCGCTTCTGGCGCAATGCGCGCACCCACACGCTGCATGACCCGGTGCGCTGGAAGTACCACCTGGTCGGCAATCAACTGCTCAACGGCATCAAGCCGCAACGTCATTCCTGGAACTGAGGTCACGCCA  >  NC_002947/3140420‑3140701

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