Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 823,380 G→A 21.6% 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,3800GA21.6% 74.5 / 16.1 37P261L (CCG→CTG) ybhOcardiolipin synthase 2
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (2/6);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGT  >  W3110S.gb/823315‑823449
                                                                 |                                                                     
cAAAAACCTGAACGCCGGCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCata                                                                  >  1:2410885/1‑71 (MQ=255)
  aaaaCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCata                                                                  <  1:1169642/69‑1 (MQ=255)
    aaCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCg                                                              >  1:2684488/1‑71 (MQ=255)
    aaCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCg                                                              >  1:824089/1‑71 (MQ=255)
     aCCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCgg                                                             <  1:2663142/71‑1 (MQ=255)
          aaCGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGc                                                        >  1:4072797/1‑71 (MQ=255)
              ccgccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATc                                                                >  1:1439644/1‑59 (MQ=255)
              ccgccTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTg                                                    <  1:1564394/71‑1 (MQ=255)
                     tAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGc                                                                     >  1:4196729/1‑47 (MQ=255)
                     tAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGc                                                        >  1:4461437/1‑60 (MQ=255)
                     tAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGAt                                              >  1:2025142/1‑70 (MQ=255)
                     tAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGAt                                              >  1:1161875/1‑70 (MQ=255)
                     tAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCACTGAATGAt                                              >  1:3210313/1‑70 (MQ=255)
                      aaCCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGa                                                   <  1:1401327/64‑1 (MQ=255)
                          aGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGttt                                        <  1:1232045/71‑1 (MQ=255)
                           gATAGTTATACAGCAAGCGCGCACCGACTCTGGCAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                       <  1:2293729/71‑1 (MQ=255)
                           gATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                       <  1:1426479/71‑1 (MQ=255)
                               gTTATACAGCAAGCGCGCACCGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATc                                    <  1:1521729/70‑1 (MQ=255)
                                      aGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGtt                                                           >  1:3607570/1‑40 (MQ=255)
                                        cAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCa                                            >  1:1034457/1‑53 (MQ=255)
                                        cAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCcgc                          >  1:2573184/1‑71 (MQ=255)
                                               gcACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGc                   <  1:3039675/71‑1 (MQ=255)
                                                   cGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATg                                                <  1:699104/38‑1 (MQ=255)
                                                   cGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGa                                               <  1:3088176/39‑1 (MQ=255)
                                                   cGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGa                                               <  1:4442449/39‑1 (MQ=255)
                                                   cGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGt                                          <  1:3892627/44‑1 (MQ=255)
                                                   cGACTCTGACAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCt               >  1:2707527/1‑71 (MQ=255)
                                                    gACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGAt                                              <  1:63537/39‑1 (MQ=255)
                                                    gACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCtt              <  1:2498898/71‑1 (MQ=255)
                                                     aCTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCtt              >  1:259420/1‑70 (MQ=255)
                                                      ctctGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCa                                <  1:1127162/51‑1 (MQ=255)
                                                       tctGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                       <  1:3272185/43‑1 (MQ=255)
                                                       tctGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTg                                       <  1:50711/43‑1 (MQ=255)
                                                           aCAATCGGCATATCCGGTTCCCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGc         >  1:1182055/1‑69 (MQ=255)
                                                           aCAATCAGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCa        <  1:4023156/70‑1 (MQ=255)
                                                              aTCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGa                                      <  1:4065405/37‑1 (MQ=255)
                                                              aTCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCAcc                              <  1:1535656/45‑1 (MQ=255)
                                                                 ggCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCgt  <  1:34168/70‑1 (MQ=255)
                                                                 |                                                                     
CAAAAACCTGAACGCCGCCTTTAACCAGATAGTTATACAGCAAGCGCGCACCGACTCTGACAATCGGCATATCCGGTTCGCCCTGAATGATCAGTTTGATCCGCACCCCGCGCCGTGCCGCTTTACGCAAGGCGT  >  W3110S.gb/823315‑823449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: