Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,029 G→C 35.8% Q35E (CAG→GAG)  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,662,0290GC35.8% 37.7 / 33.5 44Q35E (CAG→GAG) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base G (10/18);  new base C (7/9);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTT  >  NC_000913/2661887‑2662170
                                                                                                                                              |                                                                                                                                             
gcgcTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGACAAGAAAGGGCAATTCCCTg                                                                                                                                               >  1:109400/1‑143 (MQ=255)
gcgcTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTCAACGCAGGGGGGAAAGCACCATTCCTAGATAAGAAAGGGAAATTCCCTg                                                                                                                                               >  1:513045/1‑143 (MQ=255)
  gcTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTCGTTACAGATAAGAAAGGGAAATTCCCTGAc                                                                                                                                             >  2:561994/1‑143 (MQ=255)
    tAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATACCCGCACGt                                                                                                                                           >  1:758688/1‑143 (MQ=255)
        tACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg                                                                                                                                       >  1:90267/1‑143 (MQ=255)
                cAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCCTTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGCAAGGGAAATTCCCTGACGTTCGGACATCTCa                                                                                                                               >  2:415423/1‑143 (MQ=255)
                cAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTCTACGCAGACGGGAAAACAGTTGTTCCAGATAACAAAGGGAAATTCCATGACGATCGGAAATATCa                                                                                                                               >  2:78283/1‑143 (MQ=255)
                cAACGGCGATGCTGCCGGCATCTTTGCCTAACAGATAACCCCCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa                                                                                                                               <  1:665456/143‑1 (MQ=255)
                agacgcggATGCTGCTGGCAGCGTGGCCTAACAGATGACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTGTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa                                                                                                                               <  2:791949/136‑1 (MQ=255)
                      cGATGCTGCTGGCATCTTTGCCTAACAAATAACCGCCGCCTGGCCCACGCACACTGGAAACCAGACCATCCCTACGCAGACGGGAACACAGATGACCCACATACGAAAGGGCAACCCCCACAACAACGGAAATATCAGCCacc                                                                                                                         >  2:100830/1‑141 (MQ=255)
                      agATGCAGCCGCCATCCTAGCCTAACAGATAACCCCCGCCTGGCACACGTACGCTGCAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAAc                                                                                                                         <  1:68141/142‑1 (MQ=255)
                                 gCATCTCTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTCGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTCTTCCACCTAAGAAAGGGAAATTCCCCGACGTTCGCAAATATCAGCCAACGGTACCGGGcc                                                                                                              <  2:97630/143‑1 (MQ=255)
                                 gCATCTCTGCCCAACAGATAACCACCGCCCGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGcc                                                                                                              <  1:740/143‑1 (MQ=255)
                                      tttGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGCCCGAACATATCAGCCAACGGTCCCGCGCCCGcct                                                                                                         >  2:715868/1‑141 (MQ=255)
                                       ttGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTc                                                                                                        <  2:349109/143‑1 (MQ=255)
                                          ccTAACAGATAAGCACCGCCGGGGCCACGGACGCTGGAAACCAGACCATTTTGACGCAGACGGGAAAACAGTTGTTCCAGATAAGGAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCaga                                                                                                     <  2:138254/143‑1 (MQ=255)
                                               cAGATGACCACCGCCTGGTCCACGTACGCGGGAAACCAGACCAGTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGTCAACGGTACCGGGCCCGCTTCAGAGTTGa                                                                                                <  2:428126/143‑1 (MQ=255)
                                                 gATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGCCAACAGTGGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAgc                                                                                              >  2:137054/1‑143 (MQ=255)
                                                            ccTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCagg                                                                                   >  1:609710/1‑141 (MQ=255)
                                                                gTCCACGTACGCTGGAAACCAGACCATTTTTTTTTTTTTTTGAAAACAGTTGCTCCAGATAAGAAACGGATTTTCCCTGACGTTCGGGAATATCAGCCATCGGTACCGGGCCCGCTTCAGAGTTGACCGCAACGTCAAGCAtt                                                                               >  2:124427/1‑143 (MQ=255)
                                                                 tCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTg                                                                              >  1:357553/1‑143 (MQ=255)
                                                                  ccACGTACTCTGTAAACCAGCCCATTTTTCCGCAGTCGTGAAAACAGTTGTTTCAGATAAGAAAGGGAAATTCCTTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGTGCGTAACGTGAAGCATTGc                                                                             <  2:179353/143‑1 (MQ=255)
                                                                   cacgcacgCTCGAACCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCg                                                                            <  1:825864/143‑1 (MQ=255)
                                                                         cGCTGGAAACCAGACCATTTTCACGCAGACGGGAAAACAGTTGTACAAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACg                                                                      <  1:688365/143‑1 (MQ=255)
                                                                             ggAAACCAGACCGTTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCAt                                                                  <  1:93992/143‑1 (MQ=255)
                                                                             ggAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCAt                                                                  >  2:76007/1‑143 (MQ=255)
                                                                                   cAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGcc                                                            >  1:502108/1‑143 (MQ=255)
                                                                                    aGACCTGTTTTGGGCAGACGCGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGccc                                                           <  2:752620/143‑1 (MQ=255)
                                                                                           tgtttcGCTGACGGGAAGACAGTTGTTCCAGATAAGAAGTGGAAATTCCCTCACGTTCGGTAGTATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGAt                                                    <  2:622667/138‑1 (MQ=255)
                                                                                               aCGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCa                                                <  1:487180/143‑1 (MQ=255)
                                                                                                          aaaaCAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACTGCATAGCGCCCTTTAGATATCCGTCTCACGTcc                                     >  2:145860/1‑142 (MQ=255)
                                                                                                            aaCAGTTGCTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTa                                   <  1:291350/143‑1 (MQ=255)
                                                                                                                    ttCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCCTAGCGCCCGTTAGATGACAGTCTCGCGTCTTACTTccacc                           >  2:241763/1‑139 (MQ=255)
                                                                                                                    ttCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGGTGGGCGCCCCGTGAAGCATTGCGGGCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCt                           <  1:561994/143‑1 (MQ=255)
                                                                                                                    ttCCAGATAAGAAAGGGAAATGCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCt                           <  2:109400/143‑1 (MQ=255)
                                                                                                                    gtCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCt                           <  2:513045/142‑1 (MQ=255)
                                                                                                                          atcaGAAAGGGACTTGCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCCAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACt                     <  2:457561/140‑1 (MQ=255)
                                                                                                                               aaggttaaaTGCCCTGTTGTTCGAAAATATCAGCCAACGGTACCGTGCCCGCTTCAGAGTTGAGCGCAACGTCAAGTATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGccc                <  2:505101/137‑1 (MQ=255)
                                                                                                                               aaaGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGccc                <  1:292419/143‑1 (MQ=255)
                                                                                                                               aaaGGGAAATTCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGccc                <  1:463675/143‑1 (MQ=255)
                                                                                                                                 aGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTATCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTGACCTCAAACTCGCCGCt              >  1:726679/1‑143 (MQ=255)
                                                                                                                                  gggAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGGTCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTCCTTCACCTCACCCTCGCCCCTg             >  2:91962/1‑143 (MQ=255)
                                                                                                                                        tgcccTCACGTGCGGAACTATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCggg       <  2:650159/141‑1 (MQ=255)
                                                                                                                                         tCCCTCACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGCGTTGATGGCAACTTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCgggg      <  1:818488/143‑1 (MQ=255)
                                                                                                                                            cTGACGTCCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGATGCGCTCAACGTCAAGCATTGCGTTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGttt   <  1:780793/143‑1 (MQ=255)
                                                                                                                                             tGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGtttt  <  1:410029/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTT  >  NC_000913/2661887‑2662170

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