Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 453,395:1 +G 100% intergenic (‑4/‑301) yajG ← / → bolA predicted lipoprotein/regulator of penicillin binding proteins and beta lactamase transcription (morphogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb453,3951.G80.0% 35.4 / 4.1 15intergenic (‑4/‑301)yajG/bolApredicted lipoprotein/regulator of penicillin binding proteins and beta lactamase transcription (morphogene)
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (6/6);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGCATCCTGCAAGCATAAACAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTCCCGCCATCATAACATTGCCAA  >  W3110S.gb/453337‑453455
                                                           |                                                            
cgcATCCTGCAAGCATAAACAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCag                                                    <  1:36425/70‑1 (MQ=255)
              aTAAACAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCgtg                                     >  1:472668/1‑71 (MQ=255)
               tAAACAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCgtgt                                    <  1:1517549/71‑1 (MQ=255)
                 aaCAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGttt                                  >  1:491667/1‑71 (MQ=255)
                 aaCAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGttt                                  >  1:64454/1‑71 (MQ=255)
                          aCTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAAt                         <  1:174154/71‑1 (MQ=255)
                          aCTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAAt                         <  1:665191/71‑1 (MQ=255)
                          aCTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAAtt                        <  1:1426835/71‑1 (MQ=255)
                           cTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAAtt                        >  1:1344956/1‑71 (MQ=255)
                             aaCGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTcc                      >  1:476781/1‑71 (MQ=255)
                              aCGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTcc                      >  1:721344/1‑70 (MQ=255)
                              aCGGGAAGAGGATTTTTTTGAACATGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTc                       <  1:1656541/69‑1 (MQ=255)
                                                  aaCATGTTCGGGCTCTCAGAGACTCTTAAGCGTGtt                                   <  1:121045/36‑1 (MQ=255)
                                                     aTGTTCGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTCCCGCCATCATAACATTGCCaa  >  1:768667/1‑66 (MQ=255)
                                                      tGTTCGGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTCCCGCCATCATAACa         <  1:641788/59‑1 (MQ=255)
                                                           |                                                            
CGCATCCTGCAAGCATAAACAGAGCAACTAACGGGAAGAGGATTTTTTTGAACATGTTCGGGCTCTCAGAGACTCTTAAGCGTGTTTGGTAAAAATTCCCGCCATCATAACATTGCCAA  >  W3110S.gb/453337‑453455

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: