Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 711,262 A→C 11.3% intergenic (‑194/+95) 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,2620AC11.3% 219.2 / 32.2 97intergenic (‑194/+95)fur/fldADNA‑binding transcriptional dual regulator/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base A (41/45);  new base C (7/4);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 103 total aligned reads displayed.

ATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCG  >  W3110S.gb/711192‑711330
                                                                      |                                                                    
atatTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      <  1:1368439/71‑1 (MQ=255)
   ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      >  1:4836106/1‑68 (MQ=255)
   ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      >  1:3297084/1‑68 (MQ=255)
   ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      >  1:2000678/1‑68 (MQ=255)
   ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      >  1:3884160/1‑68 (MQ=255)
   ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTa                                                                      >  1:4341769/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:5378164/70‑1 (MQ=255)
          taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             <  1:4315501/70‑1 (MQ=255)
          taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             <  1:3528954/70‑1 (MQ=255)
          taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             <  1:1149516/70‑1 (MQ=255)
          taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             <  1:2725383/70‑1 (MQ=255)
          taaAAACGTGGCACAGCTTACATTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             <  1:924395/69‑1 (MQ=255)
            aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTacaa                                                             >  1:295664/1‑68 (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:5482477/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)
               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:3736555/1‑70 (MQ=255)
                        aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                               >  1:4229059/1‑70 (MQ=255)
                        aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                               >  1:4983198/1‑70 (MQ=255)
                        aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACa                                               >  1:2481526/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:3149791/71‑1 (MQ=255)
                                           aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                           <  1:1685306/71‑1 (MQ=255)
                                           aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                           <  1:659963/71‑1 (MQ=255)
                                           aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTg                           <  1:3839149/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:625759/70‑1 (MQ=255)
                                             ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                          <  1:2355029/70‑1 (MQ=255)
                                               ccACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGa                          <  1:4548456/68‑1 (MQ=255)
                                                 aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                      >  1:918314/1‑70 (MQ=255)
                                                 aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTgg                      >  1:4818108/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:2606219/1‑63 (MQ=255)
                                                     gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                         >  1:923944/1‑63 (MQ=255)
                                                     gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                         >  1:375850/1‑63 (MQ=255)
                                                     gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAt                         >  1:3745024/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:56839/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: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)
                                                                   tGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAAt        >  1:4616888/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: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)
                                                                    gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCg  >  1:4852599/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)
                                                                      |                                                                    
ATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCG  >  W3110S.gb/711192‑711330

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: