Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,841 C→T 5.4% N37N (AAC→AAT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8410CT 5.4% 309.3 / 4.9 92N37N (AAC→AATyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (46/41);  new base T (3/2);  total (49/43)
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.41e-01

TCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATC  >  minE/2569771‑2569904
                                                                      |                                                               
tCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                                                 <  1:3143304/71‑1 (MQ=255)
   ttCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCAt                                                              >  1:607700/1‑71 (MQ=255)
    tCTGGTGATTTTTTATGTAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                                                                <  1:556096/68‑1 (MQ=255)
    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                                                 >  1:139874/1‑67 (MQ=255)
     cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcc                                                            >  1:1418892/1‑71 (MQ=255)
     cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                                             <  1:3615630/70‑1 (MQ=255)
     cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                                             <  1:3232209/70‑1 (MQ=255)
     cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                                             <  1:77264/70‑1 (MQ=255)
        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcccc                                                          <  1:37745/70‑1 (MQ=255)
        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc                                                           <  1:4026137/69‑1 (MQ=255)
         tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                        >  1:3290944/1‑71 (MQ=255)
         tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                        <  1:621373/71‑1 (MQ=255)
         tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                        <  1:1401956/71‑1 (MQ=255)
         tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                        >  1:2614097/1‑71 (MQ=255)
          gATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAATCATcccc                                                          >  1:1274487/1‑68 (MQ=255)
          gATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcccc                                                          >  1:2158427/1‑68 (MQ=255)
            ttttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                     <  1:1132263/71‑1 (MQ=255)
             tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                                                 >  1:452812/1‑58 (MQ=255)
             tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                     >  1:3632706/1‑70 (MQ=255)
                  aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                >  1:3986902/1‑70 (MQ=255)
                  aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                >  1:462296/1‑70 (MQ=255)
                  aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                <  1:1134066/70‑1 (MQ=255)
                   tGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                <  1:1363809/69‑1 (MQ=255)
                   tGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCt                                               <  1:2703948/70‑1 (MQ=255)
                   tGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTa                                              <  1:3212536/71‑1 (MQ=255)
                     gAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAt                                             <  1:262934/70‑1 (MQ=255)
                     gAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAt                                             <  1:4057028/70‑1 (MQ=255)
                       atcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa                                                         >  1:2761807/1‑56 (MQ=255)
                       atcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa                                                         >  1:1354956/1‑56 (MQ=255)
                        tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCTCAACCATCCCCATGCCGACGCTAtt                                            <  1:1205087/68‑1 (MQ=255)
                        tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                                                                <  1:848391/48‑1 (MQ=255)
                        tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          >  1:3473913/1‑70 (MQ=255)
                        tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCAt                                         <  1:3369987/71‑1 (MQ=255)
                         catcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc                                                           >  1:1370937/1‑52 (MQ=255)
                          atcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                                                 <  1:1519897/45‑1 (MQ=255)
                          atcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTa                                              >  1:33530/1‑64 (MQ=255)
                           tcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcg                                      <  1:3889651/71‑1 (MQ=255)
                            cgcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcgg                                     >  1:3130024/1‑71 (MQ=255)
                             gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                                    >  1:2063193/1‑71 (MQ=255)
                             gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                                    >  1:2343356/1‑71 (MQ=255)
                              cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                   <  1:1529839/71‑1 (MQ=255)
                               gATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                            <  1:2515772/61‑1 (MQ=255)
                                 tACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                     <  1:1447604/50‑1 (MQ=255)
                                   cATGAGATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                              >  1:2952413/1‑71 (MQ=255)
                                   cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCt                               <  1:3340647/70‑1 (MQ=255)
                                   cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                              >  1:3311068/1‑71 (MQ=255)
                                   cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                              >  1:1367069/1‑71 (MQ=255)
                                     tGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                                                                <  1:2246928/35‑1 (MQ=255)
                                     tGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                                                                <  1:3143852/35‑1 (MQ=255)
                                     tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                            >  1:356224/1‑71 (MQ=255)
                                     tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                            <  1:2255860/71‑1 (MQ=255)
                                      gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                                             >  1:2638771/1‑37 (MQ=255)
                                      gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                           >  1:2772236/1‑71 (MQ=255)
                                       aCATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                           >  1:4112618/1‑70 (MQ=255)
                                          tCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          <  1:3492799/52‑1 (MQ=255)
                                          tCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                           <  1:193167/67‑1 (MQ=255)
                                           cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          >  1:1407391/1‑51 (MQ=255)
                                           cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          >  1:4085203/1‑51 (MQ=255)
                                            ccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTa                                              >  1:2167772/1‑46 (MQ=255)
                                             cTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTg                    <  1:2881446/71‑1 (MQ=255)
                                              ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                     <  1:2598525/37‑1 (MQ=255)
                                              ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                   <  1:3689817/39‑1 (MQ=255)
                                              ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                            <  1:1027894/46‑1 (MQ=255)
                                              ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGt                   >  1:197004/1‑71 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                   >  1:2969052/1‑37 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                   >  1:2277523/1‑37 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                   >  1:930878/1‑37 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCt                                               >  1:3012275/1‑41 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTc                                           >  1:1891499/1‑45 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:2480107/1‑48 (MQ=255)
                                                aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGtt                  >  1:990417/1‑70 (MQ=255)
                                                 tCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:3233766/1‑47 (MQ=255)
                                                 tCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGttt                 <  1:1208724/70‑1 (MQ=255)
                                                  cTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:1038683/1‑46 (MQ=255)
                                                  cTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:1762898/1‑46 (MQ=255)
                                                  cTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:2361629/1‑46 (MQ=255)
                                                  cTGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCAtgt                                       <  1:725454/47‑3 (MQ=255)
                                                  cTGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc               <  1:1610422/71‑1 (MQ=255)
                                                   tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc               >  1:1228502/1‑70 (MQ=255)
                                                   tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc               >  1:1917940/1‑70 (MQ=255)
                                                     aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                >  1:455242/1‑35 (MQ=255)
                                                     aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                            <  1:2231458/39‑1 (MQ=255)
                                                     aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          <  1:2023526/41‑1 (MQ=255)
                                                     aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTa                <  1:3277672/67‑1 (MQ=255)
                                                      ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          <  1:2592923/40‑1 (MQ=255)
                                                      ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                          >  1:1277890/1‑40 (MQ=255)
                                                      ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg           >  1:1638252/1‑71 (MQ=255)
                                                      ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg           <  1:514100/71‑1 (MQ=255)
                                                      ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg           <  1:1446602/71‑1 (MQ=255)
                                                       tGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAc                                       >  1:2306928/1‑42 (MQ=255)
                                                        gCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                                        >  1:201459/1‑40 (MQ=255)
                                                        gCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                   >  1:2166476/1‑45 (MQ=255)
                                                        gCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa         >  1:2509197/1‑71 (MQ=255)
                                                         ccAAAAAGCGCAACCATCCCCATGCCGACTCTATTCa                                          <  1:1533640/37‑1 (MQ=255)
                                                          cAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCAt        <  1:2097289/70‑1 (MQ=255)
                                                            aaaaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATg       <  1:2705871/69‑1 (MQ=255)
                                                            aaaaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGtt     <  1:1744519/71‑1 (MQ=255)
                                                              aaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTa    >  1:2790158/1‑70 (MQ=255)
                                                               cgcgcAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATc  <  1:2449483/70‑1 (MQ=255)
                                                                  gcAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGtt                  <  1:521538/52‑1 (MQ=255)
                                                                      |                                                               
TCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATC  >  minE/2569771‑2569904

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: