Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,872 A→G 82.8% V87A (GTA→GCA)  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,8720AG82.8% 54.9 / 7.9 29V87A (GTA→GCA) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (11/13);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

TTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGAT  >  NC_000913/2661730‑2662010
                                                                                                                                              |                                                                                                                                          
ttCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCACCCGCCTGTACCCTGACAACGGGCGGCACCTa                                                                                                                                            >  1:323879/1‑143 (MQ=255)
            aaCCAGTTCGCCTAAATTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAACTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCaa                                                                                                                                <  1:283285/143‑1 (MQ=255)
               cAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTTAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCATGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACgg                                                                                                                             >  1:435932/1‑143 (MQ=255)
                     gCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGCCATCTGCAGATTCGTCAACGGCGCtcg                                                                                                                       >  1:297980/1‑141 (MQ=255)
                          aaGGAATGTTGTTGAGAAAACCGGTGCGAAGGCCACTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACt                                                                                                                  <  1:329484/143‑1 (MQ=255)
                           aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACtt                                                                                                                 >  1:22359/1‑143 (MQ=255)
                                tgttgtTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCa                                                                                                            >  1:225705/1‑143 (MQ=255)
                                         gaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCTGGCATTTATCGCCGCCCTGGCAGCAGCCTTTACCCTGAGAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTACTCACTTCGCGAAAGGCGAtg                                                                                                   >  1:121416/1‑143 (MQ=255)
                                              ccGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCTGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATgctgct                                                                                              <  2:392613/143‑1 (MQ=255)
                                               cGGTGAGACGGTCGCTCAAAGCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGctgctg                                                                                             <  2:207207/143‑1 (MQ=255)
                                                  tcagccGGTCGCTCGAACCACGCCACAGCCCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGCGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCa                                                                                          <  2:159251/138‑1 (MQ=255)
                                                        ggTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTg                                                                                    >  1:425119/1‑143 (MQ=255)
                                                         gTCGCTCAAATCACGCCGCAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                                                                   <  2:89536/143‑1 (MQ=255)
                                                         ctCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                                                                                   <  1:359640/142‑1 (MQ=255)
                                                               cAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACa                                                                             <  1:113484/143‑1 (MQ=255)
                                                                   tCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTGCGCCAACGGCGATGCTGCTGGCATCTTCGCCTAACacaga                                                                         >  2:175988/1‑139 (MQ=255)
                                                                           caGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGc                                                                 <  1:100310/143‑1 (MQ=255)
                                                                              cgcgTAGGTCTGGCAGTTATCGCCTCCCTGGCAGCCGCCGTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCGAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTg                                                              <  2:404129/143‑1 (MQ=255)
                                                                                         ggCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATCACCACCGCCTGGTCCACGTAcc                                                   >  2:203182/1‑142 (MQ=255)
                                                                                            aTTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATCACCACCGCCTGGTCCCCGTACGCTg                                                >  1:337839/1‑143 (MQ=255)
                                                                                                 tCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATGACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGCCCACGGCCGCTGGCAAc                                           >  1:120957/1‑143 (MQ=255)
                                                                                                                       tACCCTGGCAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACgg                     <  2:297980/143‑1 (MQ=255)
                                                                                                                          ccTGACAACGGGTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTGTTACGCAGACGGGaa                  >  1:8722/1‑143 (MQ=255)
                                                                                                                          ccTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGCTAACCACCGCCTGGTCCCCGTACGCTGGACACCCGACCCATTTTACGCAGACGGGag                  >  2:398108/1‑142 (MQ=255)
                                                                                                                          ccTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCCTTTTTACGCAGACGGGaa                  >  2:392560/1‑143 (MQ=255)
                                                                                                                                 aCGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgtt           <  1:132870/143‑1 (MQ=255)
                                                                                                                                    ggTGGCATCTGCAGATCCGTCAACTGCGCTGATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgttgtt        <  2:425119/143‑1 (MQ=255)
                                                                                                                                     gTGGCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTc       >  1:414516/1‑143 (MQ=255)
                                                                                                                                       ggCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCTACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCa     <  2:78862/143‑1 (MQ=255)
                                                                                                                                          aTCTGCAGATTCGGCAACGTCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGAt  <  2:130262/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGAT  >  NC_000913/2661730‑2662010

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: