Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,816 T→G 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,8160TG97.8% 121.2 / NA 46T106P (ACC→CCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (25/20);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTC  >  NC_000913/2661685‑2661953
                                                                                                                                   |                                                                                                                                         
aagTGGAGCGCCGGGAGGACGCAGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCGAAAGTAAGGTTGTGGAGAAAACCGGTGAGACGGGCGCACAAATCACACCACAGCGCGTGGGGCAGGCATTTAt                                                                                                                                <  1:246857/141‑1 (MQ=255)
 ggTGGCGAGTCGTGAGTATGCTCACGACCAGACACATCGTGCACTTCCTGGTTATTAACCAGTTAGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATc                                                                                                                               <  1:992677/142‑1 (MQ=255)
     gcgcGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCCCAAATCACGCCACAGCGCGCGCGGCAGGCATTTATcgccg                                                                                                                           >  2:1262226/1‑143 (MQ=255)
       gcgTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGACTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCgccgcc                                                                                                                         >  1:659330/1‑143 (MQ=255)
                    gCTGACGACTATACACATCCAGTACTTCCTTGTTATTACCCAGTTCGCCTAAAGGAAGGTCGTGGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCt                                                                                                            <  1:1157602/143‑1 (MQ=255)
                        acgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGATAGCTCAAATAACGCCACAGCGCGTGGGGCAGGCATTGATCGCAGCCCTGGCAGCCGCCTTTAc                                                                                                        >  1:37635/1‑143 (MQ=255)
                         cgcgcAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                                       <  2:683808/139‑1 (MQ=255)
                         ccacCAGACACATCCATCTCTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                                       <  2:461570/141‑1 (MQ=255)
                               gACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCAAAGAGCGTGGAGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACa                                                                                                 >  1:997246/1‑143 (MQ=255)
                                   caTCCAGCACTTCCTGGTTATTACCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACgg                                                                                             <  2:225507/143‑1 (MQ=255)
                                    aTCCAGCACTTTCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGTGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGGGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACggg                                                                                            <  1:716557/143‑1 (MQ=255)
                                    aTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATTTCGGAACAACGCGTGGGGCAGGCATACATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACggg                                                                                            >  1:446479/1‑143 (MQ=255)
                                     tCCAGCTCTGTCGGGTTCTTAGCTAGTTCGGCTCAAGTAATGTTGTTGAGTATACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                           <  2:959295/143‑1 (MQ=255)
                                      ccAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTg                                                                                          >  2:1104734/1‑143 (MQ=255)
                                       taGCACTTTCTGGTTATTAACCTTTTCGCCTAAAGTAATGTTGTTTTGAAATCCGGTTAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGAATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTgg                                                                                         <  1:631615/142‑1 (MQ=255)
                                        aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAATCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCACTGGAAGACGCCTTTACCCTGACAACGGGTGGc                                                                                        >  1:686584/1‑143 (MQ=255)
                                                 tGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGa                                                                               >  1:754442/1‑143 (MQ=255)
                                                    ttattaTCCAGTTCTCTTAAGGTAATGTTGTTGTGAAATCCGGTGAGACGGTCGCTCAAATTACGCCACAGCGCGTGGGGCAGGCATGTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACATCGGGTGGCATCTACAGATTc                                                                            <  2:584648/143‑1 (MQ=255)
                                                      attaattaGTTCGCTTAAATTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCGAGAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGt                                                                          <  1:1202916/143‑1 (MQ=255)
                                                       ttaACCATTTCGCCTTATGTATTTTTTTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                                         <  1:1258999/143‑1 (MQ=255)
                                                       ttaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                                         <  2:754442/143‑1 (MQ=255)
                                                         aaCCAGTTCGCCTAAATTAATGTTGTTGAGCAATTCGTTTAGATGGTTGCTCAAATCACGACACAGCGCGTGGGGCAGGCATTGATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCaa                                                                       <  1:52740/143‑1 (MQ=255)
                                                         aaCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTAAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTCGCAGCCGCCGTTACCCTGACAACAGGTGGCACCTAAAGATTCGTCaa                                                                       >  2:479754/1‑143 (MQ=255)
                                                           ccAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCACACTGGCAGCCGCCTTTACACTGACAACGGGAGGCATATAAAGATTCGTCAACg                                                                     >  1:1098675/1‑143 (MQ=255)
                                                              gTTCGCCTAAAGTACTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCCTTTATCGCCGCCCAGGCAGCCGACTTTACCCTGACAACGGGTGGCATCCACAAATTCGTCAACGgcg                                                                  >  2:146748/1‑143 (MQ=255)
                                                                    atAAACTAATGTTGTTGAGAAAACCGGTTATTCGGTCGCTCTAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAtt                                                            <  1:721766/142‑1 (MQ=255)
                                                                     tAAAGTAATCATGTTGAGAAAACAAGAGAGACGGTCGCTCAAATCACGCCACAGCGCCTGAGGCAGGCATCTATCGCCGCCCTGGCAACCGCCTCTACCCTGACAACGCGTGGCATCTACAGCTTCGTCAACAACGCTaacga                                                           >  2:760722/1‑140 (MQ=255)
                                                                      aaaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAc                                                          >  2:775125/1‑143 (MQ=255)
                                                                        aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACGCTCACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACtt                                                        >  2:414437/1‑143 (MQ=255)
                                                                         gTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTAACACTAACAACGAGTGGCATCTAAAAATTCATCAACAGCGCTAATTACTTc                                                       >  1:539544/1‑143 (MQ=255)
                                                                                  ttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAATCGTCAACGGCGCTAATTACTTCGCCAACGGc                                              >  1:1223593/1‑143 (MQ=255)
                                                                                     agaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTGCGCCAACGGCGAt                                           >  2:1037216/1‑143 (MQ=255)
                                                                                     agaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATGACTCCGCCAAGGGCGAt                                           >  2:347317/1‑143 (MQ=255)
                                                                                         aaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCTTTCCCCCCGCCCCCCGCCCCCATCTACGGATCCGTCAACTGCACTAATTACTTCGAAAAAGGCGAtgata                                       >  1:908929/1‑140 (MQ=255)
                                                                                             ggTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTgg                                   >  1:812564/1‑143 (MQ=255)
                                                                                                  gaCGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAAAGAGAGGCAACAACAGATTCGTCAAACGCGCTAGGAACTTCGCCAACGGCGATGCTGCTGGCATCt                              >  1:970151/1‑143 (MQ=255)
                                                                                                   aCGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGTATTTATCGCCCTTCTATCATCCACCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCtt                             <  1:270079/143‑1 (MQ=255)
                                                                                                    gggTCGCTCAAATCACGCCACAGCGCGTGGGGCGGGAACATACCTCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt                            <  2:1223593/142‑1 (MQ=255)
                                                                                                    cGGTCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTTCCTTGATAACGTGTGGCATTTACAGATTCGTCAACGGCGCGAAGTACTGCGCCAACGGCGATGCTGCTGGCATCttt                            <  1:572202/143‑1 (MQ=255)
                                                                                                       tCGCTCAAATCACGCCACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGcc                         >  2:1236267/1‑143 (MQ=255)
                                                                                                         gCTCAAATCACGCGACAGCGCGAGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCACTAATTACTTCGCCAACGGCGAGGCAGCTGGCATCTTTGCCtc                       >  2:490095/1‑142 (MQ=255)
                                                                                                          cTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCATGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTaa                      <  1:881558/143‑1 (MQ=255)
                                                                                                                     ccACAGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTTCCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAGTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAaccac           <  1:935802/143‑1 (MQ=255)
                                                                                                                      cacaGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAGCAGATAaccacc          >  2:415751/1‑143 (MQ=255)
                                                                                                                        caGCGCGTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCTTGACAATGGGTGGTATTTAGAGATTCGTGAACGGCGCTAATTACGGGGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGc        <  1:1204849/143‑1 (MQ=255)
                                                                                                                              gTGGGGCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCTAAACAGATAACCACCGCCTGGTc  >  2:78641/1‑143 (MQ=255)
                                                                                                                                   |                                                                                                                                         
GCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTC  >  NC_000913/2661685‑2661953

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