Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,072 G→A 27.8% Y457Y (TAC→TAT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0720GA27.8% 352.4 / 147.7 219Y457Y (TAC→TATlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (87/71);  new base A (38/23);  total (125/94)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 230 total aligned reads displayed.

GGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGA  >  minE/1933002‑1933138
                                                                      |                                                                  
  tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:2170557/1‑69 (MQ=255)
  tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:2556433/1‑69 (MQ=255)
  tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:3678290/1‑69 (MQ=255)
  tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:3344681/1‑69 (MQ=255)
   cGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:4360687/1‑68 (MQ=255)
   cGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:1614062/1‑68 (MQ=255)
   cGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    >  1:393726/1‑68 (MQ=255)
    gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTa                                                                  >  1:177889/1‑69 (MQ=255)
    gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                                                 >  1:585365/1‑70 (MQ=255)
    gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                                                 >  1:3623786/1‑70 (MQ=255)
     aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                                               <  1:2935772/71‑1 (MQ=255)
     aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                                               <  1:2874446/71‑1 (MQ=255)
       aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                             >  1:3031612/1‑71 (MQ=255)
        ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGa                                                                <  1:203293/67‑1 (MQ=255)
        ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                                             <  1:5325187/70‑1 (MQ=255)
        ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                            >  1:243506/1‑71 (MQ=255)
        ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                                            >  1:787249/1‑71 (MQ=255)
        ccAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACa                                                             <  1:2339575/70‑1 (MQ=255)
         cAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                                           >  1:5193060/1‑71 (MQ=255)
         cAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                                           >  1:5018567/1‑71 (MQ=255)
           atacccAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGc                                                         >  1:5622597/4‑71 (MQ=255)
           aTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGc                                                         <  1:1390989/71‑1 (MQ=255)
           aTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGc                                                         >  1:4144127/1‑71 (MQ=255)
            ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                                           <  1:203776/68‑1 (MQ=255)
              cccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                      >  1:1049344/1‑71 (MQ=255)
              cccAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtc                                                      <  1:5442157/71‑1 (MQ=255)
               ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                      <  1:867354/70‑1 (MQ=255)
                cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                                                    >  1:3588923/1‑71 (MQ=255)
                cAGACCTGCTGTCGGCGGTAAGCCATGTTCAAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                                                    <  1:1257223/71‑1 (MQ=255)
                 aGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtcgtc                                                   <  1:642669/71‑1 (MQ=255)
                 aGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtcgtc                                                   >  1:2688439/1‑71 (MQ=255)
                  gACCTGCTGTCGGCGGTAAGCCATGTTGCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtcg                                                  >  1:2576131/1‑71 (MQ=255)
                    ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCgtcgtc                                                <  1:857411/71‑1 (MQ=255)
                     ctgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                               <  1:3368009/71‑1 (MQ=255)
                      tgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                                                    <  1:4293080/49‑1 (MQ=255)
                      tgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                                              >  1:3667897/1‑71 (MQ=255)
                          gTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                                                   >  1:4888566/1‑62 (MQ=255)
                          gTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                                                   >  1:4344229/1‑62 (MQ=255)
                             ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCCTCGTCGTCACCTGCGTc                                       >  1:2235179/1‑71 (MQ=255)
                              gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                      <  1:132381/71‑1 (MQ=255)
                                 gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAAc                                                              <  1:1524349/44‑1 (MQ=255)
                                 gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                                                      <  1:1561213/52‑1 (MQ=255)
                                 gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc                                   <  1:4807766/71‑1 (MQ=255)
                                  tAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                                                 <  1:2956462/56‑1 (MQ=255)
                                   aaGCCAAGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgc                                 >  1:2210741/1‑71 (MQ=255)
                                    aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                                              >  1:3049866/1‑57 (MQ=255)
                                    aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                                              >  1:5115191/1‑57 (MQ=255)
                                      ccATGTTCCAGTGCGGTGACATAATCTTCGTCGTAGAACATCGCTTCGTCGTCACCTGCATCt                                      <  1:1116696/63‑1 (MQ=255)
                                      ccATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtc                                                      <  1:1729618/47‑1 (MQ=255)
                                       cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                                                   >  1:2885321/1‑49 (MQ=255)
                                         tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                      <  1:3339026/60‑1 (MQ=255)
                                         tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                <  1:4510530/66‑1 (MQ=255)
                                           ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                         >  1:3531122/1‑71 (MQ=255)
                                           ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACa                                                             >  1:3405306/1‑35 (MQ=12)
                                           ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCtt                                                       >  1:2302530/1‑41 (MQ=16)
                                           ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                  >  1:3308526/1‑60 (MQ=255)
                                           ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                  >  1:1858649/1‑60 (MQ=255)
                                           ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                  >  1:2725835/1‑60 (MQ=255)
                                             ccAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgt                       >  1:1492064/1‑70 (MQ=255)
                                             ccAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgt                       >  1:4727515/1‑70 (MQ=255)
                                             ccAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgt                       >  1:3222538/1‑70 (MQ=255)
                                             ccAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgt                       >  1:1089950/1‑70 (MQ=255)
                                              cAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCtt                                                       >  1:546473/1‑38 (MQ=255)
                                              cAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCAc                                              >  1:5659632/1‑47 (MQ=255)
                                              cAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgtt                      <  1:691635/71‑3 (MQ=255)
                                              cAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTgtt                      <  1:2131891/71‑3 (MQ=255)
                                               aGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                  >  1:2415427/1‑58 (MQ=255)
                                                gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                               <  1:3693082/44‑1 (MQ=255)
                                                gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                               <  1:3699960/44‑1 (MQ=255)
                                                gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                                               <  1:3920011/44‑1 (MQ=255)
                                                 tGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCttt                                    <  1:1154704/54‑1 (MQ=255)
                                                 tGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                      <  1:4422532/52‑1 (MQ=255)
                                                 tGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                     >  1:1347867/1‑69 (MQ=255)
                                                     gTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                    <  1:208561/66‑1 (MQ=255)
                                                     gTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                    <  1:1990504/66‑1 (MQ=255)
                                                     gTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                    <  1:702811/66‑1 (MQ=255)
                                                       caCGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGaaac             >  1:3999744/2‑69 (MQ=255)
                                                        aCGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACg            >  1:108465/1‑71 (MQ=255)
                                                        aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGc             <  1:5499635/70‑1 (MQ=255)
                                                         cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCAc                                              >  1:2060578/1‑36 (MQ=25)
                                                         cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagc                                 <  1:4942817/49‑1 (MQ=255)
                                                         cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGt           <  1:344419/71‑1 (MQ=255)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                  >  1:551934/1‑45 (MQ=255)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagca                                >  1:3202087/1‑48 (MQ=19)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagca                                >  1:4627129/1‑48 (MQ=19)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagca                                >  1:5189370/1‑48 (MQ=19)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACg            <  1:251110/69‑1 (MQ=255)
                                                          gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGtt          >  1:2153822/1‑71 (MQ=255)
                                                          gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGc             <  1:5420870/68‑1 (MQ=255)
                                                           tAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGttc         <  1:145909/71‑2 (MQ=255)
                                                            aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt           <  1:2769100/68‑1 (MQ=255)
                                                             aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc        <  1:1217528/70‑1 (MQ=255)
                                                               cTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                       >  1:2245386/1‑37 (MQ=255)
                                                                ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGaaa              >  1:3260059/1‑60 (MQ=255)
                                                                ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                >  1:4230971/1‑43 (MQ=255)
                                                                ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc      >  1:5430545/1‑69 (MQ=255)
                                                                ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGcc     >  1:1685781/1‑70 (MQ=255)
                                                                  cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc      <  1:3314824/67‑1 (MQ=255)
                                                                  cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt    <  1:5634821/69‑1 (MQ=255)
                                                                  cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTg   <  1:3157432/70‑1 (MQ=255)
                                                                      |                                                                  
GGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGA  >  minE/1933002‑1933138

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: