Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,132 C→A 7.9% A437A (GCG→GCT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1320CA 7.9% 211.9 / 3.0 63A437A (GCG→GCTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base A (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.97e-01

AATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCA  >  minE/1933062‑1933201
                                                                      |                                                                     
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                       >  1:1938148/1‑71 (MQ=255)
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                       >  1:2004059/1‑71 (MQ=255)
aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                                       >  1:3215066/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                      >  1:2775405/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                      >  1:1399802/1‑71 (MQ=255)
 aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                                      >  1:343493/1‑71 (MQ=255)
    ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                                                   >  1:1560276/1‑71 (MQ=255)
     tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                                                   >  1:3326727/1‑70 (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)
        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:3380628/70‑1 (MQ=255)
                               aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                                         >  1:578901/1‑70 (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:2384922/1‑71 (MQ=255)
                                   gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                    >  1:1709237/1‑71 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                    <  1:3474159/70‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                    <  1:1423393/70‑1 (MQ=255)
                                    cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                                    <  1:295992/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:3740954/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                               >  1:3077415/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                               >  1:2465105/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                               >  1:1407380/1‑71 (MQ=255)
                                        tttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                               >  1:1606958/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:664871/70‑1 (MQ=255)
                                          tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                              >  1:1866890/1‑70 (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:1796153/1‑45 (MQ=255)
                                                                  gTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                               >  1:856508/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)
                                                                      cgcCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCa  <  1:1674580/70‑1 (MQ=255)
                                                                      |                                                                     
AATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCA  >  minE/1933062‑1933201

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: