Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,847 C→T 5.9% P39P (CCC→CCT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8470CT 5.9% 277.2 / 5.1 84P39P (CCC→CCTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (41/38);  new base T (3/2);  total (44/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 = 7.75e-01

GTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCT  >  minE/2569779‑2569915
                                                                    |                                                                    
gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcccc                                                                     <  1:37745/70‑1 (MQ=255)
gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc                                                                      <  1:4026137/69‑1 (MQ=255)
 tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                                   >  1:2614097/1‑71 (MQ=255)
 tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                                   >  1:3290944/1‑71 (MQ=255)
 tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                                   <  1:621373/71‑1 (MQ=255)
 tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                                                   <  1:1401956/71‑1 (MQ=255)
  gATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAATCATcccc                                                                     >  1:1274487/1‑68 (MQ=255)
  gATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcccc                                                                     >  1:2158427/1‑68 (MQ=255)
    ttttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                                <  1:1132263/71‑1 (MQ=255)
     tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                                                >  1:3632706/1‑70 (MQ=255)
          aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                           <  1:1134066/70‑1 (MQ=255)
          aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                           >  1:3986902/1‑70 (MQ=255)
          aTGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGc                                                           >  1:462296/1‑70 (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:4057028/70‑1 (MQ=255)
             gAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAt                                                        <  1:262934/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)
                tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                     >  1:3473913/1‑70 (MQ=255)
                tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCAt                                                    <  1:3369987/71‑1 (MQ=255)
                 catcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc                                                                      >  1:1370937/1‑52 (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)
                             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:4085203/1‑51 (MQ=255)
                                   cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                     >  1:1407391/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:930878/1‑37 (MQ=255)
                                        aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                              >  1:2277523/1‑37 (MQ=255)
                                        aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGa                                                              >  1:2969052/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:1917940/1‑70 (MQ=255)
                                           tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc                          >  1:1228502/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:1277890/1‑40 (MQ=255)
                                              ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                                                     <  1:2592923/40‑1 (MQ=255)
                                              ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                      <  1:514100/71‑1 (MQ=255)
                                              ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                      <  1:1446602/71‑1 (MQ=255)
                                              ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                      >  1:1638252/1‑71 (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)
                                                                    ccATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTc   <  1:939965/68‑1 (MQ=255)
                                                                    |                                                                    
GTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCT  >  minE/2569779‑2569915

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: