Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,559,619 A→G 16.7% 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
*AP0090481,559,6190AG16.7% 127.3 / 19.3 54L45S (TTA→TCA) ddpFD‑Ala‑D‑Ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base G (6/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACC  >  AP009048/1559551‑1559686
                                                                    |                                                                   
gTTGCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAg                                                                   <  1:3202542/71‑1 (MQ=255)
  tGCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                 >  1:1121593/1‑71 (MQ=255)
   gcagcaTACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGt                                                                 <  1:1150500/70‑1 (MQ=255)
    cagcaTACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGttt                                                               <  1:1155663/71‑1 (MQ=255)
           cccATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTcaccac                                                         <  1:4118062/70‑1 (MQ=255)
              aTTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCAcggcg                                                     >  1:4489640/1‑71 (MQ=255)
                   aaGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCt                                                <  1:2865391/71‑1 (MQ=255)
                      cTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGt                                              >  1:1190862/1‑70 (MQ=255)
                      cTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTa                                             <  1:4185922/71‑1 (MQ=255)
                         tgCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGtt                                                                <  1:2814824/49‑1 (MQ=255)
                         tgCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAAt                                          >  1:336067/1‑71 (MQ=255)
                          gCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGt                                              >  1:251648/1‑66 (MQ=255)
                           cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTa                                                                     <  1:1734476/42‑1 (MQ=255)
                           cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAg                                                                   >  1:3530624/1‑44 (MQ=255)
                           cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAg                                                                   >  1:3525427/1‑44 (MQ=255)
                           cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                         >  1:3708881/1‑70 (MQ=255)
                           cGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                         <  1:2525942/70‑1 (MQ=255)
                                  gCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCAt                                 >  1:3646649/1‑71 (MQ=255)
                                   cTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCAt                                 >  1:818632/1‑70 (MQ=255)
                                    ttttGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCAtt                                >  1:3609432/1‑70 (MQ=255)
                                    ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCAtt                                >  1:1290786/1‑70 (MQ=255)
                                    ttttGCCGCAGCCTGACTCGCCGACGATCCCTAAGATTTCACCACGGCGGATCTGTAAAt                                          <  1:3639784/60‑1 (MQ=255)
                                     tttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCAcggcg                                                     >  1:2138230/1‑48 (MQ=255)
                                     tttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTaa                                            <  1:3131424/57‑1 (MQ=255)
                                      ttGCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTaaa                                           >  1:2261784/1‑57 (MQ=255)
                                      ttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTcac                                                            >  1:4253384/1‑40 (MQ=255)
                                      ttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTaa                              >  1:3246250/1‑70 (MQ=255)
                                      ttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTa                               <  1:3670494/69‑1 (MQ=255)
                                      ttGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAAt                             >  1:3251674/1‑71 (MQ=255)
                                        gCCGCAGCCTGACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATgg                           >  1:3759855/1‑71 (MQ=255)
                                           gCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCAt                        >  1:2455597/1‑71 (MQ=255)
                                           gCAGCCTGACTCGCCGACGATCCATAAGGTTTCACCACCGCGGATCTGTAAATc                                         <  1:4496497/54‑1 (MQ=255)
                                            cAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATg                       <  1:2121503/71‑1 (MQ=255)
                                             aGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTa                                             <  1:583453/48‑1 (MQ=255)
                                             aGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTa                                             <  1:3714317/48‑1 (MQ=255)
                                               ccTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATAc                                    <  1:4129650/55‑1 (MQ=255)
                                                cTGACTCGCCGACGATCCCTAAGGTTTCACCAcggc                                                      >  1:851966/1‑36 (MQ=255)
                                                  gACTCGCCGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATg                 >  1:894097/1‑71 (MQ=255)
                                                   aCTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATg                 <  1:1864445/70‑1 (MQ=255)
                                                   aCTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGt                >  1:950113/1‑71 (MQ=255)
                                                    cTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCAtt                                <  1:4612328/54‑1 (MQ=255)
                                                     tCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATc                                         <  1:1099224/44‑1 (MQ=255)
                                                        ccGACGATCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCg            <  1:4592418/70‑1 (MQ=255)
                                                         cgacgaTCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTa          <  1:3175729/71‑1 (MQ=255)
                                                          gacgaTCCCTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAAt                             >  1:2573520/1‑51 (MQ=255)
                                                          gacgaTCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAg         <  1:3852021/71‑1 (MQ=255)
                                                            cgaTCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGtt       <  1:2876711/71‑1 (MQ=255)
                                                              aTCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGaa                     <  1:3199773/55‑1 (MQ=255)
                                                               tCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACAt                  >  1:4001565/1‑57 (MQ=255)
                                                               tCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACAt                  >  1:507630/1‑57 (MQ=255)
                                                               tCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTa    >  1:4371175/1‑71 (MQ=255)
                                                                 ccTGAGGTTTCACCACGGCGGATCTGTAAATCAATACCa                                  >  1:2888600/1‑39 (MQ=255)
                                                                  cTAAGGTTTCCCCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTAcc  >  1:3558232/1‑70 (MQ=255)
                                                                  cTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCg            <  1:1014528/60‑1 (MQ=255)
                                                                  cTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTAcc  >  1:3865362/1‑70 (MQ=255)
                                                                    |                                                                   
GTTGCAGCATACCCATTAAAAGCTGTGCGAGGGTGCTTTTGCCGCAGCCTGACTCGCCGACGATCCCTAAGGTTTCACCACGGCGGATCTGTAAATCAATACCATTAATGGCATGAACATGTTCCGTAGTTTTACC  >  AP009048/1559551‑1559686

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: