Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,833,686 G→A 21.9% P4L (CCA→CTA)  astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,833,6860GA21.9% 80.8 / 18.6 41P4L (CCA→CTA) astCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (12/20);  new base A (3/6);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAG  >  W3110S.gb/1833623‑1833743
                                                               |                                                         
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGTCTgaga                                                    <  1:916734/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:1226653/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:2747298/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:2264564/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:2058786/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:1672395/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:1561594/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:1224669/71‑1 (MQ=255)
ggTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:969113/71‑1 (MQ=255)
 gTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTgaga                                                    <  1:633014/70‑1 (MQ=255)
               ggAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTg                                                           >  1:2280771/1‑49 (MQ=255)
               ggAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATa                                                >  1:80616/1‑60 (MQ=255)
               ggAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATa                                                >  1:2483351/1‑60 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:1825322/64‑1 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:700408/64‑1 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:534957/64‑1 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:1766655/64‑1 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:1751181/64‑1 (MQ=255)
                                    cATTCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGt                       <  1:2988067/64‑1 (MQ=255)
                                      ttCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGCCCTCTACAGTAAATAATTCGtta                     >  1:3304188/1‑64 (MQ=255)
                                      ttCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGtta                     >  1:951443/1‑64 (MQ=255)
                                      ttCATCAAAGTTTTCACGCGTAATTAGCTGAGACATAGCGACCTCTACAGTAAATAATTCGtta                     >  1:2952828/1‑64 (MQ=255)
                                            aaaGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:1793060/1‑70 (MQ=255)
                                            aaaGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:691652/1‑70 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:989195/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:46318/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:717622/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:1473470/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:2512470/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:2408119/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:1982893/1‑71 (MQ=255)
                                                  ttCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAg  >  1:1984059/1‑71 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:2527669/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:1586236/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:1820405/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:795752/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:1829983/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:3283534/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:3260450/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:3010769/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:2959597/53‑1 (MQ=255)
                                                          tAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGt            <  1:2622498/53‑1 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:218476/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:325409/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:1982728/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:194038/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:1915327/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:3161538/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:3013922/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:1598212/1‑51 (MQ=255)
                                                               ggCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaa         >  1:1327218/1‑51 (MQ=255)
                                                               |                                                         
GGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAATAAG  >  W3110S.gb/1833623‑1833743

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: