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,0040TC94.4% 325.0 / 6.8 89E43G (GAA→GGA) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (40/44);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACG  >  NC_000913/2661957‑2662052
                                               |                                                
gTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTccc                                                <  1:2019571/50‑1 (MQ=255)
 tACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCa                                               <  1:2201468/50‑1 (MQ=255)
 tACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCa                                               <  1:1298202/50‑1 (MQ=255)
   cGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAg                                              <  1:1652818/49‑1 (MQ=255)
    gCTGGAAACCAGACCATTTTTACGCAGACGGGAGAACAGTTGTCCCAGAt                                            <  1:8137977/50‑1 (MQ=255)
    gCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGa                                             >  1:6213811/1‑49 (MQ=255)
    gCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                            <  1:2939656/50‑1 (MQ=255)
    gCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                            <  1:3307912/50‑1 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                            >  1:1491018/1‑49 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                            >  1:5785884/1‑49 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATa                                           <  1:7348412/50‑1 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATa                                           <  1:3959510/50‑1 (MQ=255)
     cTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGCCCCAGAt                                            >  1:6670097/1‑49 (MQ=39)
      tGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGAt                                            >  1:1756957/1‑48 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGa                                        <  1:4886056/50‑1 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATaa                                          <  1:6193539/48‑1 (MQ=255)
        gAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATaa                                          <  1:8147189/48‑1 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGaa                                       <  1:7294770/50‑1 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGaa                                       <  1:837742/50‑1 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGaa                                       <  1:8596708/50‑1 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGa                                        >  1:6960793/1‑49 (MQ=255)
         aaaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGa                                        >  1:3978993/1‑49 (MQ=255)
          aaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGa                                        >  1:8253074/1‑48 (MQ=255)
          aaCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGa                                        >  1:1242938/1‑48 (MQ=255)
           aCCAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAg                                     <  1:2592061/50‑1 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                    >  1:7981996/1‑50 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                    <  1:3278427/50‑1 (MQ=255)
            ccAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                    <  1:3802526/50‑1 (MQ=255)
             cAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAggg                                   <  1:1244147/50‑1 (MQ=255)
             cAGACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAgg                                    >  1:5277175/1‑49 (MQ=255)
              aGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAggg                                   >  1:5486938/1‑49 (MQ=255)
              aGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAggg                                   >  1:3825827/1‑49 (MQ=255)
              aGACCATTTGTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGa                                  <  1:3260966/50‑1 (MQ=255)
               gACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGaa                                 <  1:8278199/50‑1 (MQ=255)
               gACCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGa                                  >  1:5805661/1‑49 (MQ=255)
                aCCATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGaaa                                >  1:1564525/1‑50 (MQ=255)
                 ccATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAt                               <  1:3885865/50‑1 (MQ=255)
                 ccATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAt                               <  1:3379133/50‑1 (MQ=255)
                 ccATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAt                               >  1:4076846/1‑50 (MQ=255)
                  cATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAtt                              <  1:6573452/50‑1 (MQ=255)
                  cATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAtt                              >  1:1704623/1‑50 (MQ=255)
                  cATTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAAtt                              >  1:6791613/1‑50 (MQ=255)
                   aTTTTTACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTc                             >  1:507303/1‑50 (MQ=255)
                    tttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTcc                            >  1:5313101/1‑50 (MQ=255)
                     ttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                           >  1:1071393/1‑50 (MQ=255)
                     ttttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                           >  1:2879408/1‑50 (MQ=255)
                      tttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTccc                           <  1:2149442/49‑1 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCt                          >  1:2902451/1‑49 (MQ=255)
                       ttACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCt                          <  1:4853389/49‑1 (MQ=255)
                        tACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTg                         <  1:5636823/49‑1 (MQ=255)
                        tACGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGa                        >  1:1403009/1‑50 (MQ=255)
                         aCGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGAc                       >  1:1279226/1‑50 (MQ=255)
                          cGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGAc                       <  1:5840903/49‑1 (MQ=255)
                          cGCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGAc                       <  1:3805556/49‑1 (MQ=255)
                           gCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                     >  1:7084776/1‑50 (MQ=255)
                           gCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                     >  1:2587181/1‑50 (MQ=255)
                           gCAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGt                     >  1:2309269/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                    >  1:2898161/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                    >  1:3791735/1‑50 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                    <  1:5711759/50‑1 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                    <  1:3936827/50‑1 (MQ=255)
                            cAGACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGtt                    >  1:5198797/1‑50 (MQ=255)
                              gACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                  <  1:830777/50‑1 (MQ=255)
                              gACGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                  <  1:7432772/50‑1 (MQ=255)
                               aCGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg                 <  1:8167239/50‑1 (MQ=255)
                               aCGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCg                  >  1:6981235/1‑49 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCgg                 >  1:6551582/1‑49 (MQ=255)
                                cGGGAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa                <  1:1303092/50‑1 (MQ=255)
                                 gggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaa               <  1:2038651/50‑1 (MQ=255)
                                 gggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa                >  1:2850172/1‑49 (MQ=255)
                                 gggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGa                >  1:6008964/1‑49 (MQ=255)
                                  ggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa              <  1:5935289/50‑1 (MQ=255)
                                  ggAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGaaa              >  1:3871700/1‑50 (MQ=255)
                                   gAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAat             <  1:3844058/50‑1 (MQ=255)
                                   gAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAat             <  1:30647/50‑1 (MQ=255)
                                   gAAAACAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAat             <  1:545131/50‑1 (MQ=255)
                                    aaaaCAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata            >  1:1521605/1‑50 (MQ=255)
                                    aaaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata            >  1:7755525/1‑50 (MQ=255)
                                    aaaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata            <  1:1521553/50‑1 (MQ=255)
                                    aaaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata            >  1:1405520/1‑50 (MQ=255)
                                    aaaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAata            <  1:4227166/50‑1 (MQ=255)
                                     aaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat           <  1:7634865/50‑1 (MQ=255)
                                     aaaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat           >  1:1470269/1‑50 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat           >  1:5527252/1‑49 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat           >  1:4214716/1‑49 (MQ=255)
                                      aaCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATc          <  1:6309787/50‑1 (MQ=255)
                                       aCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAatat           >  1:4274531/1‑48 (MQ=255)
                                       aCAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa         <  1:3760574/50‑1 (MQ=255)
                                        cAGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAg        >  1:5977145/1‑50 (MQ=255)
                                         aGTTGTCCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGc       <  1:3684820/50‑1 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc      <  1:8038191/50‑1 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc      <  1:6450789/50‑1 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc      <  1:1578853/50‑1 (MQ=255)
                                          gttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc      <  1:8462396/50‑1 (MQ=255)
                                           ttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGcc      >  1:4211533/1‑49 (MQ=255)
                                           ttgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCa     <  1:6135789/50‑1 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa    >  1:6790940/1‑50 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa    >  1:5782451/1‑50 (MQ=255)
                                            tgtcCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCaa    >  1:5537837/1‑50 (MQ=255)
                                              ttCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACg  >  1:2391687/1‑50 (MQ=255)
                                               |                                                
GTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACG  >  NC_000913/2661957‑2662052

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