Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,639 G→A 31.6% S268S (TCC→TCT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6390GA31.6% 227.1 / 134.1 171S268S (TCC→TCTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (54/63);  new base A (25/29);  total (79/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 173 total aligned reads displayed.

CGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACG  >  minE/1933569‑1933707
                                                                      |                                                                    
   taaGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                    <  1:2749264/68‑1 (MQ=255)
   tCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                    >  1:5617105/1‑70 (MQ=255)
     aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:4235832/1‑70 (MQ=255)
          gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                                                            >  1:4159908/4‑71 (MQ=255)
          gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                                                            >  1:5217119/4‑71 (MQ=255)
            tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCt                                                             <  1:4508666/68‑1 (MQ=255)
              tAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGt                                                         <  1:3745819/70‑1 (MQ=255)
                aTCCGCATACGCCATGTAGAGTTACATCATTGTGAACTCAGGATTATGGCGAACAGaa                                                                   >  1:4041847/1‑58 (MQ=255)
                   tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg                                                    >  1:3010245/1‑70 (MQ=255)
                        acgCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                    <  1:850492/47‑1 (MQ=255)
                        acgCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                    <  1:3607038/47‑1 (MQ=255)
                        aaGCCATGTAGATTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa                                              >  1:3969084/1‑71 (MQ=255)
                        aaGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGt                                                <  1:5298862/69‑1 (MQ=255)
                        aaGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa                                              >  1:2642069/1‑71 (MQ=255)
                         cgCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTa                                              <  1:2422989/69‑1 (MQ=255)
                         cgCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTa                                              <  1:1962743/69‑1 (MQ=255)
                         cgCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTa                                              <  1:2043105/69‑1 (MQ=255)
                         aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATa                                                                >  1:1549854/1‑52 (MQ=255)
                         aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa                                              <  1:125228/70‑1 (MQ=255)
                         aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc                                             <  1:2954998/71‑1 (MQ=255)
                         aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc                                             >  1:4836025/1‑71 (MQ=255)
                         aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc                                             <  1:489138/71‑1 (MQ=255)
                          gCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACg                                            <  1:1362355/71‑1 (MQ=255)
                          gCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACg                                            <  1:4122351/71‑1 (MQ=255)
                            cATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGtt                                               <  1:3032491/66‑1 (MQ=255)
                             aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:1890924/71‑1 (MQ=255)
                             aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:165428/71‑1 (MQ=255)
                             aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:5021681/71‑1 (MQ=255)
                             aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:4256620/71‑1 (MQ=255)
                             aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:3501760/71‑1 (MQ=255)
                              tGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTa                                                       >  1:3692779/1‑56 (MQ=255)
                              tGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGa                                                   <  1:1720756/60‑1 (MQ=255)
                              tGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGa                                                   <  1:4393245/60‑1 (MQ=255)
                               gTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGaaa                                                                  <  1:4776305/44‑1 (MQ=16)
                               gTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGaaa                                                                  <  1:2916993/44‑1 (MQ=16)
                               gTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                    <  1:3789488/42‑1 (MQ=16)
                               gTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                    <  1:4228840/42‑1 (MQ=16)
                               gTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:5561247/1‑44 (MQ=255)
                               gTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:4270924/1‑44 (MQ=255)
                               gTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:2373620/1‑44 (MQ=255)
                               gTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:1446979/1‑44 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                  >  1:4208479/1‑43 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAcc                                                              <  1:1531179/47‑1 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTg                                        >  1:3513224/1‑69 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:1783731/70‑1 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:4059263/70‑1 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                      >  1:5302943/1‑71 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                      >  1:3599301/1‑71 (MQ=255)
                                tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                      >  1:176459/1‑71 (MQ=255)
                                  gagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCt                                                             >  1:5110416/1‑46 (MQ=255)
                                   agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTc                                   >  1:565703/1‑71 (MQ=255)
                                    gTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                  <  1:3840340/71‑1 (MQ=255)
                                     ttCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGa                                 <  1:671639/71‑1 (MQ=255)
                                       ccATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGaa                                <  1:2076415/70‑1 (MQ=255)
                                       ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGaa                                                  >  1:1003256/1‑52 (MQ=255)
                                       ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       >  1:3729672/1‑63 (MQ=255)
                                       ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       >  1:2682370/1‑63 (MQ=255)
                                       ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAttt                                    >  1:2290000/1‑66 (MQ=255)
                                        catcatTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGtt                                               <  1:4010804/54‑1 (MQ=255)
                                        catcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:5653391/62‑1 (MQ=255)
                                        catcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:5546947/62‑1 (MQ=255)
                                        catcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:4511490/62‑1 (MQ=255)
                                        catcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGa                                       <  1:5484193/62‑1 (MQ=255)
                                          tcatTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGa                                       >  1:774495/1‑60 (MQ=255)
                                          tcatTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATacc                            >  1:2523856/1‑70 (MQ=255)
                                            atGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACg                          >  1:5339757/1‑71 (MQ=255)
                                              tgTGAACTCAGGATTATGGCGAACAGAAATACCTTc                                                           <  1:3865925/35‑1 (MQ=12)
                                              ggTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGa                                 <  1:2793424/62‑1 (MQ=255)
                                               gTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCgtt                        >  1:1080564/1‑68 (MQ=255)
                                                tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                      <  1:1146341/71‑1 (MQ=255)
                                                tGAACTCAGGATTATGTCGAACAGAAATACCTTCATTACgg                                                    <  1:480592/41‑1 (MQ=255)
                                                tGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTca                      <  1:554714/71‑2 (MQ=255)
                                                 gAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                      >  1:4582207/1‑70 (MQ=255)
                                                  aaCTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                         <  1:91075/50‑1 (MQ=255)
                                                  aaCTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                    <  1:4827354/71‑1 (MQ=255)
                                                  aaCTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                    <  1:4527674/71‑1 (MQ=255)
                                                  aaCTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTcaa                     >  1:3317304/1‑68 (MQ=255)
                                                    ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                  <  1:2892083/71‑1 (MQ=255)
                                                    ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                  >  1:4649937/1‑71 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACg                                            <  1:5368982/42‑1 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGa                     <  1:675965/65‑1 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccac               >  1:1408252/1‑71 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccac               >  1:5437985/1‑71 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccac               >  1:2777119/1‑71 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                  <  1:1602887/68‑1 (MQ=255)
                                                       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAAGGTTGATTTCGAATACACGCTCGAAGccac               >  1:2157608/1‑71 (MQ=255)
                                                       tGGGTTATCACGTACGGAAATACCTTCGTTACGGAAGTTACg                                            <  1:1613899/42‑1 (MQ=255)
                                                         ggTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAg                   >  1:4676412/1‑65 (MQ=255)
                                                         gattATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAAccgcc              >  1:4934209/3‑70 (MQ=255)
                                                          attATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCg                                  >  1:2875523/2‑49 (MQ=37)
                                                             aTGGCGAACAGAAATACCTTCATTACGGAAGTTACg                                            >  1:5447087/1‑36 (MQ=19)
                                                             aTGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaa                   >  1:277397/1‑60 (MQ=255)
                                                             aTGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCgcctacaa         >  1:748797/1‑69 (MQ=255)
                                                             aTGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCgcctaca          <  1:1084518/70‑2 (MQ=255)
                                                                aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac        >  1:1326860/1‑69 (MQ=255)
                                                                aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAAccacca      <  1:2818589/71‑1 (MQ=255)
                                                                 cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAg     <  1:1516433/71‑1 (MQ=255)
                                                                   tACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                      <  1:1002749/36‑1 (MQ=255)
                                                                    aCGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACg  >  1:1033978/1‑71 (MQ=255)
                                                                     cGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                    <  1:4410153/52‑1 (MQ=255)
                                                                      |                                                                    
CGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACG  >  minE/1933569‑1933707

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: