Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,821 C→G 100% C104S (TGC→TCC)  iscR ← DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,8210CG100.0% 47.2 / NA 16C104S (TGC→TCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (7/9);  total (7/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGT  >  NC_000913/2661685‑2661958
                                                                                                                                        |                                                                                                                                         
gcgttgtgcgttgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTAt                                                                                                                                     <  1:169000/143‑1 (MQ=255)
             tgtgcaTGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCa                                                                                                                        <  2:3245/138‑1 (MQ=255)
                         cgaaCAGACACATCCATCACTTCCTGGTTATTAACCAGTTCGCCTAAATTAATGTTGTTGAGAAAACCTGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                                            <  1:164693/143‑1 (MQ=255)
                             gaGACACAGCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                                                                                                        <  1:261171/142‑1 (MQ=255)
                             cAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                                                                                                        <  2:286853/143‑1 (MQ=255)
                                   caTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACgg                                                                                                  <  1:187288/143‑1 (MQ=255)
                                                       ttaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAGACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                                              <  2:84987/143‑1 (MQ=255)
                                                          aCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCGTGGCAGCCGCCTTTAACCTGACAACGGGAGGCATCTAGAGGTTCGTCAAc                                                                           >  2:212021/1‑143 (MQ=255)
                                                              gTTCGCCCAACGTAATGTTGTTGCGAAACCCGGTGAGACGGTCGCTCACATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGgcg                                                                       <  2:234780/143‑1 (MQ=255)
                                                                              gttgttGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCACCGGCGCTAACTACTTCACCaa                                                       >  1:43443/1‑143 (MQ=255)
                                                                                    gagaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGTAGCCGCCTTTACCCTGACAACGGGTGACATCTACAGATTCGTCAACAGCGCTAAGTACTTCGCCAACGGCGa                                                 >  1:81687/1‑143 (MQ=255)
                                                                                                         gCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTa                            >  1:160496/1‑143 (MQ=255)
                                                                                                           tgacgtCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc                          <  2:177262/138‑1 (MQ=255)
                                                                                                                 cACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCCTCTACCGATTCATCAACTCCGCTAATTACTTCGCCAACAGCGACGCTGCTGGCAACATTGCCTAACAGATaa                    >  1:221365/1‑143 (MQ=255)
                                                                                                                       acaGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTACCAGCTaccacccg             >  2:252886/1‑141 (MQ=255)
                                                                                                                                   tCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTACTTACTTCGCCAACGGCGCTGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACga  >  1:12728/1‑142 (MQ=255)
                                                                                                                                        |                                                                                                                                         
GCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGT  >  NC_000913/2661685‑2661958

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