Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,341,041 C→T 9.2% A41T (GCG→ACG)  wcaB ← predicted acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,341,0410CT 9.2% 155.2 / 4.6 54A41T (GCG→ACG) wcaBpredicted acyl transferase
Reads supporting (aligned to +/- strand):  ref base C (32/17);  new base T (2/3);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAAC  >  minE/1340972‑1341110
                                                                     |                                                                     
cgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGgccgcc                                                                      <  1:1128082/71‑1 (MQ=255)
    aGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGTCca                                                                    >  1:3000089/1‑69 (MQ=255)
         ggATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCca                                                                    >  1:2810851/1‑64 (MQ=255)
         ggATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGAttg                                                             >  1:1635396/1‑71 (MQ=255)
         ggATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGAttg                                                             >  1:2863009/1‑71 (MQ=255)
         ggATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGAttg                                                             >  1:2959170/1‑71 (MQ=255)
           aTTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGAttgtt                                                           <  1:2845087/71‑1 (MQ=255)
             ttCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGccc                                                                     >  1:1884652/1‑59 (MQ=255)
             ttCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGAttgtt                                                           >  1:338496/1‑69 (MQ=255)
              tCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAg                                                        >  1:17072/1‑71 (MQ=255)
              tCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAg                                                        >  1:1415351/1‑71 (MQ=255)
              tCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAg                                                        >  1:2367434/1‑71 (MQ=255)
               cATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCcaca                                                                  <  1:1983058/60‑1 (MQ=255)
               cATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAg                                                        >  1:2791130/1‑70 (MQ=255)
                aTAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGa                                                      >  1:3950227/1‑71 (MQ=255)
                  aaCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCTGGGCCGCCCACAGATTGTTGAGGAc                                                     >  1:2285820/1‑70 (MQ=255)
                  aaCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACg                                                    >  1:931588/1‑71 (MQ=255)
                   aCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGTCCACAGATTGTTGAGGACGt                                                   >  1:795322/1‑71 (MQ=255)
                   aCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGa                                                      <  1:550278/68‑1 (MQ=255)
                   aCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGt                                                   >  1:3482903/1‑71 (MQ=255)
                   aCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGt                                                   <  1:3251640/71‑1 (MQ=255)
                   aCCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGt                                                   >  1:3019548/1‑71 (MQ=255)
                     cGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTc                                                 >  1:744928/1‑71 (MQ=255)
                      gAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCt                                                >  1:251125/1‑71 (MQ=255)
                        aaaaGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGTCCACAGATTGTTGAGGACGTTCttt                                              <  1:694125/71‑1 (MQ=255)
                           aGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTg                                           <  1:2565183/71‑1 (MQ=255)
                            gCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTgc                                          <  1:916125/71‑1 (MQ=255)
                            gCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTgc                                          <  1:3811508/71‑1 (MQ=255)
                             cATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTgcg                                         >  1:2981222/1‑71 (MQ=255)
                                tCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGTCCACAGATTGTTGAGGACGTTCttt                                              <  1:3599104/63‑1 (MQ=255)
                                tCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGcc                                       >  1:591596/1‑70 (MQ=255)
                                tCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCca                                      >  1:3671804/1‑71 (MQ=255)
                                 cGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCtttt                                             >  1:218250/1‑63 (MQ=255)
                                      ataataCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCg                                 >  1:2293329/1‑70 (MQ=255)
                                      aTAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCacac                                   <  1:3381759/68‑1 (MQ=255)
                                      aTAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCg                                 <  1:2061262/70‑1 (MQ=255)
                                        aaTGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGa                                                         >  1:2874163/1‑44 (MQ=255)
                                         aTGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACa                             <  1:3933950/71‑1 (MQ=255)
                                          tGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTg                                           >  1:1906246/1‑56 (MQ=255)
                                           gCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaa                            >  1:2309560/1‑70 (MQ=255)
                                            cGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaa                           <  1:1859448/70‑1 (MQ=255)
                                            cGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaa                           >  1:1930727/1‑70 (MQ=255)
                                             gATACAGCACCCGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaaa                         >  1:3843166/1‑71 (MQ=255)
                                             gATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGt                                                   <  1:2169418/45‑1 (MQ=255)
                                              aTACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAAt                        <  1:2982762/71‑1 (MQ=255)
                                                 cAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAgg                                                       >  1:3909611/1‑37 (MQ=255)
                                                  aGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCttttt                                            >  1:3930550/1‑47 (MQ=255)
                                                  aGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCttttt                                            >  1:198718/1‑47 (MQ=255)
                                                    cACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCaa                  >  1:116428/1‑71 (MQ=255)
                                                      ccAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAaca                <  1:1222376/71‑1 (MQ=255)
                                                       cagcagCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTgcgc                                        <  1:3394471/46‑1 (MQ=255)
                                                       cagcagCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCAGACCGAACa                             <  1:2630294/57‑1 (MQ=255)
                                                         gcagcGGGGCCGTCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa             <  1:2151065/71‑1 (MQ=255)
                                                          cagcGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa             >  1:838008/1‑70 (MQ=255)
                                                          cagcGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt            <  1:2536519/71‑1 (MQ=255)
                                                              ggggCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAgg         >  1:190002/1‑70 (MQ=255)
                                                                 gccgccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt            <  1:3576936/64‑1 (MQ=255)
                                                                  ccgccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGa    <  1:230993/71‑1 (MQ=255)
                                                                    gccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAAt                        >  1:184924/1‑49 (MQ=255)
                                                                     cccACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAAc  >  1:527498/1‑70 (MQ=255)
                                                                     |                                                                     
CGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAAC  >  minE/1340972‑1341110

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: