Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,434,349 T→A 100% intergenic (+112/‑197) fklB → / → cycA FKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,434,3490TA84.8% 76.5 / 8.6 33intergenic (+112/‑197)fklB/cycAFKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base A (16/12);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAT  >  W3110S.gb/4434289‑4434410
                                                            |                                                             
aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                     >  1:2334574/1‑71 (MQ=255)
aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                     >  1:2153581/1‑71 (MQ=255)
 aCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGt                                                    <  1:1293384/71‑1 (MQ=255)
         gTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGtt                                            <  1:1478354/71‑1 (MQ=255)
           ttGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                          >  1:1145807/1‑71 (MQ=255)
           ttGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                          >  1:805661/1‑71 (MQ=255)
            tGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTAt                                          <  1:2671359/70‑1 (MQ=255)
                 aaCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGt                                    <  1:2008613/71‑1 (MQ=255)
                   cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                  <  1:1165740/71‑1 (MQ=255)
                        cTATTTTGTGTTGCAAAGACCCCGTAATCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                             <  1:78610/71‑1 (MQ=255)
                           ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGtt                                            <  1:479200/53‑1 (MQ=255)
                           ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                  <  1:841568/63‑1 (MQ=255)
                            tttGTGTTCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                         >  1:797347/1‑70 (MQ=255)
                             ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                          >  1:2593026/1‑53 (MQ=255)
                             ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                         >  1:2142139/1‑70 (MQ=255)
                             ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                         >  1:230582/1‑70 (MQ=255)
                             ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                         <  1:818147/70‑1 (MQ=255)
                             ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                         >  1:1223368/1‑70 (MQ=255)
                                  ttGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                             <  1:987396/45‑1 (MQ=255)
                                       aaGACCCCGTAAGCGTGTATTATTGTGAGCTTTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACa              <  1:561080/71‑1 (MQ=255)
                                          aCCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa           <  1:1691974/71‑1 (MQ=255)
                                            cccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATc                                         >  1:2436367/1‑39 (MQ=255)
                                            cccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATc                                         >  1:2180297/1‑39 (MQ=255)
                                            cccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATc                                         >  1:1847933/1‑39 (MQ=255)
                                                  aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:2474049/1‑71 (MQ=255)
                                                  aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:487786/1‑71 (MQ=255)
                                                  aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:1437760/1‑71 (MQ=255)
                                                  aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:1144009/1‑71 (MQ=255)
                                                  aGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:427095/1‑71 (MQ=255)
                                                  aGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:687752/1‑71 (MQ=255)
                                                  aGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:1716569/1‑71 (MQ=255)
                                                  aGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACa   >  1:1379942/1‑71 (MQ=255)
                                                  aGCGTGTATATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAt  <  1:2389722/71‑1 (MQ=255)
                                                        tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAt  >  1:2673495/1‑66 (MQ=255)
                                                        tattatTGTGAGCTGTTTCGCGTTATCACAGTGATATGACACTCACTTTAAACATAAAATTAACAt  >  1:2456617/1‑66 (MQ=255)
                                                            |                                                             
AACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAT  >  W3110S.gb/4434289‑4434410

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: