Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,559,619 A→G 28.6% 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,6190AG28.6% 48.6 / 25.3 35L45S (TTA→TCA) ddpFD‑Ala‑D‑Ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (13/12);  new base G (4/6);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACCAAG  >  W3110S.gb/1559554‑1559689
                                                                 |                                                                      
gcagcaTACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGtt                                                                   >  1:1322850/1‑71 (MQ=255)
   gcaTACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTca                                                                >  1:900134/1‑71 (MQ=255)
           aTTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCAcggcg                                                        <  1:1331788/71‑1 (MQ=255)
            ttAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCAcggcg                                                        >  1:173293/1‑70 (MQ=255)
             tAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGa                                                      >  1:3730603/1‑71 (MQ=255)
             tAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGa                                                      >  1:4603117/1‑71 (MQ=255)
              aaaaGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGa                                                      >  1:4161123/1‑70 (MQ=255)
              aaaaGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGAt                                                     >  1:946336/1‑71 (MQ=255)
               aaaGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTcac                                                               <  1:3237602/60‑1 (MQ=255)
                    tgtgCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTcac                                                               <  1:226711/55‑1 (MQ=255)
                      tgCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAAt                                             >  1:4307612/1‑71 (MQ=255)
                        cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                            >  1:4602496/1‑70 (MQ=255)
                        cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                            <  1:2417855/70‑1 (MQ=255)
                        cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                            <  1:1058167/70‑1 (MQ=255)
                              tGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATc                                            <  1:708815/64‑1 (MQ=255)
                              tGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATc                                            <  1:1997883/64‑1 (MQ=255)
                              tGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCa                                     <  1:4401949/71‑1 (MQ=255)
                                 ttttTCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                    <  1:4758631/37‑1 (MQ=255)
                                 ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                    <  1:845344/37‑1 (MQ=255)
                                 ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                    <  1:3993366/37‑1 (MQ=255)
                                 ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                    <  1:2993828/37‑1 (MQ=255)
                                 ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTcac                                                               >  1:2803187/1‑42 (MQ=255)
                                          aGCCTTACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATg                          >  1:1424545/1‑70 (MQ=255)
                                          aGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGa                         <  1:2591952/71‑1 (MQ=255)
                                          aGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGa                         <  1:2161548/71‑1 (MQ=255)
                                           gCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGaa                        <  1:3001824/71‑1 (MQ=255)
                                           gCCTGACTCGCCGACGATCCATGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGaa                        <  1:3490916/71‑1 (MQ=255)
                                                 cTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGt                                                 >  1:4237356/1‑40 (MQ=255)
                                                 cTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGt                                                 >  1:51809/1‑40 (MQ=255)
                                                    gCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCg               <  1:921013/71‑1 (MQ=255)
                                                    gCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCg               >  1:4131194/1‑71 (MQ=255)
                                                    gCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCg               >  1:3232587/1‑71 (MQ=255)
                                                       gacgaTCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAAc                       <  1:1821480/60‑1 (MQ=255)
                                                             cccTAAGGTTTCACCACTGCGGATCTGTAAATCaa                                          >  1:2397280/1‑35 (MQ=255)
                                                             cccTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGt           >  1:2122027/1‑66 (MQ=255)
                                                                 aaGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACCaa   >  1:3943476/1‑70 (MQ=255)
                                                                 aaGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACCAAg  >  1:632189/1‑71 (MQ=255)
                                                                 aaGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACCAAg  >  1:4307298/1‑71 (MQ=255)
                                                                 |                                                                      
GCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACCAAG  >  W3110S.gb/1559554‑1559689

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