Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 478,343 C→T 47.4% A45T (GCT→ACT)  ylaC ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb478,3430CT47.4% 2.4 / 51.7 38A45T (GCT→ACT) ylaCpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (6/14);  new base T (5/13);  total (11/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

CAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTTGTTGTC  >  W3110S.gb/478274‑478412
                                                                     |                                                                     
cAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGc                                                                       <  1:1197904/70‑1 (MQ=255)
 aaCCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGc                                                                       <  1:3757719/69‑1 (MQ=255)
   ccAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACg                                                                   <  1:2389856/71‑1 (MQ=255)
    cAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGt                                                                  <  1:355067/71‑1 (MQ=255)
     aGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTa                                                                 <  1:205817/71‑1 (MQ=255)
     aGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTa                                                                 <  1:4822172/71‑1 (MQ=255)
        aGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                               <  1:2026268/70‑1 (MQ=255)
               cagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGAt                                                        >  1:759372/1‑70 (MQ=255)
                             gTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCGACGTa                                                                 <  1:3746817/47‑1 (MQ=255)
                                tCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGtt                                       >  1:4143341/1‑70 (MQ=255)
                                      tGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATa                                                       <  1:3106801/48‑1 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTAAATACCGATAAACAGCCACGGATGTTTACGGAt                                <  1:337971/69‑1 (MQ=255)
                                          aCATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCgg                                            <  1:2109662/55‑1 (MQ=255)
                                          aCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                            <  1:2737331/71‑1 (MQ=255)
                                          aCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaa                             >  1:3919325/1‑70 (MQ=255)
                                           cATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                               <  1:3340856/35‑1 (MQ=255)
                                             tCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACt                          >  1:4476065/1‑70 (MQ=255)
                                                     aGGGTGGAAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGtt                                       <  1:5028612/49‑1 (MQ=255)
                                                       ggTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGc               <  1:4873021/71‑1 (MQ=255)
                                                       ggTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGc               <  1:1112600/71‑1 (MQ=255)
                                                       ggTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGc               <  1:1228430/71‑1 (MQ=255)
                                                         tGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCgg             >  1:1326144/1‑71 (MQ=255)
                                                           gCAAAAAAAGTAACGTACATACCGATAAACAGCCACGGATGTTTACGGATAAAACTGATACTAAAGCgggg           >  1:2008627/1‑71 (MQ=255)
                                                           gCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                            >  1:2786818/1‑54 (MQ=255)
                                                             aaaaaaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGa                        >  1:2846488/1‑56 (MQ=255)
                                                             aaaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGtt         <  1:3930007/71‑1 (MQ=255)
                                                              aaaaaaGTAACGTACATACCGATAAACAGCCCCGGATGtt                                       <  1:3415980/40‑1 (MQ=255)
                                                              aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                            >  1:990323/1‑51 (MQ=255)
                                                              aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGAt                       >  1:3818881/1‑56 (MQ=255)
                                                              aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGttt        <  1:706424/71‑1 (MQ=255)
                                                              aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGtt         <  1:692527/70‑1 (MQ=255)
                                                               aaaaaGTAACGTACATACCGATAAACAGCCCCGGATGtt                                       <  1:1389420/39‑1 (MQ=255)
                                                               aaaaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg       <  1:3046823/71‑1 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg       <  1:3817972/71‑1 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg       <  1:1106562/71‑1 (MQ=255)
                                                                 aaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgtt     <  1:2018263/71‑1 (MQ=255)
                                                                 aaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgtt     <  1:3502288/71‑1 (MQ=255)
                                                                  aaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                            >  1:2579406/1‑47 (MQ=255)
                                                                   aGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGTGTTtgttgt   <  1:3620127/71‑1 (MQ=255)
                                                                    gCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTTGTTGTc  <  1:2036790/71‑1 (MQ=255)
                                                                     |                                                                     
CAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTTGTTGTC  >  W3110S.gb/478274‑478412

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: