Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,952 C→G 10.0% L497L (CTG→CTC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9520CG10.0% 226.9 / 12.8 80L497L (CTG→CTClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (42/30);  new base G (2/6);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGA  >  minE/1932883‑1933021
                                                                     |                                                                     
aaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACa                                                                      <  1:3511213/71‑1 (MQ=255)
  aaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAg                                                                     >  1:1444503/1‑70 (MQ=255)
   aCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGaa                                                                   >  1:3369038/1‑71 (MQ=255)
   aCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGa                                                                    >  1:3018266/1‑70 (MQ=255)
     cccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATa                                                                 <  1:3796069/71‑1 (MQ=255)
      ccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATaa                                                                >  1:1444569/1‑71 (MQ=255)
        gACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACg                                                              >  1:3567785/1‑71 (MQ=255)
          cATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTc                                                            >  1:2945439/1‑71 (MQ=255)
           aTTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACa                                                                      <  1:173788/60‑1 (MQ=255)
             ttGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcgcg                                                         >  1:3018044/1‑71 (MQ=255)
               gCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACa                                                                      <  1:3562910/56‑1 (MQ=255)
                 cGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg                                                     >  1:2432956/1‑71 (MQ=255)
                 cGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg                                                     >  1:1225247/1‑71 (MQ=255)
                  ggggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGa                                                       <  1:287509/68‑1 (MQ=255)
                  ggggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGa                                                       <  1:2862970/68‑1 (MQ=255)
                   gggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGt                                                   <  1:41688/71‑1 (MQ=255)
                    ggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTa                                                  >  1:3629386/1‑71 (MQ=255)
                    ggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTa                                                  >  1:1318654/1‑71 (MQ=255)
                     gTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                 <  1:2428158/71‑1 (MQ=255)
                         tgAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGc                                              <  1:540989/70‑1 (MQ=255)
                          gAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTg                                            >  1:2583037/1‑71 (MQ=255)
                            gCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAAt                                                                  <  1:1629071/47‑1 (MQ=255)
                            gCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATgg                                               >  1:1979411/1‑66 (MQ=255)
                            gCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATgg                                               >  1:950560/1‑66 (MQ=255)
                               tAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcg                                                           <  1:155957/51‑1 (MQ=255)
                               tAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcg                                                           <  1:286593/51‑1 (MQ=255)
                                aaCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGt                                       >  1:2621741/1‑70 (MQ=255)
                                 aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGaa                                                                   >  1:1927194/1‑41 (MQ=255)
                                 aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGaa                                                                   >  1:589158/1‑41 (MQ=255)
                                 aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                      <  1:2531095/70‑1 (MQ=255)
                                  cGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                    <  1:2635929/71‑1 (MQ=255)
                                    tAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGt                                           >  1:3772398/1‑62 (MQ=255)
                                     aaTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATg                                                >  1:653092/1‑56 (MQ=255)
                                     aaTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                 >  1:645167/1‑71 (MQ=255)
                                      aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                 >  1:3127157/1‑70 (MQ=255)
                                      aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                 >  1:3537682/1‑70 (MQ=255)
                                          ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                            >  1:1149017/1‑71 (MQ=255)
                                           tATTTTACCGGACGCATCGCCGGGAACAGAATAACGt                                                             <  1:1271524/37‑1 (MQ=255)
                                             ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                          <  1:2293919/54‑1 (MQ=255)
                                             ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAcc                          <  1:1169249/70‑1 (MQ=255)
                                             ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                         >  1:878366/1‑71 (MQ=255)
                                             ttttACCGGACGCATCGCCGCGAACAGCATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                         >  1:2314682/1‑71 (MQ=255)
                                               ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                 <  1:2323056/45‑1 (MQ=255)
                                               ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                 <  1:2990411/45‑1 (MQ=255)
                                               ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                              >  1:3650037/1‑64 (MQ=255)
                                                tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                       <  1:1280495/70‑1 (MQ=255)
                                                tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                      >  1:2324401/1‑71 (MQ=255)
                                                  ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCGGTTGGTGAACAGCATTACCATACg                     <  1:3049851/70‑1 (MQ=255)
                                                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                          >  1:1446317/1‑47 (MQ=255)
                                                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                  <  1:1995375/71‑1 (MQ=255)
                                                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                  <  1:142934/71‑1 (MQ=255)
                                                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                  <  1:3554281/71‑1 (MQ=255)
                                                     gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                  >  1:2499762/1‑70 (MQ=255)
                                                        gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                              >  1:3378994/1‑55 (MQ=255)
                                                        gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                              >  1:1704730/1‑55 (MQ=255)
                                                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                         <  1:255271/43‑1 (MQ=255)
                                                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                         <  1:2212840/43‑1 (MQ=255)
                                                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                         <  1:2511229/43‑1 (MQ=255)
                                                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGt                                           >  1:3600651/1‑40 (MQ=255)
                                                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                         >  1:342351/1‑58 (MQ=255)
                                                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc             >  1:399179/1‑70 (MQ=255)
                                                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc             >  1:1589153/1‑70 (MQ=255)
                                                           tCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                          <  1:998539/40‑1 (MQ=255)
                                                            cGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                        >  1:1545448/1‑41 (MQ=255)
                                                             gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                    >  1:2507471/1‑44 (MQ=21)
                                                             gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa          >  1:389074/1‑70 (MQ=255)
                                                              ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATa              <  1:2249842/65‑1 (MQ=255)
                                                              ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta        <  1:1387639/71‑2 (MQ=255)
                                                              ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt         >  1:2309855/1‑70 (MQ=255)
                                                               cGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta        <  1:3545116/70‑2 (MQ=255)
                                                                gggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc      <  1:2036599/71‑1 (MQ=255)
                                                                 ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                        <  1:2573551/52‑1 (MQ=255)
                                                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTTGTGat                                    <  1:113343/40‑2 (MQ=39)
                                                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGa                <  1:957812/60‑1 (MQ=255)
                                                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc      >  1:2420740/1‑70 (MQ=255)
                                                                  gAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa    <  1:3836490/71‑1 (MQ=255)
                                                                  gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc      >  1:51047/1‑69 (MQ=255)
                                                                  gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc      >  1:446491/1‑69 (MQ=255)
                                                                   aaCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc      >  1:2017540/1‑68 (MQ=255)
                                                                    aCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGa  >  1:2874504/1‑71 (MQ=255)
                                                                     |                                                                     
AAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGA  >  minE/1932883‑1933021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: