Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,782,074 T→G 8.0% T98P (ACC→CCC)  iscS‑II ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,782,0740TG 8.0% 147.9 / 5.2 50T98P (ACC→CCC) iscS‑IIcysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base G (2/2);  total (25/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 = 6.56e-01

CTGCCAGGCGGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAC  >  NC_002947/2781927‑2782217
                                                                                                                                                    |                                                                                                                                               
cTGCCAGGCGGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTc                                                                                                                                                <  8:264287/149‑1 (MQ=255)
cTGCCAGGCGGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTc                                                                                                                                                <  3:78364/149‑1 (MQ=255)
         gggTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAg                                                                                                                                       <  1:94183/149‑1 (MQ=255)
                ggCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATa                                                                                                                                >  7:249509/1‑149 (MQ=255)
                                         gcgaCGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGAc                                                                                                       <  8:143093/146‑1 (MQ=255)
                                                tgtCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCg                                                                                                >  2:212775/1‑149 (MQ=255)
                                                    gAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACCTCGGGCGCGACCGa                                                                                                                                        >  6:144747/1‑106 (MQ=255)
                                                    gAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACCTCGGGCGCGACCGa                                                                                                                                        <  5:144747/106‑1 (MQ=255)
                                                            ccTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTgg                                                                                    >  4:232853/1‑149 (MQ=255)
                                                            ccTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTgg                                                                                    >  8:15536/1‑149 (MQ=255)
                                                             ctggtgtTCCAGCTGGCGGGTGATCAGGTGGCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGgc                                                                                   <  8:249509/146‑1 (MQ=255)
                                                             cTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGTGCCGACGCACTCGGc                                                                                              >  3:27614/1‑138 (MQ=255)
                                                             cTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGTGCCGACGCACTCGGc                                                                                              <  4:27614/138‑1 (MQ=255)
                                                                      cAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGc                                                                                              <  8:20426/129‑1 (MQ=255)
                                                                      cAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGc                                                                                              >  7:20426/1‑129 (MQ=255)
                                                                              tggTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCg                                                                  >  6:187178/1‑149 (MQ=255)
                                                                                  gATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGAc                                                              >  1:262217/1‑149 (MQ=255)
                                                                                  gATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGAc                                                              >  4:337016/1‑149 (MQ=255)
                                                                                     cAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCt                                                                          <  5:201868/134‑1 (MQ=255)
                                                                                     cAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCt                                                                          >  6:201868/1‑134 (MQ=255)
                                                                                       ggTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCa                                                         <  1:212775/149‑1 (MQ=255)
                                                                                        gtgtCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTg                                                                >  5:40142/1‑141 (MQ=255)
                                                                                        gtgtCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTg                                                                <  6:40142/141‑1 (MQ=255)
                                                                                         tgtCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGTTGGCGCGCCTGCTCCACGGCGTGACGTGCAGc                                                       <  7:25519/149‑1 (MQ=255)
                                                                                             cGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTgg                                                   >  4:135822/1‑149 (MQ=255)
                                                                                             cGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTgg                                                   <  2:262217/149‑1 (MQ=255)
                                                                                                gCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGcc                                                 >  6:222296/1‑148 (MQ=255)
                                                                                                    gCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGaa                                              <  7:15536/147‑1 (MQ=255)
                                                                                                     ccGATGCCTTGGGCAATGCCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGc                                           <  1:262821/149‑1 (MQ=255)
                                                                                                     ccGATGCCTTGGGCAATGCCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGc                                           >  2:262821/1‑149 (MQ=255)
                                                                                                             ttGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCTACGCACTCGGCGACTTGc                                                                                       >  5:250439/1‑97 (MQ=255)
                                                                                                             ttGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGc                                                                                       <  6:250439/97‑1 (MQ=255)
                                                                                                             ttGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGc                                           <  6:274506/141‑1 (MQ=255)
                                                                                                             ttGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGc                                           >  5:274506/1‑141 (MQ=255)
                                                                                                               gggCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGc                                   <  5:222296/147‑1 (MQ=255)
                                                                                                               gggCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTg                                 >  8:213432/1‑149 (MQ=255)
                                                                                                                   aaTACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCg                                                                                                <  8:307950/82‑1 (MQ=255)
                                                                                                                        ccTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCggg                        >  1:46530/1‑149 (MQ=255)
                                                                                                                          ttGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTg                                                    <  8:18289/119‑1 (MQ=255)
                                                                                                                          ttGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGtt                      >  1:189079/1‑149 (MQ=255)
                                                                                                                          gtGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGGCGAAGGCGTGGCCGCTGGCGGCCGGGtt                      >  5:284072/2‑149 (MQ=255)
                                                                                                                               cgcCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCg                  >  4:289843/1‑148 (MQ=255)
                                                                                                                                 ccAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAg               <  3:337016/149‑1 (MQ=255)
                                                                                                                                   aGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGtt             <  2:46530/149‑1 (MQ=255)
                                                                                                                                         ttgGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTg       <  3:289843/149‑1 (MQ=255)
                                                                                                                                          tgGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTgg      >  2:305526/1‑149 (MQ=255)
                                                                                                                                          tgGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTg       <  5:5445/148‑1 (MQ=255)
                                                                                                                                          tgGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTg       >  6:5445/1‑148 (MQ=255)
                                                                                                                                            gaccTCGGGCGCGACCGAGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTg                                                    >  2:139428/4‑101 (MQ=255)
                                                                                                                                            gaccTCGGGCGCGACCGAGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTg                                                    <  1:139428/98‑1 (MQ=255)
                                                                                                                                                    tCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGCCGCCGGGTTGCCGAAGTTGGCCTGGTTAc  >  5:350086/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTGCCAGGCGGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAC  >  NC_002947/2781927‑2782217

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