Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,050,271 A→G 11.6% L28L (CTT→CTC glnE ← fused deadenylyltransferase and adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,050,2710AG11.6% 182.6 / 12.0 69L28L (CTT→CTCglnEfused deadenylyltransferase and adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (3/5);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACG  >  minE/2050202‑2050335
                                                                     |                                                                
ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTaa                                                                 <  1:772090/71‑1 (MQ=255)
ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTaa                                                                 <  1:3391693/71‑1 (MQ=255)
 tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAg                                                                >  1:100386/1‑71 (MQ=255)
 tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAg                                                                >  1:619103/1‑71 (MQ=255)
     gATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGt                                                               <  1:576242/68‑1 (MQ=255)
     gATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGa                                                             >  1:3764625/1‑70 (MQ=255)
     gATGCGCAATCACGCTGTCATGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTGAGt                                                               <  1:1419301/68‑1 (MQ=255)
       tGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTc                                                          <  1:3033642/71‑1 (MQ=255)
         cgcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTcc                                                         <  1:283945/70‑1 (MQ=255)
         cgcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                        <  1:2695980/71‑1 (MQ=255)
         cgcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                        <  1:1485093/71‑1 (MQ=255)
          gcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                       >  1:1457580/1‑71 (MQ=255)
           cAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                       <  1:2756460/70‑1 (MQ=255)
           cAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCg                                                      <  1:1348576/71‑1 (MQ=255)
            aaTCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCgg                                                     >  1:229823/1‑71 (MQ=255)
              tCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGc                                                    >  1:302604/1‑70 (MQ=255)
              acACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGc                                                    >  1:3591640/2‑70 (MQ=255)
                 cGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTaaa                                                <  1:771804/71‑1 (MQ=255)
                  gCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAg                                               >  1:195551/1‑71 (MQ=255)
                      tCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                           >  1:2399702/1‑71 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTa                                                                  <  1:752157/44‑1 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGctc                                           >  1:3280140/1‑67 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTg                                                              >  1:3167678/1‑48 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTg                                                              >  1:3263152/1‑48 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTgg                                        >  1:1391159/1‑70 (MQ=255)
                          gCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTgg                                        >  1:2249431/1‑70 (MQ=255)
                           cacaAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                        <  1:2820736/53‑1 (MQ=255)
                           cacaAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGc                                       <  1:928263/70‑1 (MQ=255)
                             caAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                    <  1:609910/71‑1 (MQ=255)
                              aaaaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                   >  1:2545204/1‑71 (MQ=255)
                              aaaaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                   >  1:699289/1‑71 (MQ=255)
                              aaaaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                    >  1:704041/1‑70 (MQ=255)
                               aaaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTGAGTGATTc                                                          <  1:3638997/47‑1 (MQ=255)
                               aaaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                   <  1:2038484/70‑1 (MQ=255)
                                aaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCg                                  >  1:1497799/1‑70 (MQ=255)
                                aaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGc                                 <  1:3293650/71‑1 (MQ=255)
                                aaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGc                                 >  1:1746218/1‑71 (MQ=255)
                                aaTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGc                                 <  1:1698481/71‑1 (MQ=255)
                                 aTCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCt                                <  1:141072/71‑1 (MQ=255)
                                  tCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTc                                                          <  1:1958948/44‑1 (MQ=255)
                                  tCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTc                                                          <  1:2597044/44‑1 (MQ=255)
                                   cACTAAAAGTAAGTACTGACTTAGCCTGTGCGCTAAGTGATTCCTCGGCTaaa                                                >  1:1068793/1‑53 (MQ=255)
                                     cTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaac                            <  1:344358/71‑1 (MQ=255)
                                      tAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGt                                                               <  1:292202/35‑1 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCaacaac                         >  1:2762493/1‑71 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                        >  1:667247/1‑40 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                    >  1:2765142/1‑60 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                    >  1:1226133/1‑60 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaacaa                          <  1:2401014/70‑1 (MQ=255)
                                        aaaGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaaca                           <  1:2449733/69‑1 (MQ=255)
                                         aagtaagtACTGACTTCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACg                        <  1:64636/71‑1 (MQ=255)
                                         aagtaagtACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACg                        <  1:2984599/71‑1 (MQ=255)
                                             aagtaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAg                                     >  1:1680665/1‑54 (MQ=255)
                                             aagtaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                     <  1:124719/70‑1 (MQ=255)
                                              agtaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACg                        >  1:2393610/1‑66 (MQ=255)
                                              agtaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCt                    >  1:2493331/1‑70 (MQ=255)
                                                taCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                  >  1:811226/1‑70 (MQ=255)
                                                taCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                  >  1:3859748/1‑70 (MQ=255)
                                                taCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                  >  1:2691137/1‑70 (MQ=255)
                                                  cTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                   <  1:815599/51‑1 (MQ=255)
                                                   tGACTTCGCCTGTGCGCTAAGTGATTCCTCTGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCt                    <  1:1314841/65‑1 (MQ=255)
                                                     aCTTCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc            <  1:2502602/71‑1 (MQ=255)
                                                     aCTTCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc            <  1:2497974/71‑1 (MQ=255)
                                                     aCTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa             <  1:13359/70‑1 (MQ=255)
                                                       ttCGCCTGTGCGCTGAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg          >  1:382593/1‑71 (MQ=255)
                                                       ttCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg          >  1:449724/1‑71 (MQ=255)
                                                         cGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgct        >  1:1030492/1‑71 (MQ=255)
                                                         cGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgct        >  1:1267234/1‑71 (MQ=255)
                                                         cGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg          <  1:1187611/69‑1 (MQ=255)
                                                          gCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg       >  1:3702292/1‑71 (MQ=255)
                                                            cTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg               <  1:1767293/61‑1 (MQ=255)
                                                               tgCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACg  <  1:929777/71‑1 (MQ=255)
                                                                     |                                                                
CTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACG  >  minE/2050202‑2050335

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: