Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,662,004 T→C E43G (GAA→GGA)  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,0040TC96.5% 539.1 / 6.5 144E43G (GAA→GGA) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (73/65);  total (77/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 159 total aligned reads displayed.

GTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAA  >  NC_000913/2661957‑2662050
                                               |                                              
gTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTccc                                              <  1:4351000/50‑1 (MQ=255)
 tACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCa                                             <  1:1192626/50‑1 (MQ=255)
 tACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCa                                             <  1:1780199/50‑1 (MQ=255)
   cGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGa                                           <  1:374324/50‑1 (MQ=255)
   cGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAg                                            <  2:1748527/49‑1 (MQ=255)
   cGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAg                                            <  2:2352720/49‑1 (MQ=255)
   cGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAg                                            <  2:4578668/49‑1 (MQ=255)
    gCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          <  1:574618/50‑1 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGAt                                          >  2:2541622/1‑49 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  1:1237630/1‑49 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  1:4368995/1‑49 (MQ=255)
      tGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  2:3816580/1‑48 (MQ=255)
      tGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  2:1953350/1‑48 (MQ=255)
      tGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  2:2785223/1‑48 (MQ=255)
       ggAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  1:2020588/1‑47 (MQ=255)
       ggAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                          >  2:1239838/1‑47 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATaa                                        <  2:742144/48‑1 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATaa                                        <  1:3431290/48‑1 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAg                                       >  1:1362634/1‑49 (MQ=255)
         aaaCCAGACCATTTTTACGCAGCCGGGACAACAGTTGTACCAGATAAGa                                      >  2:1822067/1‑49 (MQ=38)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGaa                                     <  1:5488754/50‑1 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGaa                                     <  2:3228191/50‑1 (MQ=255)
           aCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAg                                   <  1:3202644/50‑1 (MQ=255)
           aCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAg                                   <  2:4034893/50‑1 (MQ=255)
           aCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAg                                   >  1:1680365/1‑50 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAgg                                  >  1:1071828/1‑50 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                  >  1:1033414/1‑50 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                  <  1:1601248/50‑1 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                  >  2:2709035/1‑50 (MQ=255)
             cAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAggg                                 <  2:823725/50‑1 (MQ=255)
             cAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                  >  1:4164013/1‑49 (MQ=255)
             cAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                  >  2:2650693/1‑49 (MQ=255)
              aGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAggg                                 >  2:1201971/1‑49 (MQ=255)
               gACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGaa                               <  2:5532062/50‑1 (MQ=255)
               gACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGa                                >  2:2872401/1‑49 (MQ=255)
                aCCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGaaa                              <  1:4659069/50‑1 (MQ=255)
                  cATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAtt                            <  1:339324/50‑1 (MQ=255)
                  cATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAtt                            >  2:4373841/1‑50 (MQ=255)
                    tttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTcc                          >  2:5040308/1‑50 (MQ=255)
                     ttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                         >  2:3157703/1‑50 (MQ=255)
                     ttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                         >  1:2967045/1‑50 (MQ=255)
                      tttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                         <  2:2058358/49‑1 (MQ=255)
                      tttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                         <  1:2503537/49‑1 (MQ=255)
                      tttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                         <  2:5501806/49‑1 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCt                        <  1:4965083/49‑1 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCt                        <  1:4829211/49‑1 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCt                        <  1:2498883/49‑1 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGGCCCAGATAAGAAAGGGAAATTCCCt                        <  1:5108983/49‑1 (MQ=39)
                        tACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGa                      >  2:2953163/1‑50 (MQ=255)
                        tACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGa                      >  2:425521/1‑50 (MQ=255)
                         aCGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGa                      <  1:3196545/49‑1 (MQ=255)
                         aCGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGAc                     >  2:1344533/1‑50 (MQ=255)
                          cGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGAc                     <  2:3596706/49‑1 (MQ=255)
                           ggaGACGGGCAACCCGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                   >  2:5008122/3‑50 (MQ=255)
                           gCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                   <  1:5392700/50‑1 (MQ=255)
                           gCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                   >  2:1850832/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGCTAAGAAAGGGAAATTCCCTGACGtt                  >  2:879793/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                  <  1:4159673/50‑1 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                  >  1:220291/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                  <  1:2318626/50‑1 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                  <  1:5017309/50‑1 (MQ=255)
                             aGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTc                 <  1:796035/50‑1 (MQ=255)
                             aGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTc                 <  2:71747/50‑1 (MQ=255)
                             aGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTc                 <  2:3134763/50‑1 (MQ=255)
                             aGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTc                 >  2:2553716/1‑50 (MQ=255)
                             aGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTc                 <  2:1110911/50‑1 (MQ=255)
                              gACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                <  1:2839864/50‑1 (MQ=255)
                              gACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                >  1:3727003/1‑50 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg               >  2:4091669/1‑49 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg               >  1:1926464/1‑49 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa              <  1:477747/50‑1 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa              <  2:1408620/50‑1 (MQ=255)
                                 gggAAAGCAGATGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa              >  2:3794261/1‑49 (MQ=255)
                                 gggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa              >  1:4524692/1‑49 (MQ=255)
                                  ggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa            >  2:1560720/1‑50 (MQ=255)
                                  ggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa            <  2:1040528/50‑1 (MQ=255)
                                  ggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa            >  1:4815068/1‑50 (MQ=255)
                                   gAAAACAGTTGTTCCAGATAAGAAAGGGAAATTACCTGACGTTCGGaaa            >  2:975280/1‑49 (MQ=255)
                                   gAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa            >  2:2408036/1‑49 (MQ=255)
                                   gAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa            >  1:5166601/1‑49 (MQ=255)
                                    aaaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata          <  2:4300956/50‑1 (MQ=255)
                                     aaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat         >  2:4445131/1‑50 (MQ=255)
                                     aaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat         >  2:673543/1‑50 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat         >  2:2276179/1‑49 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat         >  2:3863029/1‑49 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATc        <  2:1595183/50‑1 (MQ=255)
                                       aCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat         >  2:1009145/1‑48 (MQ=255)
                                       aCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa       <  1:4389071/50‑1 (MQ=255)
                                       aCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa       <  1:1314586/50‑1 (MQ=255)
                                        cAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAg      <  1:1167443/50‑1 (MQ=255)
                                        cAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAg      >  2:5327264/1‑50 (MQ=255)
                                         aGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGc     >  1:4956457/1‑50 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc    <  2:3559780/50‑1 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGc     >  2:3965992/1‑49 (MQ=255)
                                          ggtgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGc     >  2:3995214/3‑49 (MQ=255)
                                           ttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc    >  2:5492055/1‑49 (MQ=255)
                                           ttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc    >  2:5552794/1‑49 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa  >  1:3061683/1‑50 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa  >  1:32250/1‑50 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa  >  2:1444093/1‑50 (MQ=255)
                                               |                                              
GTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAA  >  NC_000913/2661957‑2662050

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: