Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,131 G→C 7.9% Q438E (CAA→GAA)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1310GC 7.9% 186.3 / 5.8 63Q438E (CAA→GAA) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (34/24);  new base C (3/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTC  >  minE/1933062‑1933200
                                                                     |                                                                     
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                      >  1:1938148/1‑71 (MQ=255)
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                      >  1:2004059/1‑71 (MQ=255)
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                      >  1:3215066/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                     >  1:343493/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                     >  1:2775405/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                     >  1:1399802/1‑71 (MQ=255)
    ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                                                  >  1:1560276/1‑71 (MQ=255)
     tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                                                  >  1:3326727/1‑70 (MQ=255)
     tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                       >  1:1112987/1‑65 (MQ=255)
     tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                       >  1:1108282/1‑65 (MQ=255)
      cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGcc                                                                   <  1:1098252/68‑1 (MQ=255)
       aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                                               >  1:2998414/1‑71 (MQ=255)
        tcatAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                                                  <  1:3612344/64‑1 (MQ=255)
        tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                      >  1:3085872/1‑63 (MQ=255)
        tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                                       >  1:878957/1‑62 (MQ=255)
        tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                                              <  1:2753595/71‑1 (MQ=255)
         cGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGAGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCt                                                             >  1:2037954/1‑71 (MQ=255)
          atAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCtt                                                            <  1:1614705/70‑1 (MQ=255)
             gAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGAt                                                               >  1:507916/1‑65 (MQ=255)
                cATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                                                      <  1:1095778/71‑1 (MQ=255)
                    gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAtt                                                  >  1:133817/1‑71 (MQ=255)
                    gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAt                                                   <  1:4038586/70‑1 (MQ=255)
                      ttCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTc                                                           >  1:511336/1‑60 (MQ=255)
                         gtcgtcACCTGCGTCTTTCGCGGCAACATGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt                                             >  1:2566904/1‑71 (MQ=255)
                          tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTc                                            >  1:1773734/1‑71 (MQ=255)
                           cgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                                                      <  1:1239521/60‑1 (MQ=255)
                            gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAg                                               >  1:3465021/1‑66 (MQ=255)
                            gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGc                                          >  1:2524913/1‑71 (MQ=255)
                              cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                                                  >  1:4053302/1‑45 (MQ=255)
                              cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCt                                         >  1:113009/1‑70 (MQ=255)
                               aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTaa                                       <  1:144523/71‑1 (MQ=255)
                               aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                        >  1:578901/1‑70 (MQ=255)
                               aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                        <  1:3380628/70‑1 (MQ=255)
                                 cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTg                                                                 <  1:1579981/43‑1 (MQ=255)
                                 cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                                     >  1:3917514/1‑71 (MQ=255)
                                 cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                                     >  1:3774672/1‑71 (MQ=255)
                                  tGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc                                    >  1:4093441/1‑71 (MQ=255)
                                   gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   >  1:1709237/1‑71 (MQ=255)
                                   gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   >  1:2384922/1‑71 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   <  1:3474159/70‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   <  1:295992/70‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   <  1:1423393/70‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc                                    <  1:2446426/69‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCg                                  >  1:2379679/1‑71 (MQ=255)
                                      tCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGt                                 >  1:2540938/1‑70 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:3077415/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:1606958/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:1407380/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:3740954/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:2465105/1‑71 (MQ=255)
                                         ttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAg                                               <  1:3769836/53‑1 (MQ=255)
                                         ttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                        <  1:1485610/60‑1 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTaga                             >  1:970288/1‑68 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                             <  1:3442265/70‑1 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                             >  1:3708823/1‑70 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                             >  1:1866890/1‑70 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                             <  1:664871/70‑1 (MQ=255)
                                               gCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                   <  1:1474792/59‑1 (MQ=255)
                                                 aaCCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCAc                     <  1:1239501/71‑1 (MQ=255)
                                                      gATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt                                             >  1:1557447/1‑42 (MQ=255)
                                                       aTCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGc                                              >  1:3046573/1‑40 (MQ=255)
                                                        tCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGaccacc              <  1:1941729/71‑1 (MQ=255)
                                                            ggAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAAt           >  1:858607/1‑70 (MQ=255)
                                                             gAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGa                   >  1:942031/1‑61 (MQ=255)
                                                               aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAGAGCCGTTACCGATTTCACGACCACCAATgaag       <  1:1008066/71‑1 (MQ=255)
                                                               aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATtaag       <  1:949388/71‑1 (MQ=255)
                                                               aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATgaag       <  1:2645024/71‑1 (MQ=255)
                                                               aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATg          >  1:1155013/1‑68 (MQ=255)
                                                               aGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATg          >  1:699282/1‑68 (MQ=255)
                                                                  gTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:856508/1‑45 (MQ=255)
                                                                  gTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                              >  1:1796153/1‑45 (MQ=255)
                                                                   ttGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                             <  1:731182/45‑1 (MQ=255)
                                                                   ttGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCAc                     <  1:3467781/53‑1 (MQ=255)
                                                                     gcgcCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGa                           >  1:1170541/1‑45 (MQ=255)
                                                                     gcgcCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTc  <  1:226877/70‑1 (MQ=255)
                                                                     |                                                                     
AATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTC  >  minE/1933062‑1933200

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