Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,117 A→T 8.4% D442E (GAT→GAA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1170AT 8.4% 231.4 / 11.8 83D442E (GAT→GAAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (49/27);  new base T (2/5);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCAC  >  minE/1933049‑1933181
                                                                    |                                                                
aGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                >  1:782468/1‑71 (MQ=255)
 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTTGCGGCAACCTGATcc                                                               <  1:3507696/71‑1 (MQ=255)
 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                               <  1:470044/71‑1 (MQ=255)
 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                               <  1:1430577/71‑1 (MQ=255)
 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                >  1:3010349/1‑70 (MQ=255)
        gACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                >  1:1316063/1‑63 (MQ=255)
         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                       >  1:3104385/1‑71 (MQ=255)
         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                       <  1:515665/71‑1 (MQ=255)
         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                       >  1:3797359/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)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGa                                                                  >  1:699377/1‑56 (MQ=255)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                >  1:3720045/1‑58 (MQ=255)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                >  1:1490997/1‑58 (MQ=255)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                   >  1:3215066/1‑71 (MQ=255)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                   >  1:1938148/1‑71 (MQ=255)
             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                   >  1:2004059/1‑71 (MQ=255)
              aTCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGttcag                                                  <  1:2250979/71‑4 (MQ=255)
              aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                      >  1:3357430/1‑67 (MQ=255)
              aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                  >  1:343493/1‑71 (MQ=255)
              aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                  >  1:1399802/1‑71 (MQ=255)
              aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                  >  1:2775405/1‑71 (MQ=255)
                 ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                               >  1:1560276/1‑71 (MQ=255)
                 ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                                                             >  1:194897/1‑57 (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)
                     tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                                <  1:229421/50‑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)
                                 gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                               <  1:2317202/39‑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:578901/1‑70 (MQ=255)
                                            aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                     <  1:3380628/70‑1 (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:3474159/70‑1 (MQ=255)
                                                 cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                <  1:295992/70‑1 (MQ=255)
                                                 cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                <  1:1423393/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:2465105/1‑71 (MQ=255)
                                                     tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc           >  1:3077415/1‑71 (MQ=255)
                                                     tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc           >  1:1606958/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:664871/70‑1 (MQ=255)
                                                       tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc          >  1:1866890/1‑70 (MQ=255)
                                                       tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc          >  1:3708823/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)
                                                                    |                                                                
AGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCAC  >  minE/1933049‑1933181

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