Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,477,260 G→A 54.1% intergenic (‑45/+236) pyrL ← / ← yjgH pyrBI operon leader peptide/predicted mRNA endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,477,2600GA54.1% ‑1.2 / 27.9 23intergenic (‑45/+236)pyrL/yjgHpyrBI operon leader peptide/predicted mRNA endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base A (4/8);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGT  >  W3110S.gb/4477197‑4477325
                                                               |                                                                 
gTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGgcgc                                                             >  1:53119/1‑70 (MQ=255)
 tCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCa                                                            <  1:964811/70‑1 (MQ=255)
          cTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAacggcg                                                              <  1:252305/59‑1 (MQ=255)
                  tACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAt                                                       <  1:1145297/58‑1 (MQ=255)
                  tACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATAtt                                          <  1:1624025/71‑1 (MQ=255)
                      tcctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGt                                      <  1:883204/71‑1 (MQ=255)
                              gCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCaa                              <  1:1518938/71‑1 (MQ=255)
                               aaaaTTGTCCGGCATTATACTCATCGTCAGAATCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCaaa                             <  1:506030/70‑1 (MQ=255)
                                aaaTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGa                                                   >  1:1161768/1‑48 (MQ=255)
                                aaaTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAAc                            >  1:401591/1‑71 (MQ=255)
                                        cGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAt                     <  1:194923/70‑1 (MQ=255)
                                             ttATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGAGCAAACGGTTTAttt                   <  1:904934/67‑1 (MQ=255)
                                             ttATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGGCATCGCAAACGGTTTAttt                   <  1:4053/67‑1 (MQ=255)
                                               ataCTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAAt               >  1:1305610/1‑69 (MQ=255)
                                                taCTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaa            <  1:963273/71‑1 (MQ=255)
                                                   tCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                               <  1:1761905/49‑1 (MQ=255)
                                                   tCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                               <  1:1052517/49‑1 (MQ=255)
                                                   tCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAg         >  1:439541/1‑71 (MQ=255)
                                                     aTCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                          >  1:1044644/1‑52 (MQ=255)
                                                      tCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTa             <  1:790682/64‑1 (MQ=255)
                                                           aGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGt  <  1:1374395/70‑1 (MQ=255)
                                                            gAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaaaa          >  1:471014/1‑61 (MQ=255)
                                                            gAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTa             <  1:1124429/58‑1 (MQ=255)
                                                            gAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTa             <  1:1079863/58‑1 (MQ=255)
                                                               |                                                                 
GTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGT  >  W3110S.gb/4477197‑4477325

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: