Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,136,857 A→T 5.1% intergenic (‑80/+69) rpsO ← / ← truB 30S ribosomal subunit protein S15/tRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,136,8570AT 5.1% 288.0 / 5.7 98intergenic (‑80/+69)rpsO/truB30S ribosomal subunit protein S15/tRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base A (56/36);  new base T (3/2);  total (59/39)
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.72e-01

Only 100 of 106 total aligned reads displayed.

CAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAT  >  minE/2136788‑2136926
                                                                     |                                                                     
cAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTa                                                                       >  1:310945/1‑70 (MQ=255)
  aaGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAAc                                                                    <  1:1924985/71‑1 (MQ=255)
     tatataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAt                                                                 <  1:3145547/71‑1 (MQ=255)
     tatataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAt                                                                 <  1:551907/71‑1 (MQ=255)
      atataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                                >  1:571201/1‑71 (MQ=255)
       tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAttt                                                               <  1:3568197/71‑1 (MQ=255)
       tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                                <  1:106420/70‑1 (MQ=255)
       tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                                <  1:940605/70‑1 (MQ=255)
          aGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAc                                                             >  1:1300620/1‑70 (MQ=255)
           gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACTTTGcgt                                                           <  1:2934059/71‑2 (MQ=37)
           gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACTTTGcgt                                                           <  1:3008278/71‑2 (MQ=37)
           gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcgc                                                           <  1:1493192/71‑1 (MQ=255)
           gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcg                                                            >  1:2141521/1‑70 (MQ=255)
           gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcg                                                            <  1:846311/70‑1 (MQ=255)
              tGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGAc                                                        >  1:3107607/1‑71 (MQ=255)
                aaaGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcg                                                            >  1:548241/1‑65 (MQ=255)
                aaaGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGt                                                      >  1:1814420/1‑71 (MQ=255)
                 aaGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGt                                                      >  1:1774086/1‑70 (MQ=255)
                 aaGGACGACGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAAGCGACGt                                                      >  1:1913191/1‑70 (MQ=255)
                  aGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:1065076/71‑1 (MQ=255)
                   ggACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGtt                                                     >  1:1164704/1‑69 (MQ=255)
                    gACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcgc                                                           >  1:657847/1‑62 (MQ=255)
                    gACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAcgc                                                           >  1:950517/1‑62 (MQ=255)
                    gACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGAc                                                        <  1:3672877/65‑1 (MQ=255)
                      ctccGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACaa                                                                  <  1:839327/51‑1 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:1358636/1‑70 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:3425046/1‑70 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:1365212/1‑70 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:824382/1‑70 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:734205/1‑70 (MQ=255)
                           aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:700327/1‑70 (MQ=255)
                             ctctAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:649635/1‑68 (MQ=255)
                             ctctAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:649637/1‑68 (MQ=255)
                               ctAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGa                                                         >  1:112193/1‑53 (MQ=255)
                               ctAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAACAATTTACGCGACGTTAAGCGGCa                                             <  1:3719325/64‑1 (MQ=255)
                                tAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTAt                                         <  1:3480607/68‑1 (MQ=255)
                                 aaTTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACttt                                                               >  1:2084192/1‑45 (MQ=21)
                                 aaTTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACttt                                                               >  1:3032617/1‑45 (MQ=21)
                                 aaTTCAGAGATCCCAGTGTACGTTACGTAAATTGTTTAACACttt                                                               >  1:3027260/1‑45 (MQ=18)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:1744635/55‑1 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:2498458/55‑1 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:959648/55‑1 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:3653967/55‑1 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            >  1:1071269/1‑63 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTa                                          >  1:2389565/1‑65 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTc                                       <  1:930057/68‑1 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCt                                      >  1:2824995/1‑69 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTa                                     >  1:1741347/1‑70 (MQ=255)
                                  aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTAc                                    >  1:3707561/1‑71 (MQ=255)
                                     cAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTc                                  <  1:1956606/70‑1 (MQ=255)
                                      cgCGATACCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGa                                                         >  1:474675/2‑46 (MQ=255)
                                      aGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            <  1:1026229/59‑1 (MQ=255)
                                      aGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGt                                <  1:3603841/71‑1 (MQ=255)
                                       gCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTa                                          >  1:1798180/1‑60 (MQ=255)
                                       gCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTa                                          >  1:1230292/1‑60 (MQ=255)
                                         gATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            <  1:223756/56‑1 (MQ=255)
                                          aTCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                                    <  1:335425/47‑1 (MQ=255)
                                            cccAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAg                                                  >  1:1191493/1‑47 (MQ=255)
                                            cccAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAg                                                  >  1:3623713/1‑47 (MQ=255)
                                            cccAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGAc                          >  1:1385670/1‑71 (MQ=255)
                                            ccaaGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGAc                          >  1:2125968/4‑71 (MQ=255)
                                             ccAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACTTTAAGc                                                 <  1:482163/47‑1 (MQ=255)
                                              cAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                            <  1:760908/51‑1 (MQ=255)
                                                gtgtACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGt                      >  1:3019671/1‑71 (MQ=255)
                                                gtgtACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGt                      >  1:2189302/1‑71 (MQ=255)
                                                   tACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTAt                   >  1:2943063/1‑71 (MQ=255)
                                                    aCGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATc                  >  1:1952580/1‑71 (MQ=255)
                                                     cGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCg                 >  1:2446745/1‑71 (MQ=255)
                                                      gTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCa                                             <  1:1095974/42‑1 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                           >  1:1724092/1‑43 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                           >  1:2409155/1‑43 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                           >  1:685018/1‑43 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                           >  1:2885871/1‑43 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                           >  1:598693/1‑43 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGc                >  1:3712383/1‑70 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCa               <  1:395454/71‑1 (MQ=255)
                                                       ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGAGGCAGt                                           >  1:2101742/1‑43 (MQ=255)
                                                        tACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTa                               >  1:3399040/1‑54 (MQ=255)
                                                         actcAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGc                >  1:3602217/4‑68 (MQ=255)
                                                          cGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCt                        >  1:242419/1‑59 (MQ=255)
                                                          cGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCt                        >  1:2994281/1‑59 (MQ=255)
                                                           gCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGAc           <  1:388077/71‑1 (MQ=255)
                                                           gCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGAc           >  1:225774/1‑71 (MQ=255)
                                                             aaaGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGa            <  1:2341763/68‑1 (MQ=255)
                                                              taatgtTAACCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCa               <  1:2503040/61‑1 (MQ=255)
                                                               aGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACg          >  1:2273553/1‑68 (MQ=255)
                                                               aGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGtt       <  1:1024666/71‑1 (MQ=255)
                                                               aGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGtt       <  1:2375306/71‑1 (MQ=255)
                                                                gtgtTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCATGACGGTTa      <  1:1245749/71‑1 (MQ=255)
                                                                gtgtTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTa      <  1:1276068/71‑1 (MQ=255)
                                                                gtgtTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTa      <  1:1965784/71‑1 (MQ=255)
                                                                  gtTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCa               <  1:1858545/60‑1 (MQ=255)
                                                                   ttAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACa   >  1:1535568/1‑71 (MQ=255)
                                                                   ttAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACa   >  1:2959410/1‑71 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGAc                          >  1:2895039/1‑46 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAt  <  1:1761972/70‑1 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAt  >  1:460365/1‑70 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAt  >  1:830768/1‑70 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAt  <  1:2057907/70‑1 (MQ=255)
                                                                     aaaCAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAt  >  1:2752866/1‑70 (MQ=255)
                                                                     |                                                                     
CAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGACGGTTAACAT  >  minE/2136788‑2136926

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: