Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,859 T→C 8.6% A43A (GCT→GCC yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8590TC 8.6% 262.5 / 13.8 93A43A (GCT→GCCyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base C (3/5);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGC  >  minE/2569790‑2569927
                                                                     |                                                                    
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)
     tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCTCAACCATCCCCATGCCGACGCTAtt                                                                   <  1:1205087/68‑1 (MQ=255)
     tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                                 >  1:3473913/1‑70 (MQ=255)
     tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCAt                                                                <  1:3369987/71‑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:2343356/1‑71 (MQ=255)
          gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                                                           >  1:2063193/1‑71 (MQ=255)
           cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                                          <  1:1529839/71‑1 (MQ=255)
            gATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                                                   <  1:2515772/61‑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)
                  tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                                                   >  1:356224/1‑71 (MQ=255)
                  tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                                                   <  1:2255860/71‑1 (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)
                           ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                                                   <  1:1027894/46‑1 (MQ=255)
                           ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGt                                          >  1:197004/1‑71 (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)
                                  aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                                                                   <  1:2231458/39‑1 (MQ=255)
                                  aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                                 <  1:2023526/41‑1 (MQ=255)
                                  aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTa                                       <  1:3277672/67‑1 (MQ=255)
                                   ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                                 >  1:1277890/1‑40 (MQ=255)
                                   ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                                 <  1:2592923/40‑1 (MQ=255)
                                   ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                  <  1:1446602/71‑1 (MQ=255)
                                   ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                  >  1:1638252/1‑71 (MQ=255)
                                   ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                  <  1:514100/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)
                                                  atcATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCg                   >  1:2212166/3‑71 (MQ=255)
                                                    cATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGt                  <  1:3314130/70‑1 (MQ=255)
                                                      tCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTc               <  1:1419141/71‑1 (MQ=255)
                                                        cccATGCCGACGCTATTCATACGGCGGTCTGGGTGAGCCTGTTTACTCTGCATGtt                            >  1:3788282/1‑56 (MQ=255)
                                                        cccATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCt              <  1:708545/70‑1 (MQ=255)
                                                         tcATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGt                             <  1:3151205/53‑1 (MQ=255)
                                                         ccATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTc               <  1:939965/68‑1 (MQ=255)
                                                          cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGttt                >  1:4115960/1‑66 (MQ=255)
                                                          cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:3910992/1‑71 (MQ=255)
                                                          cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:1575045/1‑71 (MQ=255)
                                                          cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:1727044/1‑71 (MQ=255)
                                                          cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt            >  1:2191511/1‑70 (MQ=255)
                                                          aaTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:502231/2‑71 (MQ=255)
                                                           aTGCCGACGCTATTCATACGGCGGGCTGGGTTAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:194766/1‑70 (MQ=255)
                                                           aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGt                             <  1:1738963/52‑1 (MQ=255)
                                                           aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg          <  1:633390/71‑1 (MQ=255)
                                                           aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg          <  1:427964/71‑1 (MQ=255)
                                                            tGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg           >  1:1332663/1‑69 (MQ=255)
                                                              ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGttt                >  1:1139567/1‑62 (MQ=255)
                                                              ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc       >  1:210277/1‑71 (MQ=255)
                                                              ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc       >  1:2339690/1‑71 (MQ=255)
                                                               cGACTCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCt      <  1:2711162/71‑1 (MQ=255)
                                                                gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa                                >  1:2551776/1‑44 (MQ=255)
                                                                gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg     <  1:2407866/71‑1 (MQ=255)
                                                                gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg     <  1:225772/71‑1 (MQ=255)
                                                                 aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc       >  1:2201061/1‑68 (MQ=255)
                                                                 aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc       >  1:1413131/1‑68 (MQ=255)
                                                                 aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc       >  1:3326168/1‑68 (MQ=255)
                                                                 aCGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg     <  1:664420/70‑1 (MQ=255)
                                                                 aCGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg     <  1:2047322/70‑1 (MQ=255)
                                                                   gCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtgag         <  1:2844112/64‑3 (MQ=255)
                                                                   gCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGc  >  1:1555806/1‑71 (MQ=255)
                                                                   gCTATTCATACGGCGGGCTGGCGGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtgag         <  1:811604/64‑3 (MQ=255)
                                                                     tATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGt                             <  1:599802/42‑1 (MQ=255)
                                                                     |                                                                    
TGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGC  >  minE/2569790‑2569927

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: