Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 823,380 G→A 41.2% P261L (CCG→CTG)  ybhO ← cardiolipin synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb823,3800GA41.2% 18.5 / 36.5 34P261L (CCG→CTG) ybhOcardiolipin synthase 2
Reads supporting (aligned to +/- strand):  ref base G (7/13);  new base A (6/8);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGC  >  W3110S.gb/823313‑823442
                                                                   |                                                              
cTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGc                                                              >  1:1503882/1‑70 (MQ=255)
     aaaCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATcc                                                        <  1:317065/71‑1 (MQ=255)
                 cgccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGa                                            <  1:3253767/71‑1 (MQ=255)
                 cgccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTg                                             <  1:2987529/70‑1 (MQ=255)
                 cgccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGa                                            <  1:771343/71‑1 (MQ=255)
                  gccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGaa                                           >  1:4010104/1‑71 (MQ=255)
                                aGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATc                             >  1:121968/1‑71 (MQ=255)
                                   tataCAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCg                           <  1:1024972/70‑1 (MQ=255)
                                   tataCAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCg                           >  1:3285093/1‑70 (MQ=255)
                                    ataCAGCAAGCGCGCACCGACTCTGACAATCGGCa                                                             >  1:157610/1‑35 (MQ=255)
                                    ataCAGCAAGCGCGCACCGACTCTGACAATCGGCa                                                             >  1:3881248/1‑35 (MQ=255)
                                    ataCAGCAAGCGCGCACCGACTCTGACAATCAGCa                                                             >  1:3809150/1‑35 (MQ=255)
                                     taCAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                <  1:3580936/63‑1 (MQ=255)
                                     taCAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                <  1:4033044/63‑1 (MQ=255)
                                         gCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCg                           <  1:764051/64‑1 (MQ=255)
                                         gCAAGCGCGCACCGACTATGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCg                           <  1:27799/64‑1 (MQ=255)
                                           aaGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCcgcg                  >  1:3018361/1‑71 (MQ=255)
                                           aaGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCcgcg                  >  1:1548518/1‑71 (MQ=255)
                                            aGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGaa                                           >  1:339924/1‑45 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                               <  1:725142/56‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGcc                <  1:2276860/71‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGcc                <  1:102798/71‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGcc                <  1:107628/71‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGcc                <  1:1647229/71‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                               <  1:2546259/56‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                               <  1:2065083/56‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                               <  1:4011769/56‑1 (MQ=255)
                                             gcgcgcACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                               <  1:1174020/56‑1 (MQ=255)
                                                    ccGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATg                                         <  1:2537046/39‑1 (MQ=255)
                                                    ccGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATg                                         <  1:2270419/39‑1 (MQ=255)
                                                     cGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCt        >  1:606050/1‑71 (MQ=255)
                                                      gACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGc         >  1:1631089/1‑69 (MQ=255)
                                                             aCAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGc  <  1:3637360/69‑1 (MQ=255)
                                                             aCAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCa                         >  1:521866/1‑46 (MQ=255)
                                                                   |                                                              
CTCAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGC  >  W3110S.gb/823313‑823442

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: