Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,126 G→A 9.9% R439R (CGC→CGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1260GA 9.9% 243.0 / 12.7 81R439R (CGC→CGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (47/26);  new base A (3/5);  total (50/31)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAG  >  minE/1933058‑1933195
                                                                    |                                                                     
aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                     <  1:515665/71‑1 (MQ=255)
aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                     >  1:3797359/1‑71 (MQ=255)
aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                     >  1:3104385/1‑71 (MQ=255)
 cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGt                                                                    <  1:2941025/71‑1 (MQ=255)
 cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                    <  1:2302727/71‑1 (MQ=255)
 cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                    >  1:3679064/1‑71 (MQ=255)
  gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGtt                                                                   <  1:381791/71‑1 (MQ=255)
  gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGtt                                                                   >  1:3037191/1‑71 (MQ=255)
  gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                    <  1:3778205/70‑1 (MQ=255)
    aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                 >  1:1938148/1‑71 (MQ=255)
    aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                 >  1:2004059/1‑71 (MQ=255)
    aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                 >  1:3215066/1‑71 (MQ=255)
     aTCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGttcag                                                                <  1:2250979/71‑4 (MQ=255)
     aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                    >  1:3357430/1‑67 (MQ=255)
     aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                >  1:2775405/1‑71 (MQ=255)
     aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                >  1:343493/1‑71 (MQ=255)
     aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                >  1:1399802/1‑71 (MQ=255)
        ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                                             >  1:1560276/1‑71 (MQ=255)
         tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                                             >  1:3326727/1‑70 (MQ=255)
         tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                  >  1:1112987/1‑65 (MQ=255)
         tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                  >  1:1108282/1‑65 (MQ=255)
          cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGcc                                                              <  1:1098252/68‑1 (MQ=255)
           aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                                          >  1:2998414/1‑71 (MQ=255)
            tcatAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                                             <  1:3612344/64‑1 (MQ=255)
            tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                 >  1:3085872/1‑63 (MQ=255)
            tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                  >  1:878957/1‑62 (MQ=255)
            tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                                         <  1:2753595/71‑1 (MQ=255)
             cGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGAGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCt                                                        >  1:2037954/1‑71 (MQ=255)
              atAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCtt                                                       <  1:1614705/70‑1 (MQ=255)
                 gAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGAt                                                          >  1:507916/1‑65 (MQ=255)
                    cATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                                                 <  1:1095778/71‑1 (MQ=255)
                        gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAtt                                             >  1:133817/1‑71 (MQ=255)
                        gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAt                                              <  1:4038586/70‑1 (MQ=255)
                          ttCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTc                                                      >  1:511336/1‑60 (MQ=255)
                             gtcgtcACCTGCGTCTTTCGCGGCAACATGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt                                        >  1:2566904/1‑71 (MQ=255)
                              tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                     <  1:3866983/41‑1 (MQ=255)
                              tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTc                                       >  1:1773734/1‑71 (MQ=255)
                               cgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                                                 <  1:1239521/60‑1 (MQ=255)
                                gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAg                                          >  1:3465021/1‑66 (MQ=255)
                                gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGc                                     >  1:2524913/1‑71 (MQ=255)
                                  cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                                             >  1:4053302/1‑45 (MQ=255)
                                  cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCt                                    >  1:113009/1‑70 (MQ=255)
                                   aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTaa                                  <  1:144523/71‑1 (MQ=255)
                                   aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                   <  1:3380628/70‑1 (MQ=255)
                                   aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                   >  1:578901/1‑70 (MQ=255)
                                     cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTg                                                            <  1:1579981/43‑1 (MQ=255)
                                     cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                                >  1:3917514/1‑71 (MQ=255)
                                     cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                                >  1:3774672/1‑71 (MQ=255)
                                      tGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc                               >  1:4093441/1‑71 (MQ=255)
                                       gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              >  1:1709237/1‑71 (MQ=255)
                                       gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              >  1:2384922/1‑71 (MQ=255)
                                        cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              <  1:1423393/70‑1 (MQ=255)
                                        cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              <  1:295992/70‑1 (MQ=255)
                                        cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              <  1:3474159/70‑1 (MQ=255)
                                        cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc                               <  1:2446426/69‑1 (MQ=255)
                                        cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCg                             >  1:2379679/1‑71 (MQ=255)
                                          tCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGt                            >  1:2540938/1‑70 (MQ=255)
                                            tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                         >  1:1407380/1‑71 (MQ=255)
                                            tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                         >  1:3740954/1‑71 (MQ=255)
                                            tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                         >  1:1606958/1‑71 (MQ=255)
                                            tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                         >  1:3077415/1‑71 (MQ=255)
                                            tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                         >  1:2465105/1‑71 (MQ=255)
                                             ttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAg                                          <  1:3769836/53‑1 (MQ=255)
                                             ttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                   <  1:1485610/60‑1 (MQ=255)
                                              tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTaga                        >  1:970288/1‑68 (MQ=255)
                                              tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                        <  1:3442265/70‑1 (MQ=255)
                                              tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                        >  1:3708823/1‑70 (MQ=255)
                                              tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                        <  1:664871/70‑1 (MQ=255)
                                              tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                        >  1:1866890/1‑70 (MQ=255)
                                                   gCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                              <  1:1474792/59‑1 (MQ=255)
                                                     aaCCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCAc                <  1:1239501/71‑1 (MQ=255)
                                                          gATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt                                        >  1:1557447/1‑42 (MQ=255)
                                                           aTCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGc                                         >  1:3046573/1‑40 (MQ=255)
                                                            tCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGaccacc         <  1:1941729/71‑1 (MQ=255)
                                                                ggAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAAt      >  1:858607/1‑70 (MQ=255)
                                                                 gAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGa              >  1:942031/1‑61 (MQ=255)
                                                                   aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAGAGCCGTTACCGATTTCACGACCACCAATgaag  <  1:1008066/71‑1 (MQ=255)
                                                                   aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATtaag  <  1:949388/71‑1 (MQ=255)
                                                                   aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATgaag  <  1:2645024/71‑1 (MQ=255)
                                                                   aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATg     >  1:1155013/1‑68 (MQ=255)
                                                                   aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATg     >  1:699282/1‑68 (MQ=255)
                                                                    |                                                                     
ACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAG  >  minE/1933058‑1933195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: