Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,028,050 G→A 5.2% I228I (ATC→ATT parC ← DNA topoisomerase IV, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,028,0500GA 5.2% 305.5 / 4.1 96I228I (ATC→ATTparCDNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base G (42/49);  new base A (2/3);  total (44/52)
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.11e-01

CGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCT  >  minE/2027983‑2028120
                                                                   |                                                                      
cgcCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtc                                                                     >  1:3652996/1‑71 (MQ=255)
 gcaaTCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtct                                                                    <  1:3059955/68‑1 (MQ=255)
 gcCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCTTAGATTTTACGGAtct                                                                    <  1:1817566/71‑1 (MQ=255)
 gcCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtct                                                                    <  1:765495/71‑1 (MQ=255)
  ccATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg                                                                        <  1:2630985/66‑1 (MQ=255)
  ccATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtctc                                                                   >  1:360961/1‑71 (MQ=255)
  ccATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtctc                                                                   <  1:1802292/71‑1 (MQ=255)
     tcttctTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg                                                                        <  1:4343179/63‑1 (MQ=255)
     tcttctTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgc                                                                >  1:4133536/1‑71 (MQ=255)
       ttctttctTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGa                                                                       >  1:1088265/1‑62 (MQ=255)
       ttctttctTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAt                                                                      >  1:3840424/1‑63 (MQ=255)
       ttctttctTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgcgc                                                              >  1:809426/1‑71 (MQ=255)
       ttctttctTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgcgc                                                              >  1:248189/1‑71 (MQ=255)
       ttctttctTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgc                                                                >  1:400177/1‑69 (MQ=255)
        tctttcttCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgcgc                                                              >  1:1484122/1‑70 (MQ=255)
            tcttcCACACCGCGCGCATACGCACTGAACCACTTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGa                                                          <  1:2702278/70‑1 (MQ=255)
            tcttcCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg                                                           >  1:4675969/1‑69 (MQ=255)
            tcttcCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg                                                           <  1:4411750/69‑1 (MQ=255)
            tcttcCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg                                                           <  1:2032721/69‑1 (MQ=255)
             cttcCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAg                                                        >  1:1272233/1‑71 (MQ=255)
                ccACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTg                                                      <  1:1797134/70‑1 (MQ=255)
                 cacaCCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAg                                                        <  1:1482304/67‑1 (MQ=255)
                  acacCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGAATCTCGGCGCGCGAAGTGAt                                                    >  1:4089924/1‑70 (MQ=255)
                  acacCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGAATCTCGGCGCGCGAAGTGAt                                                    >  1:4075234/1‑70 (MQ=255)
                   cacCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtctct                                                                  >  1:3015152/1‑54 (MQ=255)
                   cacCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATaa                                                  >  1:756882/1‑71 (MQ=255)
                      cgcgcgCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAt                                                                      <  1:1284228/48‑1 (MQ=255)
                      cgcgcgCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAttt                                               <  1:1911660/71‑1 (MQ=255)
                      cgcgcgCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAtt                                                <  1:3042746/70‑1 (MQ=255)
                      cgcgcgCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAtt                                                >  1:376701/1‑70 (MQ=255)
                        cgcgcATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtct                                                                    >  1:2754368/1‑48 (MQ=255)
                           gcATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGAATCTCGGCGCGCGAAGTGATAATTTCCGCt                                          <  1:4261676/71‑1 (MQ=255)
                             aTACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCtt                                         >  1:2148494/1‑70 (MQ=255)
                             aTACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGAATCTCGGCGCGCGAAGt                                                       <  1:218165/56‑1 (MQ=255)
                              tACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:4343049/71‑1 (MQ=255)
                              tACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:4529776/71‑1 (MQ=255)
                              tACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGAATCTCGGCGCGCGAAGt                                                       <  1:3871986/55‑1 (MQ=255)
                                cGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAg                                      <  1:4189866/70‑1 (MQ=255)
                                cGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAg                                      >  1:3781655/1‑70 (MQ=255)
                                cGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGt                                     >  1:1126606/1‑71 (MQ=255)
                                cGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGt                                     >  1:2437685/1‑71 (MQ=255)
                                  cACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCg                                   >  1:4296029/1‑71 (MQ=255)
                                   aCTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCgg                                  <  1:1353586/71‑1 (MQ=255)
                                   aCTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCgg                                  <  1:636099/71‑1 (MQ=255)
                                   aCTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCg                                   >  1:901336/1‑70 (MQ=255)
                                   aCTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCCCGAAGTGATAATTTCCGCTTCAGTCgg                                  <  1:3642776/71‑1 (MQ=255)
                                    cTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg                                                           <  1:1896858/45‑1 (MQ=255)
                                    cTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg                                                           <  1:1705908/45‑1 (MQ=255)
                                    cTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAtt                                                <  1:2482555/56‑1 (MQ=255)
                                    cTGAACCACGTCCGTTCTCGTAGATTTCACGGAtct                                                                    <  1:3663786/36‑1 (MQ=255)
                                    cTGAACCACGTCCGTTCACGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGa                                 <  1:2017788/71‑1 (MQ=255)
                                     tGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAttt                                               >  1:3903760/1‑56 (MQ=255)
                                     tGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:1781974/64‑1 (MQ=255)
                                     tGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGa                                 >  1:2301736/1‑70 (MQ=255)
                                      gAACCACGTCCGTTCTCGTAGATTTTACGGATCTGGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGAt                                >  1:2758607/1‑70 (MQ=255)
                                      gAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATaa                                                  >  1:3513738/1‑52 (MQ=255)
                                      gAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGAt                                >  1:22978/1‑70 (MQ=255)
                                      gAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGAt                                >  1:3647858/1‑70 (MQ=255)
                                       aaCCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGAt                                <  1:746565/69‑1 (MQ=255)
                                       aaCCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGAt                                >  1:2422578/1‑69 (MQ=255)
                                       aaCCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATaa                              <  1:2075690/71‑1 (MQ=255)
                                         ccACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCTCGAAGTGAt                                                    <  1:3046334/47‑1 (MQ=255)
                                         ccACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTc                                    >  1:3503804/1‑63 (MQ=255)
                                          cACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTct                                   <  1:1623125/63‑2 (MQ=255)
                                           aCGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGaa                                                         >  1:1620070/1‑40 (MQ=255)
                                           aCGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCg                          >  1:2586279/1‑71 (MQ=255)
                                            cGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCg                          >  1:1714799/1‑70 (MQ=255)
                                             gTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGa                                                          <  1:4384547/37‑1 (MQ=255)
                                             gTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTg                                                      <  1:2755522/41‑1 (MQ=255)
                                              tCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGAt                                                    <  1:3703570/42‑1 (MQ=255)
                                              tCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGcc                       <  1:1580229/71‑1 (MQ=255)
                                              tCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGcc                       <  1:1889194/71‑1 (MQ=255)
                                               ccGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTACGCtt                                         >  1:430644/1‑52 (MQ=255)
                                                cGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTc                                              <  1:2355970/46‑1 (MQ=255)
                                                cGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGcccc                     >  1:4673053/1‑71 (MQ=255)
                                                cGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGccc                      >  1:1361984/1‑70 (MQ=255)
                                                 gTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCtt                                         <  1:3342892/50‑1 (MQ=255)
                                                  ttCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCt                                          <  1:1913873/48‑1 (MQ=255)
                                                  ttCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGa                                 <  1:2951894/57‑1 (MQ=255)
                                                  ttCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTg                   <  1:1129412/71‑1 (MQ=255)
                                                  ttCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTg                   <  1:3391220/71‑1 (MQ=255)
                                                   tctcGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAAtt                                                <  1:3247817/41‑1 (MQ=255)
                                                   tctcGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTg                   >  1:918468/1‑70 (MQ=255)
                                                    ctcGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATc                            >  1:2842316/1‑60 (MQ=255)
                                                    ctcGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATc                            >  1:1375577/1‑60 (MQ=255)
                                                       gTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGa              <  1:942587/71‑1 (MQ=255)
                                                       gTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGa              >  1:225297/1‑71 (MQ=255)
                                                        tAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:4126314/45‑1 (MQ=255)
                                                        tAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:83642/45‑1 (MQ=255)
                                                        tAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCa                                       <  1:3065378/45‑1 (MQ=255)
                                                        tAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGcg                       <  1:491108/61‑2 (MQ=255)
                                                          gATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGt                                     >  1:295206/1‑45 (MQ=255)
                                                             tttACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATcc                           <  1:4456133/52‑1 (MQ=255)
                                                             tttACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATc          >  1:2344916/1‑69 (MQ=255)
                                                               tACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCcag       >  1:1292703/1‑70 (MQ=255)
                                                                aCGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCcagc      <  1:4655860/70‑1 (MQ=255)
                                                                aCGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCcagc      <  1:43382/70‑1 (MQ=255)
                                                                aCGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCcagc      <  1:4018429/70‑1 (MQ=255)
                                                                  ggATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTc                                    >  1:3882894/1‑38 (MQ=255)
                                                                   gATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCt  >  1:2437513/1‑71 (MQ=255)
                                                                   |                                                                      
CGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCGAAGTGATAATTTCCGCTTCAGTCGGATAATCCGGCCCCTGCACGATATCCAGCAGCT  >  minE/2027983‑2028120

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: