Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 478,343 C→T 47.6% 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.6% ‑2.3 / 26.8 21A45T (GCT→ACT) ylaCpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/8);  new base T (5/5);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGT  >  W3110S.gb/478277‑478404
                                                                  |                                                             
ccAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACg                                                           <  1:1390614/71‑1 (MQ=255)
ccAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACg                                                           >  1:1613986/1‑71 (MQ=255)
     aGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACa                                                       >  1:494908/1‑70 (MQ=255)
      gCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACAt                                                      <  1:384305/70‑1 (MQ=255)
       ccAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATAc                                                    <  1:1317336/71‑1 (MQ=255)
           acagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGt                                                          <  1:1082396/61‑1 (MQ=255)
           acagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGa                                                 >  1:1544521/1‑70 (MQ=255)
           acagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGAt                                                <  1:843656/71‑1 (MQ=255)
            cagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATa                                               <  1:904085/71‑1 (MQ=255)
                               ggACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAAc                                                            >  1:719740/1‑39 (MQ=255)
                                     cAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATa                       >  1:1741016/1‑70 (MQ=255)
                                     cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaaa                                             >  1:1569965/1‑48 (MQ=255)
                                       aCATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                    <  1:695384/71‑1 (MQ=255)
                                                 cAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTa           <  1:1740823/70‑1 (MQ=255)
                                                         cAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGt  <  1:1337937/71‑1 (MQ=255)
                                                           aaaaaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGa                <  1:1404281/55‑1 (MQ=255)
                                                           aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGtt                               <  1:979056/40‑1 (MQ=255)
                                                           aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAg        <  1:1724951/63‑1 (MQ=255)
                                                           aaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCggg    <  1:1213170/67‑1 (MQ=255)
                                                              aaaGTAACGTACATACCGATAAACAGCCCCGGATGatt                              >  1:1614967/1‑38 (MQ=38)
                                                              aaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCg      >  1:1352680/1‑62 (MQ=255)
                                                                  |                                                             
CCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGT  >  W3110S.gb/478277‑478404

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: