Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,136,848 C→T 5.1% intergenic (‑71/+78) 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,8480CT 5.1% 328.6 / 5.0 97intergenic (‑71/+78)rpsO/truB30S ribosomal subunit protein S15/tRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base C (59/33);  new base T (3/2);  total (62/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

TAAAACTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGAC  >  minE/2136780‑2136917
                                                                    |                                                                     
tAAAACTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCa                                                                      >  1:1554637/1‑70 (MQ=255)
tAAAACTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCa                                                                      >  1:1190510/1‑70 (MQ=255)
 aaaaCTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCa                                                                      <  1:1576504/69‑1 (MQ=255)
  aaaCTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCaaa                                                                    <  1:161270/70‑1 (MQ=255)
     cTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAgtgt                                                                >  1:1607663/1‑71 (MQ=255)
      tCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAgtgt                                                                >  1:213391/1‑70 (MQ=255)
       ccAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGtt                                                               >  1:1106301/1‑70 (MQ=255)
        cAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGtt                                                               <  1:2405003/69‑1 (MQ=255)
        cAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGtt                                                               >  1:1403445/1‑69 (MQ=255)
        cAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTa                                                              >  1:310945/1‑70 (MQ=255)
         aaaGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGtt                                                               <  1:3702071/68‑1 (MQ=255)
          aaGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAAc                                                           <  1:1924985/71‑1 (MQ=255)
             tatataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAt                                                        <  1:551907/71‑1 (MQ=255)
             tatataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAt                                                        <  1:3145547/71‑1 (MQ=255)
              atataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                       >  1:571201/1‑71 (MQ=255)
               tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAttt                                                      <  1:3568197/71‑1 (MQ=255)
               tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                       <  1:940605/70‑1 (MQ=255)
               tataGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAAtt                                                       <  1:106420/70‑1 (MQ=255)
                  aGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTAc                                                    >  1:1300620/1‑70 (MQ=255)
                   gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACTTTGcgt                                                  <  1:3008278/71‑2 (MQ=37)
                   gAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGTAAATTGTTTAACACTTTGcgt                                                  <  1:2934059/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:2866612/1‑62 (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:3425046/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:700327/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:1358636/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:1365212/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:824382/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:57251/1‑70 (MQ=255)
                                   aTCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:734205/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:3653967/55‑1 (MQ=255)
                                          aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                           <  1:959648/55‑1 (MQ=255)
                                          aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                           <  1:2498458/55‑1 (MQ=255)
                                          aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTa                                           <  1:1744635/55‑1 (MQ=255)
                                          aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:1071269/1‑63 (MQ=255)
                                          aTTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAg                                   >  1:3484737/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:2104992/70‑1 (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:151853/1‑60 (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:3623713/1‑47 (MQ=255)
                                                    cccAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAg                                         >  1:1191493/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:598693/1‑43 (MQ=255)
                                                               ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                  >  1:685018/1‑43 (MQ=255)
                                                               ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                  >  1:1724092/1‑43 (MQ=255)
                                                               ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                  >  1:2885871/1‑43 (MQ=255)
                                                               ttACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGt                                  >  1:2409155/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:225774/1‑71 (MQ=255)
                                                                   gCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGAc  <  1:388077/71‑1 (MQ=255)
                                                                    |                                                                     
TAAAACTCCAAAGTATATAGAATGAAAGGACGCCGATCTCTAATTCAGCGATCCCAGTGTACGTTACGCAAAGTGTTAAACAATTTACGCGACGTTAAGCGGCAGTATTCTACTCGTAGCGACCTGTTATCGCAAGAC  >  minE/2136780‑2136917

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: