Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,044 G→A 12.1% H467Y (CAT→TAT)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0440GA12.1% 593.3 / 52.2 223H467Y (CAT→TAT) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (112/84);  new base A (15/12);  total (127/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 223 total aligned reads displayed.

GCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAA  >  minE/1932976‑1933112
                                                                    |                                                                    
gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGtt                                                                    >  1:259517/1‑71 (MQ=255)
gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGtt                                                                    >  1:4266918/1‑71 (MQ=255)
 cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc                                                                   <  1:2287609/71‑1 (MQ=255)
    ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 <  1:2188498/70‑1 (MQ=255)
      ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg                                                              >  1:409699/1‑71 (MQ=255)
       gTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg                                                              <  1:4421421/70‑1 (MQ=255)
        tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                            <  1:594551/71‑1 (MQ=255)
        tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                            >  1:837025/1‑71 (MQ=255)
         aaaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAgcgcag                                                           >  1:2897954/2‑69 (MQ=255)
          aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt                                                          <  1:4710148/71‑1 (MQ=255)
          aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt                                                          <  1:244594/71‑1 (MQ=255)
          aaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAgcgca                                                            >  1:3310219/1‑68 (MQ=12)
          aaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGt                                                          >  1:1117808/1‑71 (MQ=255)
             aGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 >  1:4575428/1‑61 (MQ=255)
              gCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                            >  1:1486122/1‑65 (MQ=255)
                  tACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTaa                                                   >  1:1038074/1‑70 (MQ=255)
                   aCCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 >  1:4662066/1‑55 (MQ=255)
                   aCCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 >  1:106668/1‑55 (MQ=255)
                    ccATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                 <  1:4487654/70‑1 (MQ=255)
                    ccATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                 <  1:2725785/70‑1 (MQ=255)
                      aTACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCtt                                               >  1:4060993/1‑70 (MQ=255)
                      aTACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCtt                                               >  1:4324786/1‑70 (MQ=255)
                       tACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCa                                             >  1:4403023/1‑71 (MQ=255)
                       tACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCa                                             >  1:3330490/1‑71 (MQ=255)
                         cGGTCGATACCAATTCCCAGACCTGCTGTCGTCGGTAAGCCATGTTCCa                                                                 >  1:1506123/1‑49 (MQ=255)
                         cGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATc                                           >  1:2249762/1‑71 (MQ=255)
                         cGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTcgt                                            <  1:1316448/70‑3 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc                                                                   <  1:4285014/44‑1 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:2170557/1‑69 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:2556433/1‑69 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:5486413/1‑69 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:3344681/1‑69 (MQ=255)
                            tCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:3678290/1‑69 (MQ=255)
                            tCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATa                                        <  1:310847/71‑1 (MQ=255)
                             cGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                          >  1:4360687/1‑68 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTa                                        >  1:177889/1‑69 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                       >  1:715726/1‑70 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                       >  1:3623786/1‑70 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                       >  1:344547/1‑70 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                       >  1:364375/1‑70 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                       <  1:4864166/70‑1 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTAGGCGGTAAGCCATGTTCCa                                                                 >  1:5280513/1‑44 (MQ=255)
                               aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                     <  1:2935772/71‑1 (MQ=255)
                                tACCAATTCCCAGACCTGCTGTCGGCGGTAAGGCATGTTCCAGTGCGGTGa                                                        >  1:3754797/1‑51 (MQ=255)
                                tACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc                                                                   <  1:3074998/40‑1 (MQ=255)
                                tACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGa                                                        >  1:1707581/1‑51 (MQ=255)
                                tACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTc                                              >  1:4768398/1‑61 (MQ=255)
                                 aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCa                                             >  1:276887/1‑61 (MQ=255)
                                 aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCAt                                            >  1:2142107/1‑62 (MQ=255)
                                 aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                   >  1:3031612/1‑71 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGa                                      <  1:203293/67‑1 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                   <  1:4927099/70‑1 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                   <  1:5325187/70‑1 (MQ=255)
                                  ccAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACa                                   <  1:2339575/70‑1 (MQ=255)
                                   cAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                 >  1:5193060/1‑71 (MQ=255)
                                   cAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                 >  1:5018567/1‑71 (MQ=255)
                                    aaTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 >  1:1326279/1‑38 (MQ=255)
                                    aaTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                                                                 >  1:876057/1‑38 (MQ=255)
                                      ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                 <  1:203776/68‑1 (MQ=255)
                                       tcacAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCtt                             >  1:227026/4‑71 (MQ=255)
                                       tCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCgg                                                           >  1:2886909/1‑41 (MQ=255)
                                       tCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCtt                             >  1:2434104/1‑71 (MQ=255)
                                        cccAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtc                            <  1:5442157/71‑1 (MQ=255)
                                         ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcg                           <  1:2704678/71‑1 (MQ=255)
                                         ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                            <  1:1557915/70‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                 <  1:5576144/48‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTc                                              <  1:3379930/51‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTc                                              <  1:3304704/51‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCAt                                            >  1:679348/1‑53 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCAt                                            >  1:4010927/1‑53 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                          >  1:156595/1‑71 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCAAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                          <  1:1257223/71‑1 (MQ=255)
                                           aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                         <  1:2018189/71‑1 (MQ=255)
                                           aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                         <  1:2215396/71‑1 (MQ=255)
                                           aGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtcgtc                         <  1:642669/71‑1 (MQ=255)
                                           aGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtcgtc                         >  1:3930457/1‑71 (MQ=255)
                                             aCCTTCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                       <  1:5632619/71‑1 (MQ=255)
                                             aCCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTa                                                    <  1:1659818/42‑1 (MQ=255)
                                               ctgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCa                     <  1:3368009/71‑1 (MQ=255)
                                                tgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                    >  1:3667897/1‑71 (MQ=255)
                                                  ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCt                  >  1:1990703/1‑71 (MQ=255)
                                                    gTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                         >  1:4888566/1‑62 (MQ=255)
                                                    gTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                         >  1:4159028/1‑62 (MQ=255)
                                                     tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTg                 <  1:282302/69‑1 (MQ=255)
                                                        gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt            <  1:132381/71‑1 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc         <  1:4807766/71‑1 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc         <  1:3262295/71‑1 (MQ=255)
                                                            tAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                       <  1:863496/56‑1 (MQ=255)
                                                            tAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                       <  1:141426/56‑1 (MQ=255)
                                                            tAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg        <  1:5427428/71‑1 (MQ=255)
                                                             aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                       >  1:2502760/1‑55 (MQ=255)
                                                              aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                          >  1:5667727/1‑51 (MQ=255)
                                                              aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAc                    >  1:5115191/1‑57 (MQ=255)
                                                                ccATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTtc                            <  1:1729618/47‑1 (MQ=255)
                                                                ccATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagcatt    <  1:2202938/71‑4 (MQ=255)
                                                                 cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc    >  1:2967475/1‑70 (MQ=255)
                                                                 cATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACa                                   >  1:1107454/1‑39 (MQ=12)
                                                                 cATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACa                                   >  1:3964936/1‑39 (MQ=12)
                                                                  aTGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCaa  >  1:1084495/1‑71 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg      <  1:3178810/66‑1 (MQ=255)
                                                                    |                                                                    
GCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAA  >  minE/1932976‑1933112

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: