Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 920,084 C→T 9.3% W172* (TGG→TAG)  mpaA ← murein peptide amidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE920,0840CT 9.3% 186.8 / 7.1 64W172* (TGG→TAG) mpaAmurein peptide amidase A
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base T (2/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACC  >  minE/920018‑920154
                                                                  |                                                                      
caTTCGCCTAATTCGCTGTGTCTGGTATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggc                                                                    <  1:4557444/71‑1 (MQ=255)
caTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcagg                                                                     >  1:520589/1‑70 (MQ=255)
 aTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggca                                                                   <  1:1572174/71‑1 (MQ=255)
 aTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcagg                                                                     >  1:1281290/1‑69 (MQ=255)
  ttCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggcag                                                                  <  1:3494691/71‑1 (MQ=255)
    cGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggcagg                                                                 <  1:1703677/70‑1 (MQ=255)
     gCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGt                                                                >  1:467916/1‑70 (MQ=255)
      ccTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTg                                                              <  1:88765/71‑1 (MQ=255)
      ccTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTg                                                              <  1:887021/71‑1 (MQ=255)
       cTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggcag                                                                  <  1:1349134/66‑1 (MQ=255)
          aTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATg                                                          >  1:3843588/1‑71 (MQ=255)
          aTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATg                                                          >  1:1403013/1‑71 (MQ=255)
           ttCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggca                                                                   <  1:3059580/61‑1 (MQ=255)
           ttCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGc                                                         >  1:1961311/1‑71 (MQ=255)
            tCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCg                                                        >  1:1003633/1‑71 (MQ=255)
                    tCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGtata                                                 <  1:1099733/70‑1 (MQ=255)
                      tGGGATCTTCAATACAGGCCAGTTGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGc                                                         >  1:368652/1‑60 (MQ=255)
                       gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTg                                                              <  1:146915/54‑1 (MQ=255)
                       gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCt                                             >  1:2679324/1‑71 (MQ=255)
                        ggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTg                                            >  1:3638861/1‑71 (MQ=255)
                         gATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGTATAAGCTgg                                           <  1:2703481/71‑1 (MQ=255)
                          aTCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTgg                                           >  1:758982/1‑70 (MQ=255)
                           tCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGaa                                                            >  1:3357586/1‑52 (MQ=255)
                           tCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGtata                                                 <  1:4438473/63‑1 (MQ=255)
                             ttCAATACAGGCCAGTGGATCGTGGAATGAGACCACCca                                                                       >  1:3808817/1‑39 (MQ=255)
                             ttCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGaa                                                            >  1:994451/1‑50 (MQ=255)
                               cAATACAGGCCAGTGGATCGTGGAGTGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGca                                         <  1:4615977/67‑1 (MQ=255)
                               cAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGta                                                   >  1:2784143/1‑57 (MQ=255)
                               cAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGta                                                   >  1:1221011/1‑57 (MQ=255)
                               cAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGca                                         <  1:1297782/67‑1 (MQ=255)
                               cAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGc                                          <  1:1339821/66‑1 (MQ=255)
                               cAATACAGGCCAGCGGATCGTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGta                                                   >  1:3935631/1‑57 (MQ=255)
                               cAATACAGGCCAGCGGATCGTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGta                                                   >  1:1721390/1‑57 (MQ=255)
                                aaTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAgcg                                    <  1:1649287/71‑1 (MQ=255)
                                 aTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAgcgc                                   <  1:169856/71‑1 (MQ=255)
                                 aTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAgcgc                                   <  1:3593169/71‑1 (MQ=255)
                                  tACAGGCCAGTGGATCGTGGAATGAGACCACCcaggcagg                                                                 <  1:900735/40‑1 (MQ=255)
                                  tACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGcc                                  >  1:1620520/1‑71 (MQ=255)
                                  tACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGcc                                  <  1:2017667/71‑1 (MQ=255)
                                    cAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTg                                                              <  1:2485701/41‑1 (MQ=255)
                                    cAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCTGTGTATAAGCTGGCACAGCGCCTg                                >  1:4194651/1‑71 (MQ=255)
                                    cAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGgt                                                      >  1:2595916/1‑49 (MQ=255)
                                    cAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCt                                 <  1:3928920/70‑1 (MQ=255)
                                       gCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGGGTATAAGCTGGCACAGCGCCTGGGt                             <  1:906530/71‑1 (MQ=255)
                                          aGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTc                          <  1:2337777/71‑1 (MQ=255)
                                           gTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGc                                                         >  1:720237/1‑39 (MQ=255)
                                           gTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTc                          <  1:3099066/70‑1 (MQ=255)
                                             ggATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAg                                      >  1:3977083/1‑56 (MQ=255)
                                             ggATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAgg                       >  1:3570945/1‑71 (MQ=255)
                                              gATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAgcg                                    >  1:4565443/1‑57 (MQ=255)
                                              gATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGt                      >  1:2255275/1‑71 (MQ=255)
                                              gATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGt                      >  1:5638/1‑71 (MQ=255)
                                                tCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGgt                                                      >  1:3921121/1‑37 (MQ=255)
                                                tCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTggg                              >  1:1563783/1‑61 (MQ=255)
                                                tCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTggg                              >  1:1952731/1‑61 (MQ=255)
                                                tCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTc                    <  1:3795953/71‑1 (MQ=255)
                                                 cGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAg                                               >  1:3102913/1‑43 (MQ=255)
                                                  gTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGTATAAGTTgg                                           <  1:1620090/46‑1 (MQ=255)
                                                  gTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGTATAAGTTgg                                           <  1:3919904/46‑1 (MQ=255)
                                                  gTGGAATGAGACCACCTAGGCAGGTTGAATGCGGTGTATAAGTTgg                                           <  1:1673820/46‑1 (MQ=255)
                                                     gAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGc               <  1:3302712/71‑1 (MQ=255)
                                                      aaTGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGcc              <  1:3335136/71‑1 (MQ=255)
                                                       aTGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCa             <  1:1674359/71‑1 (MQ=255)
                                                        cgagaCCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAgg                       <  1:3372881/59‑1 (MQ=255)
                                                          agaCCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGc                                          >  1:2589339/1‑39 (MQ=255)
                                                          agaCCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGt          >  1:213297/1‑71 (MQ=255)
                                                          agaCCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGt          <  1:1268309/71‑1 (MQ=255)
                                                            accaccCATGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGttt        <  1:2112003/71‑1 (MQ=255)
                                                            accaccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGttt        <  1:1432632/71‑1 (MQ=255)
                                                            accaccCAGGAAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGttt        >  1:3948716/1‑71 (MQ=255)
                                                             ccaccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGttt        <  1:1121751/70‑1 (MQ=255)
                                                              caccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGt      <  1:4107779/71‑1 (MQ=255)
                                                              caccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGt      >  1:830108/1‑71 (MQ=255)
                                                              caccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGt      >  1:2478356/1‑71 (MQ=255)
                                                                cccAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCa    >  1:2049317/1‑71 (MQ=255)
                                                                cccAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCa    <  1:3900450/71‑1 (MQ=255)
                                                                 ccAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc   <  1:1274796/71‑1 (MQ=255)
                                                                 ccAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc   <  1:4540971/71‑1 (MQ=255)
                                                                  caggcaggTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAcc  >  1:2234420/1‑71 (MQ=255)
                                                                  |                                                                      
CATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACC  >  minE/920018‑920154

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