Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,782,079 C→A 8.0% G96V (GGC→GTC)  iscS‑II ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,782,0790CA 8.0% 156.0 / 4.4 50G96V (GGC→GTC) iscS‑IIcysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGG  >  NC_002947/2781936‑2782222
                                                                                                                                                |                                                                                                                                               
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                                                                               >  6:201868/1‑134 (MQ=255)
                                                                            cAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCt                                                                               <  5:201868/134‑1 (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                                                        <  2:262217/149‑1 (MQ=255)
                                                                                    cGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTgg                                                        >  4:135822/1‑149 (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                                                >  5:274506/1‑141 (MQ=255)
                                                                                                    ttGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGc                                                <  6:274506/141‑1 (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)
                                                                                                                                           gcgcgCCCGAGGGCCATCCCAGCTGCGCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAc       <  6:350086/143‑1 (MQ=255)
                                                                                                                                             gcgcCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCcagg  >  4:46616/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGGTGACGGCCCAGCCCTGACGTTCCAGCTCGGCCACGGTGTCGAGCACGGCCTTGTGTTCCAGCTGGCTGGTGATCAGGTGTCCGGGCTGGCCGATGCCTTGGGCAATACCCTTGAGCGCCAGGTTGTTGGACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGG  >  NC_002947/2781936‑2782222

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: