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,8210CG98.4% 59.2 / ‑6.7 21C104S (TGC→TCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (10/10);  total (10/11)
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.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGT  >  NC_000913/2661690‑2661958
                                                                                                                                   |                                                                                                                                         
gcgcGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCGCCGACCACGCCACAGCGCGTGGGTCAGGGATTTATcgccg                                                                                                                                <  2:313819/143‑1 (MQ=255)
 ctcgtcgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTCGTTATTAACCAGTTCGCCTAATGTAATGTTGTTCAGAAAACCGGTGAGACGGTCGCTCCAATCACCCCACAGCGCCTGGCTCAGGGATCTCAcgccgc                                                                                                                               >  1:224352/3‑143 (MQ=255)
  gcgTCGTTAGTATGCTGACGACCAGACACATCCAGCACCTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCgccgcc                                                                                                                              <  1:245503/143‑1 (MQ=255)
       gtGAGCATGCTCACGACCAGATACATCCGTCATCTCCTGGTCATTCACCAGTTCGCCTAAAGTAATGTTGGTGAGAACCCCGGTGAGACGGTCCCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGc                                                                                                                         <  2:95127/143‑1 (MQ=255)
             aTGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAgccgc                                                                                                                   <  1:255503/143‑1 (MQ=255)
                    cgacCAGCCACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                                            <  2:105/143‑1 (MQ=255)
                           acacATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCCGACaa                                                                                                     >  2:7889/1‑143 (MQ=255)
                                 ccAGCACTTCCTGGCTAGCAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGTCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTg                                                                                               <  1:191875/143‑1 (MQ=255)
                                    gCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCa                                                                                            >  1:100543/1‑143 (MQ=255)
                                     cACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGACTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGCGGCAt                                                                                           >  2:283460/1‑143 (MQ=255)
                                                 attaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGt                                                                               >  1:225600/1‑143 (MQ=255)
                                                     aCCAGTTCGCCTCAACTAACTTTGCTGAGAAAACCGGTGAGCCGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCCGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAAc                                                                           <  1:237238/143‑1 (MQ=255)
                                                       cAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCTGCTTCTGCCCTGACAACGGGTGGCATCTACAGATTCGTCAACgg                                                                         >  1:26580/1‑143 (MQ=255)
                                                                       aTGTTGTTGAGAAAACCGGTGAGGCGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGcc                                                         <  1:181431/143‑1 (MQ=255)
                                                                         gttgttGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCaa                                                       <  2:172789/143‑1 (MQ=255)
                                                                          ttgttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAGTCGTCAACGGCGCTAATTACCTCGCCCAc                                                      >  2:19153/1‑143 (MQ=255)
                                                                                               cGGTCGTTCCAAGCACGCCACAGCGCGTGGTTCAGTGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt                                 <  2:266664/143‑1 (MQ=255)
                                                                                                         aaTCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCATGGCAGCCGCCATTCCCCTGACAACGGGTGGCATCGAGAGATTCGTAACCGGAGATAAGTGCCTCGCCATCGGAGACGCGACTGGCACCTTTTCCTACCAgc                       >  2:259007/1‑142 (MQ=255)
                                                                                                             aCGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAac                   >  2:40490/1‑143 (MQ=255)
                                                                                                               gCCCCAGCGCGTGGGTCAGGGATTTAGCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAacca                 <  2:73819/143‑1 (MQ=255)
                                                                                                                              tCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACACATAACACCCGCCTGGTCCACGt  >  1:321292/1‑143 (MQ=255)
                                                                                                                                   |                                                                                                                                         
GCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGT  >  NC_000913/2661690‑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: