Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,084,437 Δ1 bp coding (131/1977 nt) speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,084,4340C.100.0% 132.6 / NA 31coding (134/1977 nt)speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/18);  total (13/18)

CGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAT  >  W3110S.gb/3084375‑3084489
                                                           |                                                       
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       >  1:2970152/1‑62 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       <  1:996905/62‑1 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       >  1:555829/1‑62 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       <  1:3312689/62‑1 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       <  1:2939510/62‑1 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       <  1:2270286/62‑1 (MQ=255)
cGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       >  1:198928/1‑62 (MQ=255)
 gtccgggtccggGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGccc                                                       <  1:749834/61‑1 (MQ=255)
           gggCACACGCTAATGTGGCCCCGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAata                                           >  1:2027652/1‑62 (MQ=255)
           gggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCaa                                             >  1:736486/1‑60 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:1650514/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:898363/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:672813/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:3619728/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:3215242/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:3139629/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:216581/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:2675515/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:2390941/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:1883764/62‑1 (MQ=255)
                  cGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAg                                    <  1:2286379/62‑1 (MQ=255)
                         gtgGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACgcag                             >  1:143596/1‑62 (MQ=255)
                         gtgGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACgcag                             >  1:2350272/1‑62 (MQ=255)
                         gtgGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACgca                              >  1:808647/1‑61 (MQ=255)
                              ccAGCTCGTTAACGTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCt                        >  1:2881803/1‑62 (MQ=255)
                                           gTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCt            <  1:3438883/61‑1 (MQ=255)
                                           gTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCt            <  1:4219112/61‑1 (MQ=255)
                                           gTCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTg           <  1:1635263/62‑1 (MQ=255)
                                            tCATAGTAGTTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCt            <  1:2584820/60‑1 (MQ=255)
                                                 gtagtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGcg     <  1:4245991/62‑2 (MQ=255)
                                                 gtagtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCt     <  1:3887797/62‑1 (MQ=255)
                                                 gtagtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCt     <  1:1715728/62‑1 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCa   >  1:2718178/1‑61 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:459790/1‑62 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:2321157/1‑62 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:3529135/1‑62 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:788352/1‑62 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:3523478/1‑62 (MQ=255)
                                                    gtTATTGCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAt  >  1:1442158/1‑62 (MQ=255)
                                                           |                                                       
CGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCAT  >  W3110S.gb/3084375‑3084489

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