Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,252 A→C 10.1% V397V (GTT→GTG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2520AC10.1% 141.7 / 11.2 59V397V (GTT→GTGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base C (3/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATC  >  minE/1933184‑1933321
                                                                    |                                                                     
ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                                     <  1:2176779/71‑1 (MQ=255)
 caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                    <  1:3398231/71‑1 (MQ=255)
 caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                    <  1:188691/71‑1 (MQ=255)
  accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                                   >  1:1417476/1‑71 (MQ=255)
     aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTg                                                                >  1:3490247/1‑71 (MQ=255)
       tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                    <  1:517299/65‑1 (MQ=255)
         aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTc                                                                  >  1:821439/1‑65 (MQ=255)
         aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                                             <  1:3500476/70‑1 (MQ=255)
          agaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                           <  1:900427/71‑1 (MQ=255)
              cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                           <  1:3607278/67‑1 (MQ=255)
                 aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGata                                                         <  1:2608841/66‑1 (MQ=255)
                 aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAg                                                    >  1:592188/1‑71 (MQ=255)
                  aaacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGt                                                   <  1:2930847/68‑1 (MQ=255)
                      ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                               <  1:3276287/71‑1 (MQ=255)
                       gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                              <  1:804056/71‑1 (MQ=255)
                          tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                           <  1:1018364/55‑1 (MQ=255)
                          tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg                                           >  1:946646/1‑71 (MQ=255)
                            tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                            <  1:1782671/68‑1 (MQ=255)
                            tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                         >  1:3032730/1‑71 (MQ=255)
                            tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                         >  1:597976/1‑71 (MQ=255)
                            tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGacac                                                                     <  1:1363723/43‑1 (MQ=16)
                            tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGt                                          <  1:1036006/70‑1 (MQ=255)
                             gATTTCCGGTTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                                        >  1:2912796/1‑71 (MQ=255)
                              aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                          >  1:441733/1‑68 (MQ=255)
                              aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCgg                                       <  1:694078/71‑1 (MQ=255)
                              aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACtt                                                                   >  1:1314136/1‑43 (MQ=16)
                               tttCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                                       <  1:2840527/54‑1 (MQ=255)
                                 tCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCt                                     >  1:812172/1‑70 (MQ=255)
                                    gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                 <  1:548156/71‑1 (MQ=255)
                                        tAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                 >  1:2835017/1‑67 (MQ=255)
                                        tAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGat                                                          >  1:3854578/1‑42 (MQ=16)
                                         aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCa                                                     >  1:2515140/1‑46 (MQ=255)
                                         aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                            >  1:1902277/1‑71 (MQ=255)
                                           cggcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt                          <  1:595333/68‑1 (MQ=255)
                                             tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                           <  1:1600906/36‑1 (MQ=255)
                                             tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGt                                                   <  1:1739025/44‑1 (MQ=255)
                                             tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                          >  1:2519386/1‑53 (MQ=255)
                                               gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGaa                                             >  1:2960271/1‑48 (MQ=255)
                                               gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                               <  1:114563/62‑1 (MQ=255)
                                               gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                      <  1:3072780/71‑1 (MQ=255)
                                               gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                      >  1:3825365/1‑71 (MQ=255)
                                               attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTcc                                                     >  1:3453221/2‑39 (MQ=19)
                                                  acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCa                                                     >  1:194979/1‑37 (MQ=255)
                                                  acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgtg                         <  1:960515/65‑1 (MQ=255)
                                                  acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTg                   >  1:4096231/1‑71 (MQ=255)
                                                   cgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTc                     >  1:3843282/1‑68 (MQ=255)
                                                    gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                            <  1:2994512/44‑1 (MQ=255)
                                                    gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                            <  1:626717/44‑1 (MQ=255)
                                                     acgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                          <  1:1586248/45‑1 (MQ=255)
                                                      cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc               >  1:1345709/1‑71 (MQ=255)
                                                      cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc               >  1:2553201/1‑71 (MQ=255)
                                                           cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                              >  1:229937/1‑35 (MQ=255)
                                                             gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtg                           <  1:4143974/52‑1 (MQ=255)
                                                             gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg        <  1:3014788/71‑1 (MQ=255)
                                                               ggAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGa                             <  1:1030644/48‑1 (MQ=255)
                                                                gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                      >  1:265143/1‑54 (MQ=255)
                                                                gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                      >  1:1658654/1‑54 (MQ=255)
                                                                gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCGTCGAAg     <  1:536310/71‑1 (MQ=255)
                                                                  gaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg        >  1:3434060/1‑66 (MQ=255)
                                                                  gaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg        >  1:441372/1‑66 (MQ=255)
                                                                    aaCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtc  <  1:407845/70‑1 (MQ=255)
                                                                    |                                                                     
CCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATC  >  minE/1933184‑1933321

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: