Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,066,377 G→A 14.3% noncoding (67/76 nt) ileX → tRNA‑Ile

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,066,3770GA14.3% 110.5 / 9.6 42noncoding (67/76 nt)ileXtRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base G (12/24);  new base A (3/3);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTT  >  minE/2066307‑2066446
                                                                      |                                                                     
ccATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       <  1:790497/71‑1 (MQ=255)
ccATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       >  1:3847552/1‑71 (MQ=255)
   aggCCCTTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAaggg                                                                    <  1:1904831/70‑1 (MQ=255)
     gCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       >  1:1885658/1‑66 (MQ=255)
       cccTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGccac                                                                <  1:1674132/71‑1 (MQ=255)
         tttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGccacc                                                               >  1:2112059/2‑70 (MQ=255)
         cTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGccacca                                                              <  1:956522/71‑1 (MQ=255)
          ttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAggg                                                                     <  1:67619/63‑1 (MQ=255)
          ttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGccacca                                                              >  1:2498572/1‑70 (MQ=255)
           tAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAgg                                                                      >  1:3434697/1‑61 (MQ=255)
           tAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGccac                                                                >  1:2721694/1‑67 (MQ=255)
            aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAg                                                             <  1:1772515/69‑1 (MQ=255)
            aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAg                                                             <  1:4120451/69‑1 (MQ=255)
            aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAg                                                             <  1:3260932/69‑1 (MQ=255)
            aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAg                                                             <  1:3544315/69‑1 (MQ=255)
            aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:2062071/70‑1 (MQ=255)
              cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       >  1:586414/1‑57 (MQ=255)
                 aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGccaccaaa                                                            >  1:2625063/1‑63 (MQ=255)
                 aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGCCACCAAatcat                                                        <  1:2154691/69‑4 (MQ=255)
                 aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGCCACCAAatcat                                                        <  1:2449297/69‑4 (MQ=255)
                   tGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:59521/63‑1 (MQ=255)
                     gTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       >  1:1239018/1‑50 (MQ=255)
                     gTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAgg                                                                      >  1:1850312/1‑51 (MQ=255)
                     gTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:2744322/61‑1 (MQ=255)
                      ttAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCcagcag                                                                       <  1:1437258/49‑1 (MQ=255)
                      ttAGAGCAGGCGACTCATAATCGCTTGCTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:223860/60‑1 (MQ=255)
                       tAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGccac                                                                >  1:2936499/1‑55 (MQ=255)
                       tAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGc                                                >  1:2826610/1‑71 (MQ=255)
                          agCAGGCGACTCATAATCGCTTGGGCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:3356188/56‑1 (MQ=255)
                              ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCa                                                     <  1:2603009/59‑1 (MQ=255)
                                cGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaa                                       <  1:885083/71‑1 (MQ=255)
                                 gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaaa                                      >  1:4129122/1‑71 (MQ=255)
                                 gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaaa                                      >  1:3758909/1‑71 (MQ=255)
                                 gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaaa                                      <  1:3715632/71‑1 (MQ=255)
                                 gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaaa                                      >  1:418507/1‑71 (MQ=255)
                                   cTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATc                                    <  1:2256609/71‑1 (MQ=255)
                                    tCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGccacca                                                              <  1:1978814/44‑1 (MQ=255)
                                    tCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGa                                                            <  1:214694/46‑1 (MQ=255)
                                    tCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAAGGGCCACCAAatc                                                          <  1:2157524/48‑2 (MQ=38)
                                      aTAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGt                                  >  1:1697422/1‑70 (MQ=255)
                                       tAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCa                                                     <  1:3486864/50‑1 (MQ=255)
                                         aTCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGaa                                       >  1:350823/1‑62 (MQ=255)
                                         aTCGCTTGGGCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGc                              <  1:3808411/71‑1 (MQ=255)
                                            gCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctt                           >  1:3736302/1‑71 (MQ=255)
                                               tGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCcttttt                        >  1:3136172/1‑71 (MQ=255)
                                                  tCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCg                                   <  1:2520910/57‑1 (MQ=255)
                                                     cTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                  <  1:4048854/71‑1 (MQ=255)
                                                       ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                 >  1:2511947/1‑70 (MQ=255)
                                                        gTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                 <  1:1691565/69‑1 (MQ=255)
                                                        gTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                 <  1:722160/69‑1 (MQ=255)
                                                         ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGt                                  <  1:2363819/51‑1 (MQ=255)
                                                           cAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTtatata             >  1:2964346/1‑70 (MQ=255)
                                                                 cagcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                 >  1:629353/1‑60 (MQ=255)
                                                                   gcagGGGCCACCAGATATAGCAAAGGCTGACGAGaa                                       <  1:3347453/36‑1 (MQ=255)
                                                                   gcagGGGCCACCAGATATAGCAAAGGCTGACGAGaa                                       <  1:2825457/36‑1 (MQ=255)
                                                                   gcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCct                            <  1:680299/47‑1 (MQ=255)
                                                                    cagGGGCCACCAGATATAGCAAAGGCTGACGAGaa                                       <  1:4077706/35‑1 (MQ=255)
                                                                     agGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACttt  >  1:1208320/1‑71 (MQ=255)
                                                                      |                                                                     
CCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTT  >  minE/2066307‑2066446

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: