Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,264,986 T→C 17.5% E15G (GAA→GGA)  secE ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,264,9860TC17.5% 116.0 / 17.7 54E15G (GAA→GGA) secEpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base C (4/5);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTA  >  minE/2264919‑2265054
                                                                   |                                                                    
gcgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTcc                                                                    <  1:2233129/70‑1 (MQ=255)
gcgATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTc                                                                     >  1:3756059/1‑69 (MQ=255)
    tAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGcc                                                               >  1:1456226/1‑71 (MQ=255)
        gATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCgcgc                                                           <  1:195942/71‑1 (MQ=255)
         aTAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGcc                                                          >  1:1959122/1‑71 (MQ=255)
          tAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAgg                                                                 <  1:3849324/63‑1 (MQ=255)
          tAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGccc                                                         >  1:1392793/1‑71 (MQ=255)
             ttGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGccc                                                         <  1:3127290/68‑1 (MQ=255)
             ttGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCt                                                      <  1:6051/71‑1 (MQ=255)
             ttGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCt                                                      <  1:3332946/71‑1 (MQ=255)
               gCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTc                                                    >  1:271491/1‑71 (MQ=255)
                ccGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTc                                                    <  1:2411535/70‑1 (MQ=255)
                  gACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                 <  1:2936108/71‑1 (MQ=255)
                   aCAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                 >  1:1978140/1‑70 (MQ=255)
                    cAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                 >  1:2564499/1‑69 (MQ=255)
                       tCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTcc                                                                    <  1:1345801/47‑1 (MQ=255)
                       tCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                 <  1:1809272/66‑1 (MQ=255)
                              caGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGcc                                                               <  1:2821995/45‑1 (MQ=255)
                              caGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTa                                      >  1:2466121/1‑70 (MQ=255)
                               aGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                                    <  1:2604789/71‑1 (MQ=255)
                                ggAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTc                                   >  1:2453430/1‑71 (MQ=255)
                                  aGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGccc                                                         >  1:3794598/1‑47 (MQ=255)
                                  aGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGccc                                                         >  1:64020/1‑47 (MQ=255)
                                   gCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                                <  1:2569683/71‑1 (MQ=255)
                                     aaTGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACt                              <  1:281973/71‑1 (MQ=255)
                                         ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                                >  1:2410870/1‑65 (MQ=255)
                                         ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATa                          <  1:1282664/71‑1 (MQ=255)
                                          caccacAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATaa                         <  1:1231759/71‑1 (MQ=255)
                                           accacAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATaaa                        <  1:2355267/71‑1 (MQ=255)
                                            ccacAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATa                          <  1:2645782/68‑1 (MQ=255)
                                              acaACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCt                                            <  1:1397814/48‑1 (MQ=255)
                                              acaACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaacc                      >  1:2725359/1‑70 (MQ=255)
                                              acaACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTc                                   >  1:597199/1‑57 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCAc                               <  1:2809823/60‑1 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATaaa                        >  1:1384511/1‑67 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaa                    >  1:1952259/1‑71 (MQ=255)
                                               caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTAGGTATTCGCACTCATAaaccaa                    >  1:356347/1‑71 (MQ=255)
                                                aaCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTg                                                <  1:3670524/42‑1 (MQ=255)
                                                aaCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCCTCCTTGAGCTTCGGTATTCGCACTCATAaaccaac                   >  1:559705/1‑71 (MQ=255)
                                                  cGACCCACTTCATCGCTACCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAt                                     <  1:624958/51‑1 (MQ=255)
                                                   gACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTg                <  1:1635080/71‑1 (MQ=255)
                                                   gACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTg                <  1:57384/71‑1 (MQ=255)
                                                     cccACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                                    >  1:343939/1‑49 (MQ=255)
                                                      ccACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCtt                                           <  1:574151/41‑1 (MQ=255)
                                                      ccACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaac                   >  1:1848967/1‑65 (MQ=255)
                                                       cACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                                    <  1:2989762/47‑1 (MQ=255)
                                                          ttCATCGCTTCCAGGCCTCGCCCGCTTCCTTGAGCTTCGGTAtt                                    <  1:3104158/44‑1 (MQ=255)
                                                           tCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTc                                          >  1:313015/1‑37 (MQ=255)
                                                            cATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt       >  1:2864868/1‑71 (MQ=255)
                                                            cATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt       >  1:1644528/1‑71 (MQ=255)
                                                             aTCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGttt      >  1:120098/1‑71 (MQ=255)
                                                              tCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                            <  1:3486205/48‑1 (MQ=255)
                                                              tCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                            <  1:2313962/48‑1 (MQ=255)
                                                                 cTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTa  >  1:1403646/1‑71 (MQ=255)
                                                                   |                                                                    
GCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTA  >  minE/2264919‑2265054

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: