Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,783,789 C→A 28.5% D72Y (GAC→TAC)  yidI ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,783,7890CA28.5% 44.7 / 8.1 21D72Y (GAC→TAC) yidIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (5/10);  new base A (3/3);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTAGCGCAGCAATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTCCGGCGAGAAGACCTCTTTTAATT  >  W3110S.gb/3783726‑3783853
                                                               |                                                                
aaaGTAGCGCAGCAATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTACTCtttt                                                           >  1:1150126/1‑71 (MQ=255)
 aaGTAGCGCAGCAATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCtttt                                                           <  1:1045318/70‑1 (MQ=255)
        gcagcaATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCAc                                                   <  1:837753/71‑1 (MQ=255)
          agcaATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACgg                                                 <  1:815320/71‑1 (MQ=255)
              aTGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCtttttt                                                         >  1:1325287/1‑59 (MQ=255)
              aTGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCtttttt                                                         >  1:177343/1‑59 (MQ=255)
               tGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCtttttt                                                         <  1:1256160/58‑1 (MQ=255)
                     gTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTaa                                           <  1:1643021/66‑1 (MQ=255)
                     gTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTaa                                           <  1:521913/66‑1 (MQ=255)
                     gTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTACTCTTTTTTCTCCACGGTGGTAAtttt                                       <  1:1089239/70‑1 (MQ=255)
                     gTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTACTCTTTTTTCTCCACGGTGGTAAtttt                                       <  1:379345/70‑1 (MQ=255)
                            aCCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCttt                                                            <  1:1004118/42‑1 (MQ=255)
                               gCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTa                            <  1:268374/71‑1 (MQ=255)
                                   tgattgaGGATTTTATCGACATCTACGTACTCtttt                                                           <  1:13997/36‑1 (MQ=255)
                                    gattgaGGATTTTATCGACATCTACGTACTCTTTTTTCTCCACgg                                                 >  1:294131/1‑45 (MQ=255)
                                            aTTTTATCGACAGCTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTCCGGCGagaag               <  1:1765539/71‑1 (MQ=255)
                                              tttATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCtt                               <  1:811770/53‑1 (MQ=255)
                                                tATCGACATCTACGTACTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTCCGGCGAGAAGACct           >  1:55516/1‑71 (MQ=255)
                                                         cTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTa                            >  1:649601/1‑45 (MQ=255)
                                                          tACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTcc                        >  1:769530/1‑48 (MQ=255)
                                                          tACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTCCGGCGAGAAGACCTCTTTTAAtt  >  1:1194531/1‑70 (MQ=255)
                                                               |                                                                
AAAGTAGCGCAGCAATGGCCAGTGCGATACCGCTTTGATTGAGGATTTTATCGACATCTACGTCCTCTTTTTTCTCCACGGTGGTAATTTTTTCGCCTTTTATTCCGGCGAGAAGACCTCTTTTAATT  >  W3110S.gb/3783726‑3783853

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: