Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,559,619 A→G 50.0% L45S (TTA→TCA)  ddpF ← D‑Ala‑D‑Ala transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,559,6190AG50.0% ‑5.5 / 16.9 14L45S (TTA→TCA) ddpFD‑Ala‑D‑Ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (4/3);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTAC  >  W3110S.gb/1559556‑1559685
                                                               |                                                                  
agcaTACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTc                                                             <  1:1607381/71‑1 (MQ=255)
         aTTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCAcggcg                                                    >  1:1709639/1‑71 (MQ=255)
          ttAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGttt                                                              <  1:695552/60‑1 (MQ=255)
                            tGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCa                                 >  1:1260160/1‑71 (MQ=255)
                                 ttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTa                              <  1:651781/69‑1 (MQ=255)
                                  tGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTcac                                                           >  1:1100024/1‑39 (MQ=255)
                                           cTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTaaa                                          >  1:731339/1‑47 (MQ=255)
                                              aCTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAAt                            <  1:1198121/58‑1 (MQ=255)
                                              aCTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAAt                            <  1:1543711/58‑1 (MQ=255)
                                                tCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTaa                             <  1:132978/55‑1 (MQ=255)
                                                 cGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCaa                                      >  1:318641/1‑45 (MQ=255)
                                                            ccTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAAc                   <  1:1259451/53‑1 (MQ=255)
                                                            ccTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTAc  >  1:1162945/1‑70 (MQ=255)
                                                              tAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTaa                             >  1:1111257/1‑41 (MQ=255)
                                                               |                                                                  
AGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTAC  >  W3110S.gb/1559556‑1559685

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: ‑