Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,874,329 A→G 14.1% I191T (ATT→ACT)  ygcF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,874,3290AG14.1% 217.2 / 28.9 92I191T (ATT→ACT) ygcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/34);  new base G (7/6);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTCA  >  minE/1874260‑1874397
                                                                     |                                                                    
aTCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAAt                                                                     >  1:323758/1‑71 (MQ=255)
 tCGACAAACGCGAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATc                                                                    <  1:1608041/71‑1 (MQ=255)
  cGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCg                                                                   >  1:2047372/1‑71 (MQ=255)
   gACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCgg                                                                  >  1:769002/1‑71 (MQ=255)
     cAAACGCCAATTACGCGCAATGCAGGTTTAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               <  1:881480/71‑1 (MQ=255)
      aaaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               >  1:1580084/1‑71 (MQ=255)
      aaaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               <  1:988088/71‑1 (MQ=255)
      aaaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               >  1:1812777/1‑71 (MQ=255)
      aaaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               >  1:548926/1‑71 (MQ=255)
       aaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCt                                                                >  1:640291/1‑69 (MQ=255)
       aaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               >  1:1945813/1‑70 (MQ=255)
       aaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGc                                                              <  1:1172703/71‑1 (MQ=255)
       aaCGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGc                                                              <  1:2203789/71‑1 (MQ=255)
           ccAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTg                                                          <  1:2070837/71‑1 (MQ=255)
            cAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGc                                                         <  1:1136775/71‑1 (MQ=255)
               ttACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAACCGGCTGCAGTGCAAt                                                      <  1:1443476/71‑1 (MQ=255)
                 aCGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAAt                                                                     >  1:1034232/1‑54 (MQ=255)
                 aCGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATg                                                     <  1:1022664/70‑1 (MQ=255)
                 aCGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGATGCAGTGCAAt                                                      <  1:363815/69‑1 (MQ=255)
                  cgcgCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCaa                                                       <  1:848509/67‑1 (MQ=255)
                     gcAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTaa                                                                      >  1:1264026/1‑49 (MQ=255)
                     gcAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTcg                                                >  1:876365/1‑71 (MQ=255)
                     gcAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTc                                                 >  1:2041377/1‑70 (MQ=255)
                     gcAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTc                                                 >  1:252446/1‑70 (MQ=255)
                      cAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTcgc                                               <  1:780853/71‑1 (MQ=255)
                       aaTGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTcgcg                                              >  1:511683/1‑71 (MQ=255)
                         tGCAGGTTTCAATGCACAAACGTGTGGCTTCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGt                                            >  1:1392467/1‑71 (MQ=255)
                          gCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:883866/1‑71 (MQ=255)
                          gCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:650848/1‑71 (MQ=255)
                          gCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:1726879/1‑71 (MQ=255)
                           cAGGTTTCAATGCACTAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:809218/1‑70 (MQ=255)
                           cAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:886078/1‑70 (MQ=255)
                           cAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGtt                                           >  1:185638/1‑70 (MQ=255)
                              gTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTg                                                               <  1:986935/47‑1 (MQ=255)
                              gTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatc                                       <  1:311595/71‑1 (MQ=255)
                               tttCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAAt                                                      <  1:715661/55‑1 (MQ=255)
                                ttCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGAc                                                   >  1:1179463/1‑57 (MQ=255)
                                 tCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGAc                                                   >  1:2037286/1‑56 (MQ=255)
                                 tCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcatc                                    >  1:955250/1‑71 (MQ=255)
                                 tCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcatc                                    >  1:2111942/1‑71 (MQ=255)
                                 tCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcatc                                    >  1:1769327/1‑71 (MQ=255)
                                 tCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGAAATGACTCGCGGTTTatcatc                                    >  1:1368774/1‑71 (MQ=255)
                                  cAATGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCg                                   >  1:826872/1‑71 (MQ=255)
                                   aaTGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCgg                                  >  1:280515/1‑71 (MQ=255)
                                   aaTGCACAAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCg                                   >  1:1954765/1‑70 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcatc                                    <  1:1325346/67‑1 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcatc                                    <  1:990160/67‑1 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcat                                     >  1:1135183/1‑66 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatcat                                     >  1:378564/1‑66 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTatc                                       <  1:1116464/64‑1 (MQ=255)
                                     tGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGt                                 >  1:532782/1‑70 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTg                                                          >  1:326987/1‑43 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAg                              >  1:276837/1‑71 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTc                                                 >  1:130488/1‑52 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGt                                            >  1:966400/1‑57 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGt                                            >  1:1819316/1‑57 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCg                                   >  1:80131/1‑66 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCg                                   >  1:1679210/1‑66 (MQ=255)
                                       cacaAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAg                              >  1:1481507/1‑71 (MQ=255)
                                          aaaCGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTat                                        >  1:521212/1‑58 (MQ=255)
                                           aaCGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGtt                          >  1:989942/1‑71 (MQ=255)
                                           aaCGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGtt                          >  1:1511809/1‑71 (MQ=255)
                                              gtgtGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGcc                       <  1:248147/71‑1 (MQ=255)
                                              gtgtGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGcc                       <  1:1314592/71‑1 (MQ=255)
                                              gtgtGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGcc                       <  1:602157/71‑1 (MQ=255)
                                                gtgGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCag                     >  1:592173/1‑71 (MQ=255)
                                                 tgGCATCATCCTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCagt                    <  1:1061552/71‑1 (MQ=255)
                                                 tgGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCagt                    <  1:1584708/71‑1 (MQ=255)
                                                 tgGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCag                     <  1:1223095/70‑1 (MQ=255)
                                                 tgGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCag                     >  1:1547261/1‑70 (MQ=255)
                                                 tgGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCag                     <  1:356028/70‑1 (MQ=255)
                                                  ggCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCag                     >  1:128984/1‑69 (MQ=255)
                                                      tcatcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCagt                    <  1:1316194/66‑1 (MQ=255)
                                                        atcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCa              >  1:314278/1‑70 (MQ=255)
                                                        atcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCAt             >  1:1624294/1‑71 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            <  1:797987/71‑1 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            <  1:228984/71‑1 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            <  1:614528/71‑1 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            <  1:504870/71‑1 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTCCCAGTAGTTCATc            <  1:47546/71‑1 (MQ=255)
                                                         tcCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGGCAGCGTTGCCAGTAGTTCATc            <  1:853893/71‑1 (MQ=255)
                                                          ccTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTat                                        <  1:1845491/42‑1 (MQ=255)
                                                          ccTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            >  1:1964057/1‑70 (MQ=255)
                                                          ccTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATc            >  1:1562860/1‑70 (MQ=255)
                                                           cTTTTGGCTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATcc           <  1:124716/70‑1 (MQ=255)
                                                           cTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATcc           <  1:502113/70‑1 (MQ=255)
                                                           cTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATcc           <  1:739417/70‑1 (MQ=255)
                                                            ttttGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGc                             <  1:550623/51‑1 (MQ=255)
                                                              ttGGCTAATCGGATGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCagtag                  >  1:326755/1‑60 (MQ=255)
                                                                  cTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTc   <  1:2007573/71‑1 (MQ=255)
                                                                  cTAGTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTc   <  1:1519919/71‑1 (MQ=255)
                                                                  cTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTc   <  1:992114/71‑1 (MQ=255)
                                                                   tAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTc   >  1:1453954/1‑70 (MQ=255)
                                                                    aaTCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTCa  >  1:1578003/1‑70 (MQ=255)
                                                                     |                                                                    
ATCGACAAACGCCAATTACGCGCAATGCAGGTTTCAATGCACAAACGTGTGGCATCATCCTTTTGGCTAATCGGCTGCAGTGCAATGACTCGCGGTTTATCATCGGTCAGCGTTGCCAGTAGTTCATCCAGTGCTTCA  >  minE/1874260‑1874397

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: