Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,821 C→G 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,8210CG94.4% 93.5 / NA 37C104S (TGC→TCC) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (16/18);  minor base T (1/1);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCAC  >  NC_000913/2661686‑2661956
                                                                                                                                       |                                                                                                                                       
cGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATTCAGCTCTTTCTGGTTATTAATCAGTTCGTCTAAATTTTTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATc                                                                                                                                  <  1:445358/143‑1 (MQ=255)
               gTATGCTGACGGCCAGACACACCCAGCACTTGCTGTTTATTAACCAGTTCGCCTAAATTAATGTTGTTGAGAAACCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAgcc                                                                                                                   <  2:783708/143‑1 (MQ=255)
                 aTGCTGACGTCTAGACTCATCTAGAAGTTCTTTGTTATTAACGTGTTCGCCTTATGTTATGTAGTTAAGAAAGTCGTTGAGACGGTCGCTCAAATCACGCCACAGCGCTTTTGTCAGTTATTTATCGCCGCCCTGGCAgccgc                                                                                                                 <  2:689633/143‑1 (MQ=255)
                 aTGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAgccgc                                                                                                                 <  2:338843/143‑1 (MQ=255)
                    cTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTCAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCAACAGCGAGTAGGACAGGAATTTATAACCGCCCTGGCAGCAGCCtt                                                                                                              >  1:97143/1‑143 (MQ=255)
                    cTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTACCGCCGCCCTGGCCGCCGCCtt                                                                                                              >  2:160211/1‑143 (MQ=255)
                         gacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTAAGGGATGTATCCCCGCCCTGGCAGCCGCCTCTAccc                                                                                                         >  2:568629/1‑143 (MQ=255)
                         gacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGAGGGTCAGGGATTTATCGACGCCCTGGCAGCCGCCTTTAgcc                                                                                                         >  2:661847/1‑143 (MQ=255)
                          acCAGACACATCCAGCACTTCATGGTTATTAACCAGTTCGCCTAATTTTATGTTGTTGAGTAATCCGTTTTTACTGTTTCTTATTTTACGCCTCAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCt                                                                                                        <  1:113135/143‑1 (MQ=255)
                            tggACACATCCGGGCCTGTCTGGTGATTAACTGGTTCGGCTAGAGTAATGTTGTTGAGAAAACGGGTGAGACGGTCGCTCACATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                                                                                                      <  2:345913/141‑1 (MQ=255)
                                    tCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGt                                                                                              >  2:698415/1‑143 (MQ=255)
                                       aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGc                                                                                           <  2:395895/143‑1 (MQ=255)
                                          aCTTCCTGGTTATTAACCAGTTCGCCTCAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTGAAATCACGTCACAGCGCGTGGGGCAGGGTTTTATAGCGGCCCTGGCAGCGGCCTTTAGCCTGACAACGGGTCGCATc                                                                                        >  2:543725/1‑143 (MQ=255)
                                                   ttattaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACAGGTGAGACGGTCGCTCAAATCACGACACAGCGCGTGGGTCAGGGATTTATCGCCGCCATAGCAGCCGACTTTACCATGACAAAGGGTTGCATCTACAGATTc                                                                               >  2:140492/1‑143 (MQ=255)
                                                         aCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTAGCATCTACAGATTCATCAAc                                                                         >  1:536080/1‑143 (MQ=255)
                                                                  ccTAAAGTAAGGTTGTTTAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCGGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAt                                                                <  2:307015/143‑1 (MQ=255)
                                                                   cTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGACCTGGCAGCCGCCTTTACCCTGGCCACCGGTTGCATCTCCAGATTCGGCACCGGCGCTACtt                                                               >  2:97986/1‑143 (MQ=255)
                                                                    tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACAGCGCAAATTa                                                              >  1:549143/1‑143 (MQ=255)
                                                                             gttgttGAGAAAACCGGTGAGAAGGTCGCTCAAATCACAACACAGCGCGTGGGTAAGGGATGTATCGCCGCCCTGGCAACAGCCTTTAACCTGAGAACGGGAGGCATCTACAGAATCGTCAACAGCAGTAAGTACTTGTCCaa                                                     >  2:432350/1‑143 (MQ=255)
                                                                               tgttgAGAAAAACGGAGAGACGGTAGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACg                                                   <  1:327436/143‑1 (MQ=255)
                                                                                 ttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCGACGGCGCTAATTACTTAGCCaaacgc                                                 >  1:725693/1‑139 (MQ=255)
                                                                                  ttataacaCCTGTGAGACGGTCGCGCAAGACACGACACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCg                                                <  1:645882/136‑1 (MQ=255)
                                                                                     gaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGTCCTAATTACTCCGCCACCGGCGAtg                                             >  1:724676/1‑143 (MQ=255)
                                                                                      aaaaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAtgc                                            >  1:600561/1‑143 (MQ=255)
                                                                                           cGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGctgctg                                       <  1:698415/143‑1 (MQ=255)
                                                                                              tGAGACGGTCGCTCAAATCACGCCACATCGCGTCGTTGAGGTATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGCCGATGCTGCTGGcc                                    >  1:855010/1‑142 (MQ=255)
                                                                                                   cGGTCGCGCAAATCACGCCACGGCGCGTGGGTCAGGGATGTAGCGCCGCCCTGGTTGTCGCCGTCACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt                               <  1:191114/143‑1 (MQ=255)
                                                                                                     gTCGCTCAAATCGCGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                             <  2:844060/143‑1 (MQ=255)
                                                                                                     gTCGCGCAAATCACGCCACAGCGCGGGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                             <  1:144884/143‑1 (MQ=255)
                                                                                                     agcGCACAAATCACGCCACAACACAACCATCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc                             <  1:344960/141‑1 (MQ=255)
                                                                                                          ttttaTCACGCCACAGGGCGTGGGTCAGGGATTTAGCGCCGCCGTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc                        <  1:170046/139‑1 (MQ=255)
                                                                                                          tCAAATCACGCTACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCAGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc                        <  1:177394/143‑1 (MQ=255)
                                                                                                            gagtCACGCCCCAGCGCGTGGGTCCGGGATCTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTTGCATCTACGGCTTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAg                      <  2:770527/140‑1 (MQ=255)
                                                                                                              aTCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTGTGCCTACCAGac                    >  2:160846/1‑142 (MQ=255)
                                                                                                                        aGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAGCCTGCTGGCCTCTTTGCCTCACAGATAACCGCGGcc          >  2:486014/1‑143 (MQ=255)
                                                                                                                             gTGGGTCAGGGATTTATCGCCGCCCGGGCAGCCGCCTTTACCCGAACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAATACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTc     <  2:657715/143‑1 (MQ=255)
                                                                                                                                ggTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATAACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAAACACCGCACGGTcccc  >  2:322623/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCAC  >  NC_000913/2661686‑2661956

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: