Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 197,387 G→C 5.2% S281T (AGT→ACT)  yaeL → zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048197,3870GC 5.2% 233.2 / 2.4 76S281T (AGT→ACT) yaeLzinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base C (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATT  >  AP009048/197321‑197457
                                                                  |                                                                      
ccTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    <  1:2658221/71‑1 (MQ=255)
 cTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCt                                                                   >  1:5470426/1‑71 (MQ=255)
   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:7466262/1‑71 (MQ=255)
   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:7457029/1‑71 (MQ=255)
   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:6567932/1‑71 (MQ=255)
   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:4188778/1‑71 (MQ=255)
    tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 <  1:945900/70‑1 (MQ=255)
    tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                >  1:3317047/1‑71 (MQ=255)
      tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                              <  1:3834946/71‑1 (MQ=255)
      tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                              <  1:1385580/71‑1 (MQ=255)
        aTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                            <  1:4209517/71‑1 (MQ=255)
         tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                >  1:1515477/1‑66 (MQ=255)
         tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                            >  1:3452303/1‑70 (MQ=255)
          gCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:5793100/1‑71 (MQ=255)
            ttGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        >  1:7486774/1‑71 (MQ=255)
              gTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        <  1:4798441/69‑1 (MQ=255)
               tCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    <  1:2919299/56‑1 (MQ=255)
               tCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCtt                                                             <  1:2054973/63‑1 (MQ=255)
                 cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                   >  1:7789048/1‑71 (MQ=255)
                 cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                   <  1:3041545/71‑1 (MQ=255)
                    gATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATccc                                                >  1:3979165/1‑71 (MQ=255)
                     aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCg                                               <  1:3207197/71‑1 (MQ=255)
                     aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGGCCCTTGTCTTTGACATTAATCCCg                                               <  1:3427132/71‑1 (MQ=255)
                         cccGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             >  1:2190408/1‑69 (MQ=255)
                            gggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTT‑TCTTTGACATTAATCCCGGAGAGt                                         >  1:1452639/1‑69 (MQ=255)
                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             <  1:4899522/65‑1 (MQ=255)
                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:5833893/71‑1 (MQ=255)
                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:3343685/71‑1 (MQ=255)
                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                        >  1:2436857/1‑70 (MQ=255)
                              gTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaaa                                      <  1:1821578/71‑1 (MQ=255)
                               tAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAAc                                     >  1:3602359/1‑71 (MQ=255)
                                aaaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    >  1:7649139/1‑39 (MQ=255)
                                aaaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             <  1:4636894/62‑1 (MQ=255)
                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:1416368/1‑48 (MQ=255)
                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                         <  1:532508/49‑1 (MQ=255)
                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                         <  1:523108/49‑1 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTg                                                           >  1:1695216/1‑44 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTaa                                                    >  1:204824/1‑51 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:3307857/1‑70 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:897626/1‑70 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:3053792/1‑70 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:1993208/1‑70 (MQ=255)
                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:237401/1‑70 (MQ=255)
                                     cTTAGCGTTAGAAATCGAAAGGCAGTGGAGTCCCTTGTCTTTGACa                                                        >  1:7258646/1‑46 (MQ=255)
                                     cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        >  1:4291081/1‑46 (MQ=255)
                                     cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTa                               <  1:1676505/71‑1 (MQ=255)
                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                <  1:5600713/36‑1 (MQ=255)
                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                <  1:3673409/36‑1 (MQ=255)
                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                        >  1:5032656/1‑60 (MQ=255)
                                        aGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             >  1:5177872/1‑54 (MQ=255)
                                         gCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:678515/1‑40 (MQ=255)
                                           gTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                         <  1:7501580/71‑1 (MQ=255)
                                            ttAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                         >  1:6033077/1‑70 (MQ=255)
                                                aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGc                     >  1:1378508/1‑70 (MQ=255)
                                                aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCg                    <  1:3119412/71‑1 (MQ=255)
                                                  aTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATgg                           <  1:3044962/62‑1 (MQ=255)
                                                     gAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCTATTgg               <  1:1684330/71‑1 (MQ=255)
                                                      aaaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt              <  1:5635083/71‑1 (MQ=255)
                                                      aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                >  1:7653797/1‑69 (MQ=25)
                                                      aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                >  1:4534769/1‑69 (MQ=25)
                                                       aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTTAAGCGATTGGtt             >  1:6626766/1‑71 (MQ=255)
                                                       aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGtt             <  1:3489967/71‑1 (MQ=255)
                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:6244903/44‑1 (MQ=255)
                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:1576187/44‑1 (MQ=255)
                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt            <  1:2736806/71‑1 (MQ=255)
                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt            >  1:4424518/1‑71 (MQ=255)
                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt              <  1:3372853/69‑1 (MQ=255)
                                                          gCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTg          <  1:3276904/71‑1 (MQ=255)
                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGt                          >  1:7321514/1‑54 (MQ=255)
                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:5553728/71‑1 (MQ=255)
                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:4922697/71‑1 (MQ=255)
                                                            aGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:6787180/70‑1 (MQ=255)
                                                              gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg      <  1:285673/71‑1 (MQ=255)
                                                              gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg      >  1:4511659/1‑71 (MQ=255)
                                                              gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                        <  1:6114528/53‑1 (MQ=38)
                                                              gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                        <  1:1824813/53‑1 (MQ=38)
                                                                  gTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTAtt  <  1:890886/71‑1 (MQ=255)
                                                                  |                                                                      
CCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATT  >  AP009048/197321‑197457

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