Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,816 T→G 64.1% T106P (ACC→CCC)  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,8160TG64.1% 33.0 / 27.9 42T106P (ACC→CCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (11/4);  new base G (15/12);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCC  >  NC_000913/2661676‑2661948
                                                                                                                                            |                                                                                                                                    
tgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCGCTTCCTGGTTATTGACCAGTTCGCCTAAAGTAATGTTGTTGAGAACACCGGTGAGACGGTCGCTCAAATCCCGCTACAGCGCGCGGGGCa                                                                                                                                    >  1:428126/1‑143 (MQ=255)
tgtgCGGGTGCGTGGCGCGTCGTGAGTATGCGGACGACCAGACACATCCAGCACTTCCTGGTTATTAGCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCa                                                                                                                                    <  2:549971/143‑1 (MQ=255)
          cGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATc                                                                                                                          >  1:97630/1‑143 (MQ=255)
                gcgTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACCGCGCGTGGGTCAGGCATTTATCgccgcc                                                                                                                    >  2:567805/1‑143 (MQ=255)
                   tcgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTg                                                                                                                 >  1:34028/1‑143 (MQ=255)
                   tcgtGAGTATGCCGACGACCAGACACATCGAGCACTTCCTGGTTAGTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTg                                                                                                                 <  2:370671/143‑1 (MQ=255)
                    cgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTgg                                                                                                                <  1:302566/143‑1 (MQ=255)
                           ctGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAgccgc                                                                                                         <  2:228938/142‑1 (MQ=255)
                              cTGACGACCAGACACATCCAGCACTTCCTGGTTATTCCCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCtt                                                                                                      <  2:571549/143‑1 (MQ=255)
                                      cAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                                                                                              <  1:286322/143‑1 (MQ=255)
                                       aGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGAc                                                                                             >  1:622667/1‑143 (MQ=255)
                                         acacATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCCCAAATCACGCCACAGCGCGTGGGGCAGGCATCTATCGCCGCCCGGGCCGCCGCCTTTACCCTGACag                                                                                           >  2:463675/1‑142 (MQ=255)
                                            caTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACgg                                                                                        >  1:457561/1‑143 (MQ=255)
                                             aTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGCTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGCGCAGGCATTTATCGCCGCCCTGGCACCCGCCTTTACCCTGACAGCggg                                                                                       >  2:291350/1‑143 (MQ=255)
                                              tCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                      <  1:477886/143‑1 (MQ=255)
                                              tCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                      <  1:212852/143‑1 (MQ=255)
                                                    aCTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCATCCGCCTTTTCCCTTACTCCGGGTAGCAta                                                                                >  1:12149/1‑142 (MQ=255)
                                                       gccTGGTTAGTACCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTAc                                                                             <  1:492756/142‑1 (MQ=255)
                                                         cTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCTGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCCTCTACAg                                                                           >  2:426530/1‑143 (MQ=255)
                                                                ttaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGGCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGAGAACGGGTGGCATCTACAGATGCGtg                                                                    >  2:478102/1‑142 (MQ=255)
                                                                  aaCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCaa                                                                  >  2:730956/1‑143 (MQ=255)
                                                                    ccAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACCACGGGCGGCATCTCCAGATTCGTCAACg                                                                >  2:818488/1‑143 (MQ=255)
                                                                             cTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGCCCTCTACAGATGAGTCAACGgccataaat                                                       >  2:148817/1‑136 (MQ=255)
                                                                             cTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTACtt                                                       >  1:349109/1‑143 (MQ=255)
                                                                              tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGACGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACGGATTCGTAAACGGCGCTAAtga                                                      >  1:343627/1‑141 (MQ=255)
                                                                                aaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATGTATCGGCGCCCTGGCAGCCGCAGTGACCCTGAGAACGGGTGGGGTCTGCAGATGCGGCGACGGCGGTAATGAcg                                                    >  2:825864/1‑142 (MQ=255)
                                                                                 aGTAATGTTGTTGAGAAACCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACtt                                                   <  1:495924/143‑1 (MQ=255)
                                                                                 aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACCGGTTCGCCAACGGCGCTAATTACtt                                                   >  2:487180/1‑143 (MQ=255)
                                                                                           ttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAACTACTTCGCCAACggg                                         >  2:740/1‑142 (MQ=255)
                                                                                                  aaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAtgctg                                  <  1:51721/143‑1 (MQ=255)
                                                                                                  aaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTCCTTCGCCAAGGACgctgcta                                  >  1:752620/1‑142 (MQ=255)
                                                                                                    ccGGTGGGACGGTCGCTCCAATCCCGCCCCAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATgctgct                                <  2:694507/143‑1 (MQ=255)
                                                                                                            aCGGTCGCTCAAATCACGCCTCAGCTCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGCCAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCtt                        <  2:77035/143‑1 (MQ=255)
                                                                                                            aCGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCtt                        >  1:650159/1‑143 (MQ=255)
                                                                                                               gTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACCACGGGTGGCATCTACAGCTACGTCAACGGCGCTAATTACTTCGCCAACGGCGGTGCTGCCGGCACCTTTGc                     >  2:172174/1‑143 (MQ=255)
                                                                                                                tCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGcc                    <  1:158975/143‑1 (MQ=255)
                                                                                                                  gCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTa                  >  1:155656/1‑143 (MQ=255)
                                                                                                                   cTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTaa                 <  1:730956/143‑1 (MQ=255)
                                                                                                                      aaaTCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACACTGACAACGCGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAg              >  2:410029/1‑143 (MQ=255)
                                                                                                                               cacaGCGCGTGGGTCAGGCTTTTATCGCCGCCCTGGCAGCCGCCTTTACACTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCAGGAATCTTTGCCTAACAGATAaccacc     >  1:201339/1‑143 (MQ=255)
                                                                                                                                 caGCGCGTGGGGCAGGCGTTTATCGCCGCCCCACCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGc   <  1:777087/143‑1 (MQ=255)
                                                                                                                                  aGCGCGTGGGTCAGGCATTTATCGCCGCCCTGCCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGCTTCGTCAACGGCGCTCCTTACTTCGCCCACGGCGATGCTGCTGGCCTCTTTGCCTAACAGCTAACCCCCGcc  >  1:138254/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                    
TGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCC  >  NC_000913/2661676‑2661948

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: ‑