Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,069 A→G 6.2% D458D (GAT→GAC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0690AG 6.2% 323.9 / 10.6 112D458D (GAT→GAClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (54/51);  new base G (2/5);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 113 total aligned reads displayed.

CGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAT  >  minE/1933001‑1933139
                                                                    |                                                                      
cGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCAt                                                                       <  1:3354042/70‑1 (MQ=255)
cGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATc                                                                      >  1:2674786/1‑71 (MQ=255)
     gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                                                  <  1:3042518/70‑1 (MQ=255)
      aTAGCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGaa                                                                <  1:2580252/71‑1 (MQ=255)
      aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                                                <  1:2808897/71‑1 (MQ=255)
      aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                                                <  1:2234781/71‑1 (MQ=255)
        aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                              >  1:1505129/1‑71 (MQ=255)
         ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                              >  1:546528/1‑70 (MQ=255)
         ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                              <  1:3797704/70‑1 (MQ=255)
         ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                             <  1:4037924/71‑1 (MQ=255)
         ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                             <  1:377206/71‑1 (MQ=255)
         ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                             <  1:3956211/71‑1 (MQ=255)
            aTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                              >  1:2290176/1‑67 (MQ=255)
             ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                                            <  1:436537/68‑1 (MQ=255)
             ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                                                         >  1:358568/1‑71 (MQ=255)
              tCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCtt                                                        >  1:1759730/1‑71 (MQ=255)
               cccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                       >  1:2614259/1‑71 (MQ=255)
                ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                       <  1:3279410/70‑1 (MQ=255)
                ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                       <  1:3062278/70‑1 (MQ=255)
                 cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                                                     >  1:3439755/1‑71 (MQ=255)
                 cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                                                     >  1:2620179/1‑71 (MQ=255)
                  aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                                                    <  1:1029968/71‑1 (MQ=255)
                  aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                                                    <  1:3570969/71‑1 (MQ=255)
                    aCCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                                                  <  1:2796468/71‑1 (MQ=255)
                     ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtcgg                                                  >  1:2149095/1‑69 (MQ=255)
                     ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                                                  >  1:1793098/1‑70 (MQ=255)
                     ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCgtcgtc                                                 <  1:2487568/71‑1 (MQ=255)
                      ctgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCgtcgtc                                                 <  1:109452/70‑1 (MQ=255)
                        gctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAcc                                              >  1:3228973/1‑71 (MQ=255)
                        gctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAcc                                              >  1:2749653/1‑71 (MQ=255)
                         ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                                                  >  1:487836/1‑50 (MQ=255)
                         ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGaa                                                                >  1:2553085/1‑52 (MQ=255)
                         ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                                               >  1:2167631/1‑69 (MQ=255)
                            tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                              >  1:925193/1‑51 (MQ=255)
                            tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                                                         <  1:1387151/56‑1 (MQ=255)
                            tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                                                         <  1:224303/56‑1 (MQ=255)
                              ggcggTAAGCCATGTTACAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                        >  1:3760169/1‑71 (MQ=255)
                                  gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                                                  <  1:3008451/41‑1 (MQ=255)
                                  gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                                            <  1:1312081/47‑1 (MQ=255)
                                  gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                                          <  1:1282271/65‑1 (MQ=255)
                                    aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                             >  1:3995179/1‑44 (MQ=255)
                                    aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgc                                  >  1:3583245/1‑71 (MQ=255)
                                     aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                                                  >  1:371267/1‑54 (MQ=255)
                                     aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                                                  >  1:2092251/1‑54 (MQ=255)
                                     aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                 >  1:311101/1‑71 (MQ=255)
                                       ccATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc                               <  1:2123592/71‑1 (MQ=255)
                                        cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                       <  1:996445/46‑1 (MQ=255)
                                        cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                              <  1:152277/71‑1 (MQ=255)
                                        cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                              <  1:3625449/71‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                   <  1:2522804/64‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                   <  1:1349057/64‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                   <  1:994455/64‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCaa                             <  1:1067070/70‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                              <  1:1929935/69‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            <  1:1771479/71‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            <  1:363002/71‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            <  1:3492382/71‑1 (MQ=255)
                                          tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            >  1:1174749/1‑71 (MQ=255)
                                            ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                                >  1:274670/1‑49 (MQ=255)
                                            ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            >  1:1757711/1‑69 (MQ=255)
                                            ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            >  1:1157642/1‑69 (MQ=255)
                                              ccAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                                <  1:1385538/47‑1 (MQ=255)
                                              ccAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTg                         <  1:598391/70‑1 (MQ=255)
                                                aGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                      >  1:782468/1‑71 (MQ=255)
                                                 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTTGCGGCAACCTGATcc                     <  1:3507696/71‑1 (MQ=255)
                                                 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                              <  1:1720469/62‑1 (MQ=255)
                                                 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                     <  1:470044/71‑1 (MQ=255)
                                                 gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                      >  1:3010349/1‑70 (MQ=255)
                                                      gTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                              <  1:1574322/57‑1 (MQ=255)
                                                        gACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                      >  1:1316063/1‑63 (MQ=255)
                                                         aCGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                   >  1:4024266/1‑47 (MQ=255)
                                                         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg             >  1:3797359/1‑71 (MQ=255)
                                                         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg             >  1:3104385/1‑71 (MQ=255)
                                                         aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg             <  1:515665/71‑1 (MQ=255)
                                                          cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGt            <  1:2941025/71‑1 (MQ=255)
                                                          cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt            <  1:2302727/71‑1 (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:1938148/1‑71 (MQ=255)
                                                             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc         >  1:2004059/1‑71 (MQ=255)
                                                             aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc         >  1:3215066/1‑71 (MQ=255)
                                                              aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                          >  1:2074084/1‑53 (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)
                                                                cTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            >  1:1859216/1‑49 (MQ=255)
                                                                 ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt     >  1:1560276/1‑71 (MQ=255)
                                                                 ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                            <  1:1825225/48‑1 (MQ=255)
                                                                 ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                   >  1:194897/1‑57 (MQ=255)
                                                                  tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc                               >  1:127367/1‑44 (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)
                                                                    |                                                                      
CGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAT  >  minE/1933001‑1933139

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: