Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,225,730 G→A 31.1% T117I (ACT→ATT)  arpA ← regulator of acetyl CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,225,7300GA31.1% 60.1 / 35.5 45T117I (ACT→ATT) arpAregulator of acetyl CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base A (4/10);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TAAGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGTA  >  W3110S.gb/4225662‑4225795
                                                                    |                                                                 
tAAGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAgt                                                                  <  1:628380/70‑1 (MQ=255)
tAAGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAatg                                                                 >  1:185512/1‑71 (MQ=255)
 aaGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAgtg                                                                 >  1:2544938/1‑70 (MQ=255)
  aGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTc                                                               >  1:2615296/1‑71 (MQ=255)
        aTTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGc                                                          <  1:329055/70‑1 (MQ=255)
        aTTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCg                                                         <  1:3371002/71‑1 (MQ=255)
         ttCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGt                                                        >  1:290576/1‑71 (MQ=255)
           cTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGtt                                                       >  1:803324/1‑70 (MQ=255)
              gggAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGa                                                     <  1:3794895/69‑1 (MQ=255)
              gggAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGa                                                     <  1:363531/69‑1 (MQ=255)
                 aaTTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAgg                                                   >  1:308714/1‑68 (MQ=255)
                 aaTTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTa                                                 >  1:2757958/1‑70 (MQ=255)
                   ttctttttctttTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTaa                                                <  1:312161/69‑1 (MQ=255)
                     ctttttcttttCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAAtt                                            >  1:2038175/1‑71 (MQ=255)
                     ctttttcttttCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAAtt                                            >  1:151417/1‑71 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAAtt                                            <  1:858914/68‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTc                                          <  1:2979715/70‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTc                                          <  1:3536959/70‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTc                                          <  1:3677148/70‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTc                                          <  1:152103/70‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAAtt                                            <  1:3088856/68‑1 (MQ=255)
                        tttcttttctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAATTTc                                          <  1:3214364/70‑1 (MQ=255)
                               tctATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAgg                                   <  1:1352406/70‑1 (MQ=255)
                                   tCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGGTGAGGTAAAATTTCCATTAGGTAt                                <  1:3778798/69‑1 (MQ=255)
                                     tGATTTTAGGTTTAATTCTTGTTTATCAAGAgtgt                                                                >  1:3951455/1‑35 (MQ=255)
                                     tGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTAtttt                             >  1:1599208/1‑70 (MQ=255)
                                      gATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCAt                                       >  1:2820165/1‑59 (MQ=255)
                                       aTTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTc                           >  1:2693310/1‑70 (MQ=255)
                                       aTTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTc                           >  1:3880123/1‑70 (MQ=255)
                                       aTTTTAGGTTTAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCa                          <  1:666772/71‑1 (MQ=255)
                                          ttAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCaaa                        <  1:2475838/70‑1 (MQ=255)
                                          ttAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCaaa                        <  1:3548051/70‑1 (MQ=255)
                                             ggTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGTTATTTTTCAAatat                     <  1:1157355/70‑1 (MQ=255)
                                                ttAATTCTTGTTTATCAAGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTggg                 >  1:2464148/1‑71 (MQ=255)
                                                    ttCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTaaa                                               <  1:908181/37‑1 (MQ=255)
                                                    ttCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTaaa                                               <  1:3544758/37‑1 (MQ=255)
                                                    ttCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTaaa                                               <  1:3215793/37‑1 (MQ=255)
                                                      cTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAAt                                             >  1:3621997/1‑37 (MQ=255)
                                                      cTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAAt                                             >  1:2952638/1‑37 (MQ=255)
                                                      cTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCa                                        >  1:2708261/1‑42 (MQ=255)
                                                      cTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATg           >  1:3058501/1‑71 (MQ=255)
                                                      cTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATg           >  1:1842285/1‑71 (MQ=255)
                                                          tttATCAAGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGttt       <  1:1929218/71‑1 (MQ=255)
                                                              tCAAGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGt   <  1:3073521/71‑1 (MQ=255)
                                                               cAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGc               <  1:1691903/58‑1 (MQ=255)
                                                                aaGAGTGTCTGAGCGTTGAGTTAAAATTTCCATTAGGTATTTTTCAAATATTgggtc               <  1:3497958/57‑3 (MQ=255)
                                                                aaGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGTa  <  1:3532523/70‑1 (MQ=255)
                                                                aaGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGTa  <  1:2515142/70‑1 (MQ=255)
                                                                aaGAATGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGTa  <  1:1127229/70‑1 (MQ=255)
                                                                    |                                                                 
TAAGTTTAATTCTTGGGAATTCTTTTTCTTTTCTATCTGATTTTAGGTTTAATTCTTGTTTATCAAGAGTGTCTGAGCGTTGAGGTAAAATTTCCATTAGGTATTTTTCAAATATTGGGGCAATGGTTTGAGTA  >  W3110S.gb/4225662‑4225795

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: