Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 711,258 T→A 11.7% intergenic (‑190/+99) fur ← / ← fldA DNA‑binding transcriptional dual regulator/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb711,2580TA11.7% 257.7 / 23.8 94intergenic (‑190/+99)fur/fldADNA‑binding transcriptional dual regulator/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base T (34/49);  new base A (7/4);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATA  >  W3110S.gb/711190‑711325
                                                                    |                                                                   
aaaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:4184720/71‑1 (MQ=255)
aaaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:4317089/71‑1 (MQ=255)
aaaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:4498650/71‑1 (MQ=255)
  atatTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 <  1:1368439/71‑1 (MQ=255)
     ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 >  1:3884160/1‑68 (MQ=255)
     ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 >  1:3297084/1‑68 (MQ=255)
     ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 >  1:4341769/1‑68 (MQ=255)
     ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 >  1:2000678/1‑68 (MQ=255)
     ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                 >  1:4836106/1‑68 (MQ=255)
          aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                          <  1:3994131/70‑1 (MQ=255)
          aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                          <  1:1572130/70‑1 (MQ=255)
           ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                          <  1:2572018/69‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:1149516/70‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:3528954/70‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:4315501/70‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:2725383/70‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:5378164/70‑1 (MQ=255)
            taaAAACGTGGCACAGCTTACATTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        <  1:924395/69‑1 (MQ=255)
              aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        >  1:236121/1‑68 (MQ=255)
              aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        >  1:692210/1‑68 (MQ=255)
              aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                        >  1:295664/1‑68 (MQ=255)
               aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAg                                                     >  1:3830079/1‑70 (MQ=255)
               aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAg                                                     >  1:760127/1‑70 (MQ=255)
                 aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:5482477/71‑1 (MQ=255)
                 aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:1472123/71‑1 (MQ=255)
                 aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:5139644/71‑1 (MQ=255)
                 aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:4727347/71‑1 (MQ=255)
                 aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:4774956/71‑1 (MQ=255)
                 aCGTGGCAAAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  <  1:3286070/71‑1 (MQ=255)
                  cGTGGCACAGCTTACATTTACAACGGTAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  >  1:4908574/1‑70 (MQ=255)
                  cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  >  1:1932490/1‑70 (MQ=255)
                  cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  >  1:5596218/1‑70 (MQ=255)
                  cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                  >  1:5339302/1‑70 (MQ=255)
                          aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                          >  1:4983198/1‑70 (MQ=255)
                          aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                          >  1:4229059/1‑70 (MQ=255)
                          aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                          >  1:979252/1‑70 (MQ=255)
                          aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                          >  1:3736555/1‑70 (MQ=255)
                          aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                          >  1:2481526/1‑70 (MQ=255)
                           gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:3928046/44‑1 (MQ=255)
                           gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:4070649/44‑1 (MQ=255)
                           gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTg                                                                   <  1:4557991/44‑1 (MQ=255)
                                     aCAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCaa                              >  1:482686/1‑71 (MQ=255)
                                     aCAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCa                               <  1:694639/70‑1 (MQ=255)
                                     aCAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCa                               <  1:1911335/70‑1 (MQ=255)
                                          ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAattat                          <  1:3571877/70‑1 (MQ=255)
                                            cAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                       <  1:5280959/71‑1 (MQ=255)
                                            cAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                       <  1:3219146/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:3839149/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:659963/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:4484858/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:3149791/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:1685306/71‑1 (MQ=255)
                                             aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGAAGCAATTATTTTg                      <  1:2487198/71‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAat                             <  1:316835/62‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAat                             <  1:700203/62‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                     <  1:2378761/70‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                     <  1:3155109/70‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                     <  1:2355029/70‑1 (MQ=255)
                                               ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                     <  1:625759/70‑1 (MQ=255)
                                                 ccACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                     <  1:4548456/68‑1 (MQ=255)
                                                   aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                 >  1:1935389/1‑70 (MQ=255)
                                                   aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                 >  1:3817482/1‑70 (MQ=255)
                                                   aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                 >  1:918314/1‑70 (MQ=255)
                                                   aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                 >  1:4818108/1‑70 (MQ=255)
                                                   aCTTGAACACCGGGAATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCa               >  1:1064471/1‑71 (MQ=255)
                                                     ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:455759/63‑1 (MQ=255)
                                                     ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt              >  1:721986/1‑71 (MQ=255)
                                                     ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt              >  1:4097117/1‑71 (MQ=255)
                                                      tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt              >  1:3649912/1‑70 (MQ=255)
                                                      tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt              >  1:3922649/1‑70 (MQ=255)
                                                       gAACACCGGGACATTGTACAGGTACAACAGCAAAAGt                                              >  1:3709991/1‑37 (MQ=255)
                                                       gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                    >  1:3745024/1‑63 (MQ=255)
                                                       gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                    >  1:923944/1‑63 (MQ=255)
                                                       gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                    >  1:375850/1‑63 (MQ=255)
                                                       gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                    >  1:2606219/1‑63 (MQ=255)
                                                         acacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATc          <  1:4410617/71‑1 (MQ=255)
                                                           acCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTa        <  1:789155/71‑1 (MQ=255)
                                                               ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:2866655/53‑1 (MQ=255)
                                                               ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:5016123/71‑1 (MQ=255)
                                                               ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:5312979/71‑1 (MQ=255)
                                                                taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:56839/69‑1 (MQ=255)
                                                                taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:4349092/69‑1 (MQ=255)
                                                                taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:2040489/69‑1 (MQ=255)
                                                                taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa    <  1:2566901/69‑1 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAAt   >  1:5298537/1‑70 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:4785368/1‑71 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:1603550/1‑71 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:1495799/1‑71 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:81127/1‑71 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:1262089/1‑71 (MQ=255)
                                                                 aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa  >  1:950280/1‑71 (MQ=255)
                                                                   aTTGTACAGGTACAACAGCGAAAGTTACAAATTTGTAGCAATTATTTTg                      <  1:2517312/49‑1 (MQ=255)
                                                                   aTTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                       <  1:4586683/48‑1 (MQ=255)
                                                                   aTTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                       <  1:4603933/48‑1 (MQ=255)
                                                                    |                                                                   
AAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATA  >  W3110S.gb/711190‑711325

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