Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 711,264 A→C 10.7% intergenic (‑196/+93) 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,2640AC10.7% 240.0 / 32.2 103intergenic (‑196/+93)fur/fldADNA‑binding transcriptional dual regulator/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base A (46/46);  new base C (7/4);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 103 total aligned reads displayed.

AATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTA  >  W3110S.gb/711200‑711332
                                                                |                                                                    
aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                                 <  1:1572130/70‑1 (MQ=255)
aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                                 <  1:3994131/70‑1 (MQ=255)
 ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTac                                                                 <  1:2572018/69‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:1149516/70‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:4315501/70‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:5378164/70‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:3528954/70‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:2725383/70‑1 (MQ=255)
  taaAAACGTGGCACAGCTTACATTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               <  1:924395/69‑1 (MQ=255)
    aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               >  1:692210/1‑68 (MQ=255)
    aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                               >  1:236121/1‑68 (MQ=255)
     aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAg                                                            >  1:760127/1‑70 (MQ=255)
     aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAg                                                            >  1:3830079/1‑70 (MQ=255)
       aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:4727347/71‑1 (MQ=255)
       aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:5139644/71‑1 (MQ=255)
       aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:4774956/71‑1 (MQ=255)
       aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:1472123/71‑1 (MQ=255)
       aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:5482477/71‑1 (MQ=255)
       aCGTGGCAAAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         <  1:3286070/71‑1 (MQ=255)
        cGTGGCACAGCTTACATTTACAACGGTAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         >  1:4908574/1‑70 (MQ=255)
        cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         >  1:5596218/1‑70 (MQ=255)
        cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         >  1:5339302/1‑70 (MQ=255)
        cGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCaa                                                         >  1:1932490/1‑70 (MQ=255)
                aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                                 >  1:4229059/1‑70 (MQ=255)
                aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                                 >  1:2481526/1‑70 (MQ=255)
                aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                                 >  1:4983198/1‑70 (MQ=255)
                aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                                 >  1:3736555/1‑70 (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:3219146/71‑1 (MQ=255)
                                  cAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                              <  1:5280959/71‑1 (MQ=255)
                                   aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:4484858/71‑1 (MQ=255)
                                   aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:3839149/71‑1 (MQ=255)
                                   aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:659963/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:3155109/70‑1 (MQ=255)
                                     ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                            <  1:625759/70‑1 (MQ=255)
                                     ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                            <  1:2355029/70‑1 (MQ=255)
                                     ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                            <  1:2378761/70‑1 (MQ=255)
                                       ccACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                            <  1:4548456/68‑1 (MQ=255)
                                         aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                        >  1:4818108/1‑70 (MQ=255)
                                         aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                        >  1:918314/1‑70 (MQ=255)
                                         aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                        >  1:1935389/1‑70 (MQ=255)
                                         aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                        >  1:3817482/1‑70 (MQ=255)
                                         aCTTGAACACCGGGAATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCa                      >  1:1064471/1‑71 (MQ=255)
                                           ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:455759/63‑1 (MQ=255)
                                           ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt                     >  1:4097117/1‑71 (MQ=255)
                                           ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt                     >  1:721986/1‑71 (MQ=255)
                                            tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt                     >  1:3922649/1‑70 (MQ=255)
                                            tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCAt                     >  1:3649912/1‑70 (MQ=255)
                                             gAACACCGGGACATTGTACAGGTACAACAGCAAAAGt                                                     >  1:3709991/1‑37 (MQ=255)
                                             gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:3745024/1‑63 (MQ=255)
                                             gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:375850/1‑63 (MQ=255)
                                             gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:2606219/1‑63 (MQ=255)
                                             gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:923944/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:5312979/71‑1 (MQ=255)
                                                     ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa           <  1:5016123/71‑1 (MQ=255)
                                                      taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa           <  1:4349092/69‑1 (MQ=255)
                                                      taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa           <  1:2040489/69‑1 (MQ=255)
                                                      taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa           <  1:2566901/69‑1 (MQ=255)
                                                      taCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTaa           <  1:56839/69‑1 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAAt          >  1:5298537/1‑70 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:81127/1‑71 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:950280/1‑71 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:4785368/1‑71 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:1262089/1‑71 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:1495799/1‑71 (MQ=255)
                                                       aCAATGTCCCGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATa         >  1:1603550/1‑71 (MQ=255)
                                                         aTTGTACAGGTACAACAGCGAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:2517312/49‑1 (MQ=255)
                                                         aTTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                              <  1:4603933/48‑1 (MQ=255)
                                                         aTTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttt                              <  1:4586683/48‑1 (MQ=255)
                                                           tGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAAt          >  1:506884/1‑66 (MQ=255)
                                                           tGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGc     >  1:114872/1‑71 (MQ=255)
                                                           tGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGc     >  1:181634/1‑71 (MQ=255)
                                                           tGTACAGGTACAAAAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAAt          >  1:2200284/1‑66 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:2694053/1‑71 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:4852599/1‑71 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:1904993/1‑71 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:83097/1‑71 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:4331150/1‑71 (MQ=255)
                                                            gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:4457220/1‑71 (MQ=255)
                                                            gTACAGGTAAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    >  1:3700917/1‑70 (MQ=255)
                                                             tACAGGTACAACAGCAAAAGTTACAAATTTGTATCAATTATTTTGATTGGCATTATCTATTAATACGGCg    <  1:1322479/70‑1 (MQ=255)
                                                             tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                             <  1:2115294/45‑1 (MQ=255)
                                                             tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg    <  1:3333911/70‑1 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:3249782/1‑45 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:2816068/1‑45 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:2049521/1‑45 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                           >  1:4251249/1‑45 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTa  <  1:2868617/70‑1 (MQ=255)
                                                               cAGGTACAACAGCAAAAGTTACAAAATTGTAGCAATTATTTTGAt                           >  1:1930007/1‑45 (MQ=255)
                                                                |                                                                    
AATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTA  >  W3110S.gb/711200‑711332

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: