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,8210CG99.5% 67.9 / ‑7.1 26C104S (TGC→TCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (12/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

GTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCA  >  NC_000913/2661694‑2661955
                                                                                                                               |                                                                                                                                      
gtcgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCt                                                                                                                         >  1:307258/1‑143 (MQ=255)
           gCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTCCGCCTAAAGTAACGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGCCAGGGCTCTATCGCCGCCCCGgccgccgccc                                                                                                              >  1:226789/1‑142 (MQ=255)
                            tCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                             <  1:325995/143‑1 (MQ=255)
                             ccAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTg                                                                                            >  1:151436/1‑143 (MQ=255)
                             ccAGCACTACCTGGTTATTATCCAGTTCGCCTAAAGGAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTg                                                                                            <  1:660/143‑1 (MQ=255)
                              cAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTCCCCTGACAACGGGTgc                                                                                           >  2:20684/1‑142 (MQ=255)
                                    ttCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTa                                                                                     >  1:61558/1‑143 (MQ=255)
                                       cTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCCCCCTGGCCGCCGCCTTTGCCCTGCCAACGTGAGGCACCTCCAg                                                                                  >  2:393220/1‑143 (MQ=255)
                                               taACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCCGCGCTTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGACATCTACAGATTCGTCa                                                                          >  1:42025/1‑143 (MQ=255)
                                                aaCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTCGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATCCGTCaa                                                                         >  2:367424/1‑143 (MQ=255)
                                                   cggTTCTCCGAAAGTAATGTTGTTGAGAGAGCCGGTGAGACGGTCGTTCAAATTACGCCCCAGCGCGTGGGCCAGGGATTTATCGCCGCCCCGGCAGCAGCTTTTGCCCTGACAACGGGTGGCATCTACAGATTCGTCAACgg                                                                      <  2:164478/141‑1 (MQ=255)
                                                      ttCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGgcgc                                                                   >  2:422179/1‑143 (MQ=255)
                                                                           gagaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTACTTACTTCGCCAACGGCGa                                              >  1:146497/1‑143 (MQ=255)
                                                                               aaaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATAACCGCCCTGGCAGCCGCCTTTACCCTGACCACGAGCGGCATCTACAGATTCGTCAACGGCGCGAATTGCTTCCCCAAAAGCGAtagg                                          >  2:176860/1‑140 (MQ=255)
                                                                                           cGGTCGCGCACATCACGCCGCGGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt                              <  2:265974/143‑1 (MQ=255)
                                                                                             ggcgggcAAATGACGCCAGAGCGCGTGGGGCAGGGATGTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                            <  1:180563/137‑1 (MQ=255)
                                                                                             ctctcgtaaaTCACGCCGCGGCGCTGGGGTCAGGGTTTTATCGCCGCCCTGGTAGCCGCCGTTACCGTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                            <  2:54861/136‑1 (MQ=255)
                                                                                                 cTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTaa                        <  1:510362/143‑1 (MQ=255)
                                                                                                     aaTCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGa                    <  1:339226/143‑1 (MQ=255)
                                                                                                     aaTCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGa                    <  1:74938/143‑1 (MQ=255)
                                                                                                       tCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATa                  <  2:151436/143‑1 (MQ=255)
                                                                                                        cgcGCCCCAGCGCGTGGGTCAGGCATTTCTCGCCGCCCTGGCAGCCGCCTGTACCCTGACAACGGGTGGAATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATaa                 <  2:425989/141‑1 (MQ=255)
                                                                                                        cACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATaa                 <  1:354651/143‑1 (MQ=255)
                                                                                                                     gTGGGTCGGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAGCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAGGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTc    <  2:39269/143‑1 (MQ=255)
                                                                                                                     gTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTc    <  1:168565/143‑1 (MQ=255)
                                                                                                                       gggTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCa  >  1:400767/1‑143 (MQ=255)
                                                                                                                               |                                                                                                                                      
GTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCA  >  NC_000913/2661694‑2661955

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